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I. INTRODUCTION 

The Functional Renormalisation Group (FRG) in its 
continuum formulation [1-14] has proven itself as a 
powerful tool for studying both perturbative and non- 
perturbative effects in quantum field theory and statis¬ 
tical physics, for reviews see [15-23]. In this approach a 
regularisation of a quantum theory is achieved by sup¬ 
pressing part of the propagating degrees of freedom re¬ 
lated to a cut-off scale k. This results in regularised gen¬ 
erating functionals such as the effective action F^ where 
part of the modes have been integrated out. The flow 
equation describes the response of the generating func¬ 
tional to an infinitesimal variation of k, and can be used 
to successively integrate-out modes. Hence, a generating 
functional at some initial scale A together with its flow 
serve as a definition of the quantum theory. For example, 
the flow equation allows us to calculate the full effective 
action F from an initial effective action Fa if the latter 
is well under control. For an infrared momentum cut-off 
and sufficiently large A we have a good grip on Fa as it 
can be computed perturbatively. 

The main advantages of such a formulation are its flex¬ 
ibility when it comes to truncations of the full theory, as 
well as its numerical accessibility. Both properties origi¬ 
nate in the same structural aspects of such flows. Quite 
generally functional flows are differential equations that 
relate an infinitesimal /c-variation of a generating func¬ 
tional Z with some functional of Z, its derivatives and the 
regulator. The quantum theory, and hence the physics, is 
solely specified by the boundary condition of such a flow. 
Due to this structure truncations are introduced on the 
level of the generating functional itself which leads to 
self-consistent truncated flows. Moreover, a change of 
degrees of freedom also is done on the level of the gener¬ 
ating functional, and the structure of the flow stays the 
same. Last but not least, numerical stability of the flow 
for a given problem and truncation is governed by the 
choice of the specific regularisation procedure. 

In other words, the advantages are carried by the struc¬ 
tural aspects of the functional RG, whose understand¬ 
ing and further development is the main purpose of the 
present work. It is not meant as a review and for a more 
complete list of references we refer the reader to the re¬ 
views already cited above, [15-23]. We close the intro¬ 
duction with an overview over the work. 

In section II we evaluate functional equations of quan¬ 
tum field theories, such as Dyson-Schwinger equations. 


symmetry identities, such as Slavnov-Taylor identities 
(STIs), and introduce some notation. 

In section III flows are derived for general correlation 
functions including those for the effective action and the 
Schwinger functional. We present a derivation of the flow 
equation which emphasises the subtleties of renormalisa¬ 
tion. Moreover, no use of the path integral representation 
is made, the derivation solely relies on the existence of a 
finite effective action or Schwinger functional for the full 
theory. First we introduce the setting and notion of regu¬ 
larisation. This is used to derive the general flows (3.28) 
and (3.60) which comprise the main results of this part. 
The flows discussed here include those for A^-particle ir¬ 
reducible (7VPI) quantities as well as relations between 
the different formulations. For general flows one has to 
carefully study the boundary conditions. A comparison 
of results obtained for different regularisations, in partic¬ 
ular in view of optimisation, requires the study of varia¬ 
tions of the regulator. 

In section IV we discuss the fate of RG equations of 
the full theory displaying reparameterisation invariance 
in the presence of a general regularisation. This is impor¬ 
tant when matching the scale dependence of quantities in 
the presence of the regularisation to that in the full the¬ 
ory without cut-off. The key RG flows are (4.8),(4.20) 
and are basically generalisations of (3.28) and (3.60). 

The important aspect of optimisation is investigated 
in section V. In most situations one has to rely on trun¬ 
cations to the full theory. Optimised flows should lead 
to results as close as possible to the full theory within 
each order of a given systematic truncation scheme. We 
develop a functional approach to optimisation of general 
flows which allows us to systematically access and de¬ 
velop optimisation criteria. We discuss the relation be¬ 
tween different optimisation ideas used in the literature. 
The definition of an effective cut-off scale is introduced 
and a constructive optimisation criterion is put forward 
in section VD. Roughly speaking, optimal regulators are 
those, that lead to correlation functions as close as pos¬ 
sible to that in the full theory for a given effective cut-off 
scale. 

The rest of the present paper deals with structural ap¬ 
plications of these findings. In section VI we relate flows 
to other functional methods such as Dyson-Schwinger 
equations or the use of NPl effective actions. To that 
end we consider flows in the presence of composite op¬ 
erators. In particular we construct practical renormal¬ 
isation schemes, the latter being of importance for the 
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renormalisation of Dyson-Schwinger equations and A^PI 
effective actions. 

A main motivation for the development of the present 
approach resides in its application to gauge theories. In 
section VII various structural aspects of gauge theories 
are investigated. We discuss the formulation of gauge 
theories using appropriate degrees of freedom. The mod¬ 
ification of symmetry identities in the presence of the 
regularisation and their different representations are eval¬ 
uated. The latter allow for a purely algebraic represen¬ 
tation of the symmetry identities. We also outline the 
construction of gauge-invariant flows and discuss the fate 
of gauge symmetry constraints in these formulations. We 
close with a brief evaluation of anomalous symmetries in 
the presence of a regulator. 

In section VIII we discuss consequences of the func¬ 
tional optimisation criterion and the RG equations for 
the construction of truncation schemes and optimal reg¬ 
ulators. It is shown that a specific class of regulators 
preserves the RG scalings of the underlying theory. We 
discuss the use of integrated flows that constitute finite 
renormalised Dyson-Schwinger equations. These inte¬ 
grated flows can be used in asymptotic regimes or a fixed 
point analysis within the functional RG setting. The con¬ 
structive optimisation criterion developed in section V is 
put to work within a simple example. Further applica¬ 
tions are outlined. 


II. PRELIMINARIES 

We consider the finite renormalised Euclidean Schwin¬ 
ger functional W[J] of the theory under investigation, 
where we do not only allow for source terms for the fun¬ 
damental fields (p of the theory, but also for sources for 
general tensorial composite operators flfp) with 

/ ^max 

dp\0\ exp|-S'[v>] + X! / 

n=l 

jai-a„(^l, ...,a:„)| . (2.1) 

Here at comprises possible Lorentz and gauge group in¬ 
dices and species of fields. The measure dp,[(p] ensures the 
finiteness of the Schwinger functional and hence depends 
on some renormalisation scale /i, as well as on For 

the sake of simplicity, and for emphasising the structure 
of the results, we use a condensed notation with indices 
a, b that stand for an integration over space-time and a 


summation over internal indices: 

= J d‘^xJ^ix)Mx), (2.2) 

In (2.2) we have implicitly defined the ultra-local metric 
7 “a' = 5{x — x') 7 “/, leaving the internal part 7 “/ unde¬ 
termined. In case fa involves fermionic variables we have 
J°'(j)a 7 ^ The notation as well as some properties of 

the metric 7 “** are detailed in appendix A. In the general 
case (2.1) we consider the coupling of N tensorial fields 
with rank n* < rij+i to the theory. We substitute indices 
a by multi-indices a = an ■ • ■ aim, ■ ■ ■ , o-ni ■ ■ ■ ciNnjv 
with un = rzniax- In the general case, different a^ 
can carry different internal indices, e.g. different repre¬ 
sentations of a gauge group relating to different species 
of fields. This is implicitly understood and we identify 
aij = Qj from now on in a slight abuse of notation. Gon- 
tractions read 

N 

Tl"T 2 a = ^ri“^ , (2.3) 

i=l 

and the generalised metric 7 ^^ is defined as 

N 

(7"‘’)=0(G7r. (2.4) 

i=l 

The definitions in (2.3),(2.4) are nothing but the exten¬ 
sion of the field space to include composite operators 
4>ai---an- The interest in such a general setting is twofold: 
firstly, it allows us to formulate, at all scales, the the¬ 
ory in terms of physically relevant degrees of freedom. 
Secondly, it naturally includes the coupling to compos¬ 
ite operators and related flows. The source term in the 
Schwinger functional (2.1) reads 

N 

= . (2.5) 

i=l 

For Ui = 1 for all i the general source term (2.5) boils 
down to the standard source (2.2). A simple tensorial 
example is given by a = 0,0102 and fs, = {4>a,faia-i) = 
{fa, faifa^) with o = oi = 02 = x, a Scalar field and its 
two-point function. This leads to a source term 

J‘^(j)ei = j d'^xJ{x)f{x) 

+ J d'^xd'^y J{x,y)f{x)f{y). (2.6) 

The above example also emphasises that the sources J'^ 
should be restricted to those sharing the (index-) sym¬ 
metries of the fields fa_- We illustrate this within the 
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above example of a scalar field. The source term for 
^aiaa = Satisfies 

where = i(j“i “2 _|_ j 020 i^ jg symmetric part 

of J. The anti-symmetric part _ i|'joia 2 _ ja 20 i^ 

does not couple to the field, Jf“^“^Vai 02 = 0. Conse¬ 
quently we restrict the sources to the symmetric ones. 
The symmetry properties of a function or <jf - are also 
carried by its derivatives. Again we illustrate this by the 
example introduced above: derivatives w.r.t. the function 
carry its symmetry properties. This entails that 

— F ( \ — (F -iF +F ')') 12 7) 

^j(a) “ ^ yW ~ \f ’ 2 V^,ai02 ' ^ ,a2ai ) J , J 

where = J“. The basic example is the derivative of 
J w.r.t. J. It reads 


SJ^ ^ (b) ^ 




( 2 . 8 ) 


the second entry on the rhs is the identity kernel in 
the symmetric subspace. We also have J = 0 with 
J[“l = 0, and get J^j = (0, i From 

now on we suppress this detail. Derivatives are always 
taken within the appropriate spaces defined by the cor¬ 
responding projections, and carry the related symmetry 
properties. 

Within the above conventions the Schwinger functional 
(2.1) reads 


gW[J] ^ J exp{-S'[(^] + J«0a((^)} . (2.9) 


Many of the structural results presented here can be al¬ 
ready understood within a scalar theory with a single 
field. There we have a. = a = x with the ultra-local 
metric 7“a' = 6{x — x'). In these cases one can simply ig¬ 
nore the additional notational subtleties in the presence 
of fermions and tensorial fields. 

The definition (2.9) is rather formal. For most in¬ 
teracting theories it is impossible to strictly prove the 
non-perturbative existence of c?/r[i^] exp{—S'[i^]}, not to 
mention determining it in a closed form. Here we follow 
a bootstrap approach in simply assuming that a finite 
W\J] exists. This assumption is less bold than it seems 
at first sight. It is merely the statement that the classical 
action S'[(,3] admits a well-defined quantum field theory 
in terms of appropriately chosen fields Then quite 

general normalised expectation values I[J] = {I[J, (/>]) are 
defined by 

I[J] = . (2.10) 


The I include correlation functions that relate to one 
particle irreducible (IPI) as well as connected and dis¬ 
connected Green functions in (p. Subject to the definition 
of (j) this may include A^PI Green functions in the funda¬ 
mental fields As an important sub-class included in 
(2.10) we present normalised Wpoint functions 

N 

=(n'^a.), (2.11a) 

i=l 

with 

i=l 

The correlation functions (2.11) include all moments of 
the Schwinger functional and their knowledge allows the 
construction of the latter. A simple example for (2.11) is 
{(j>), the expectation value of the operator (p coupled to 
the current with = jj. We briefly illustrate the con¬ 
struction of connected or IPI Green functions by an im¬ 
portant example, the full propagator. With the 1-point 
function {(p), the propagator Waia 2 ['^] = (<^ai^a 2 )iPi fol¬ 
lows as /ifas - 

Further important examples are correlation functions 
I where /[J, generates a symmetry of the theory at 
hand. Let us first consider general Dyson-Schwinger (DS) 
equations, formally given by 

J g[0\ (dp.[0\ vI/[<^] = 0 , (2.12) 

with bosonic operator Q. For (2.12) to hold the operator 
g{0\ has to generate a symmetry of the path integral. For 
infinitesimal transformations Q, (2.12) translates into 

/[J]=0, (2.13a) 


with 




(2.13b) 


In (2.13) we have assumed (gdfj,) = 0. We emphasise 
that this can be easily achieved by reducing dfj, to the 
flat measure with dfi = dip'i’i. With dfj,' = dip and 
= v]/^v]/ we substitute The simplest 

relevant example for (2.13) is provided by the standard 
DS equations. They encode translation invariance of the 
flat measure dp. Accordingly, the standard DS equations 
are obtained with (p = ip, g[p] = dp = dp and dt = 1. 
Within this choice we arrive at 


/dse[J] = j - (^) = 0, 
op 


(2.14a) 
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with 


I-dse = J -= -if) ■ (2.14b) 

ocp 

Eq. (2.14) is the well-known functional Dyson-Schwinger 
equation. It assumes a multiplicative renormalisation 
procedure preserving all symmetries (d/i = dip, dt = 1). 
When additive renormalisation is required, or when we 
study a renormalisation procedure breaking the symme¬ 
tries of the classical action, this can be captured in a 
non-trivial 4' 

In case Q generates a symmetry of the action, QS = 0, 
the above relation simplifies. Restricting ourselves also 
to invariant functionals 4^ with Q'i> = 0 we are led to 


/“[J]=0, (2.15a) 

with 


= 4-, (2.15b) 


where a carries the group structure of the symmetry. 
In (2.15) we have used the bosonic nature of Q as well 
as assuming that the symmetry is maintained within the 
quantisation: {Q°^dii) = 0. It is often possible and helpful 
to rewrite symmetries in terms of derivative operators Q 
with = 0. This might necessitate the introduction 
of auxiliary fields. For example, in a gauge theory we 
deal with the BRST symmetry with Q = 5, the BRST 
derivative. We add source terms for Q(j) with 

The Schwinger functional IT = W[J,Q] 
is a functional of both, J and Q, and we are led to 



and 


/. = J" 


5W[J, Q] 
SQ^ 


= 0 . 


(2.16) 


We conclude that the set of I defined in (2.10) provides 
the full information about the quantum theory as it spans 
the set of all correlations functions {O}. In this context 
we emphasise again that not all correlation functions of 
interest are directly given by the correlation functions I, 
a simple example being the propagator Wab = Ikia .2 — 


r(l)dl) 
Jai J^a.2 ■ 


The key object in the present approach is the Schwin¬ 
ger functional of the theory, or some related generating 
functional. Often one concentrates on the Wilsonian ef¬ 
fective action S'eff[4>], the generating functional for am¬ 
putated connected Green functions. It is defined by 


(2.17) 


^ More precisely it is captured in a non-trivial dfi absorbed in 


where = (5^S'/((54>)^[4> = 0]. The advantage of 

working with the Schwinger functional W or Ses is that 
it allows for the most straightforward derivation of func¬ 
tional identities. However, a more tractable object is the 
effective action T, the generating function of IPI Green 
functions of (p = {(p). It is obtained as the Legendre 
transform of W, 

r[<^] =sup(J"((.a-lT[J]) . (2.18) 

j 

Eq. (2.18) includes 7VPI effective actions [161-163] for an 
appropriate choice of (pa[p\- The definition (2.18) leads 
to 


r’"[<^] = 7"bJb<(-), (2.19a) 

kE,a[J] = <(-a(J), (2.19b) 

implying that the field (p is the mean field, (p = (cp). 
In (2.19) we have used that J‘^(pa = (p^Ja = ■ 

The derivatives in (2.19) are taken with respect to the 
variables of P and W respectively, that is T’'^ = and 
kE,a[T] = 1^. Furthermore it follows that 

lT,acr’"‘’ = 7\, (2.20) 

The definition (2.10) and the relation (2.13) translate into 
the corresponding equations in terms of IPI quantities by 
using (2.19),(2.20) as well as 

W[J{<P)] = - r[().], ( 2 . 21 ) 


and 


5 

Jj^ 


W 


,ab 


S(pb ’ 


( 2 . 22 ) 


For composite fields one usually splits up the irreducible 
part of ((pa)- As an example we study 2PI scalar fields 
0a = (0a , 0 ai 0 a 2 ). There we have 0 aia 2 = ( 0 ai 0 a 2 ) = 
0aia2 +0ai0a2 with (pa = {(fa)- Here (p^la 2 the IPI part 
of (paia 2 - This extends to general composite operators 
and we parameterise r^^[0^^] := r[0(0^^)]. The 0^^- 
derivative of P^'^ reads 


rPi-'‘[0Pi] = 0e’y0Pi)7%Jb0'''), (2.23) 

where 0c’‘^(0^'^) stands for the derivative of 0 w.r.t. 
0^^. Within the above 2PI example (2.23) boils down 
to (r^'^’'‘[0^^]) = (J“i “2 ^ JO. _|_ 2J°'^(pb), where we have 
sued that We close with the remark that 

it does not make a difference in the relations of this sec¬ 
tion whether we have tensorial multi-indices a or a vector 
index a. 
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III. FLOWS 

In interacting quantum theories it is hardly possible to 
compute generating functionals, such as the Schwinger 
functional W, in a closed form. In most situations one 
resorts to systematic expansion schemes like perturbation 
theory or the I/iV-expansion that come with a small ex¬ 
pansion parameter. In strongly interacting systems trun¬ 
cations are not supported by a small expansion param¬ 
eter and have to be used with care. In general either 
case requires renormalisation [24, 25]. Renormalisation 
group invariance encodes the independence of physics un¬ 
der general reparameterisations of the theory, or, put dif¬ 
ferently, the physical equivalence of (UV) cut-off proce¬ 
dures. RG invariance can be used to resolve the mo¬ 
mentum dependence of the theory by trading RG scal¬ 
ing for momentum scaling. RG transformations always 
imply the scaling of all parameters of the theory, e.g. 
couplings and masses. In turn, the change of a physical 
parameter is related to an RG rescaling. For example, 
changing the mass-parameter of the theory leads to the 
Gallan-Symanzik equation [26, 27]. Presented as a dif¬ 
ferential equation for a generating functional, e.g. the 
Schwinger functional W or the effective action F, it con¬ 
stitutes a functional RG equation [26]. The momentum 
dependence is more directly resolved by block-spinning 
on the lattice [28]. In the continuum theory this is im¬ 
plemented with a momentum cut-off [1-14] leading to the 
Wilsonian RG. 

The strong interrelations between the different RG 
concepts as well as their physical differences become ap¬ 
parent if presented as Functional Renormalisation Group 
equations for generating functionals. FRG formulations 
are also suitable for both discussing formal aspects as well 
as practical applications. The FRG has been introduced 
with a smooth momentum cut-off for simplifying proofs 
of perturbative renormalisability and the construction of 
effective Lagrangians in [6], see also [9, 31-33]. More 
recently, there has been an increasing interest in FRG 
methods as a computational tool for accessing both per¬ 
turbative as well as non-perturbative physics, initiated 
by [10-14]. The recent success of FRG methods was 
also triggered by formal advances that led to a deeper 
understanding of the FRG, and here we aim at further 
progress in this direction. We close with a brief overview 
on the literature in view of structural aspects: general 
formal advances have been made in [34-97]. Progress in 
the construction of FRG flows in gauge theories has been 


achieved in [98-141]. FRG flows in gravity are inves¬ 
tigated [142-147]. All these formal advances have been 
successfully used within applications, see reviews [15-23]. 


A. Setting 

The starting point of our analysis is the finite renor¬ 
malised Schwinger functional W in (2.9). So far we only 
assumed its existence without offering a method of how 
to compute it. We shall turn the problem of comput¬ 
ing the path integral (2.9) into the task of successively 
integrating out modes, each step being well-defined and 
finite. To that end we modify the Schwinger functional 
as follows: 

^W[J,R] ^ ^-AS[^,R]^W[J] ^ ( 3 ^^) 

where 

= (3.2) 

n 

If used as a regulator, the operator exp—AS' in (3.2) 
should be positive (on exp IF), and AS[^,0] = 0. For 
example, the standard setting is given by a = a, 0“ = (f°' 
and 

ASIA. = 3^1 jjs- («) 

A factor 1/2 on the rhs common in the literature is ab¬ 
sorbed into R. With the restrictions a = a, and 

up to RG subtleties, (3.3) leads to a modification of the 
kinetic term S[(p] in (2.9): S[:/] ^ S[ip]+ More 

generally, (3.3) results in a modification of the propaga¬ 
tion of the field (j) which is possibly composite. Such a 
modification can be used to suppress the propagation of 
^-modes in the path integral. In particular, it allows for 
a simple implementation of a smooth momentum cut-off 
[6, 10-14], An amplitude regularisation has been put for¬ 
ward in [55-58, 77] and relates to AS ~ S' or parts of S, 
which ensures positivity. A specifically simple flow of this 
type is the functional Gallan-Symanzik flow [26, 27]. In 
specific theories, e.g. those with non-linear gauge sym¬ 
metries, more general regulator terms can prove advan¬ 
tageous. AS can also be used to construct boundary RG 
flows, in particular thermal flows [17, 121, 122]. 

General regulator terms AS according to (3.2) involve 
higher order derivatives and derivatives w.r.t. currents 
coupled to composite operators. In this general setting 
a different point of view is more fruitful: the operator 
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exp — AS adds source terms for composite operators to 
the Schwinger functional. For example, in the standard 
case with a = a and (3.3) a source term for (pa^b with 
current is introduced. For the class of positive regu¬ 
lator terms AS[(j), i?] the exponential exp —AS is a posi¬ 
tive operator with spectrum [0,1] on exp W and the cor¬ 
relation functions (2.11). Then, under mild assumptions 
the existence of W[J,R] < hF[J, 0] follows from that of 
hF[J,0] = hF[J]. Consequently exp—AS' can be used for 
suppressing degrees of freedom, more precisely J-modes, 
in the Schwinger functional hF[J]. 

We add that W[J, i?] is not well-defined for general R. 
A simple example is a mass-like R with 
for a scalar theory. Such an insertion leads to an un- 
renormalised Callan-Symanzik flow [26, 27]. The re¬ 
quired renormalisation can be added explicitly via a re¬ 
definition of R°'^ that generates appropriate subtractions. 
This amounts to an explicit construction of a BPHZ-type 
renormalisation which is one way to render the Callan- 
Symanzik flow finite. From now on such a redefinition 
of R is assumed whenever it is necessary; in most cases, 
however, the regulators R generate finite W[J, i?] from 
the outset. A necessary condition for the latter is a suf¬ 
ficiently fast decay of R in the ultraviolet. 

Within this general setting the regulators " in 
(3.2) can be (partially) fermionic, even though AS should 
be kept bosonic (even number of fermions involved). A 
simple example is provided by coupling to a fermion 
4>a- It is in general not possible to commute J-derivatives 
and regulators Due to the (anti-) commutation 

relations of the currents J® only specific tensor structures 
have to be considered for the R-. 

- . ^^a^ai+i jjai ■■■ai+iai • • a^ j-g 

where (—1)®*“^“+! is defined in appendix A. Eq. (3.4) 
expresses the fact that fermionic currents anti-commute, 
ja.ijai +1 _ _jai+ijai^ whereas bosonic currents com¬ 
mute with both, bosonic and fermionic currents, leading 
to J^'+^ = (_i)aiai+i jai+i jai. This Symmetry struc¬ 
ture carries over to derivatives of J^. Hence, in (3.2) only 
that part of R carrying the tensor structure expressed in 
(3.4) contributes. 

For illustration, we again study this setting for the 
standard regulator (3.3) providing a modification of the 
propagator. There it follows from (3.4) that for bosonic 
variables only the symmetric part of the tensor con¬ 
tributes. For the fermionic part only the anti-symmetric 
part is relevant. Here we do not allow for mixed 


(fermionic-bosonic) parts and (3.4) reduces to 


Dab _ oba 

■^bosonic ■^bosonic ’ 


(3.5a) 


and 


Dab _ _D 

-^fermionic ^ 


ba 

fermionic • 


(3.5b) 


The corresponding AS are bosonic. 

So far we have discussed a modification of the Schwin¬ 
ger functional. The Schwinger functional W[J, i?] is only 
one, if important, correlation function. We seek an ex¬ 
tension of (2.10) consistent with (3.1): it should define 
general normalised expectation values in the regularised 
theory as well as allowing for a straightforward extension 
of the symmetry relations I[J] = 0 as given in (2.13a). 
A natural extension is 

I[J, R] = i[J, A] . (3.6) 

Eq. (3.6) entails that /[J, 0] = I[J] and guarantees well- 
defined initial conditions I[J, ooj. Moreover, applying the 
extension (3.6) to a relation I[J] = 0 we are led to 

/[J] = 0 ^ /[J,i?] = 0, Vi?. (3.7) 


Hence a symmetry relation I[J] = 0 is lifted to a symme¬ 
try relation I[J, i?] = 0 in the presence of the regulator. 
Eq. (3.6) can be rewritten solely in terms of W[J,R] as 

I[J, R] = i[J, ^,R] , (3.8a) 


with 

i[J, A, ii] = /[J, A] , (3.8b) 

see also [21]. In case i[J, A] only contains a polynomial 
in J we can easily determine I[J, jj,R\ in a closed form. 
As for i? = 0, the set of all correlation functions {0[J, ii]} 
can be constructed from the set {i[J, ii]}. A general flow 
describes the response of the theory to a variation of the 
source i? and, upon integration, resolves the theory. Such 
flows are provided by derivatives w.r.t. R of correlation 
functions 0[J,R\ in the presence of the regulator 


^i^ai...a„ 


50[J, R] 


(3.9) 


Here is a small variation. Basic examples for 

correlation functions O are the Schwinger functional 
W[J,R\ and the expectation values I[J,R\ defined in 
(3.8). 

In case we define one-parameter flows R{k) that are 
trajectories in the space of regulators i? and hence in 
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theory space, the general derivatives (3.9) provide valu¬ 
able information about the the stability of the chosen 
one-parameter flows, in particular if these flows are sub¬ 
ject to truncations. Stable one-parameter flows can be 
deduced from the condition 




60[J, R] 

^Rai 


= 0 , 


•^stab 


(3.10) 


where is the set of operators that provide a reg- 

ularisation of the theory at some physical cut-off scale 
kes, and i?stab G Eq. (3.10) ensures that the flow 

goes in the direction of steepest descent in case (3.10) de¬ 
scribes a minimum. If flows are studied within given ap¬ 
proximations schemes, the stability condition (3.10) can 
be used to optimise the flow. Note that (3.10), in partic¬ 
ular in finite approximations, does not necessarily lead 
to a single i?stab- Then (3.10) defines a hypersurface of 
stable regulators. We also emphasise that (3.10) cannot 
vanish in all directions 6R except at a stable fixed point 
in theory space. Consequently one has to ensure within 
an optimisation procedure that the variations SR± con¬ 
sidered are orthogonal to the direction of the flow. If this 
is not achieved, no condition is obtained at all. We shall 
come back to the problem of optimisation in section V. 


B. One-parameter flows 

1. Derivation 

In most cases we are primarily interested in the un¬ 
derlying theory at i? = 0, that is 0[J] = (!1[J, 0], e.g. 
in W[J] = W[J,Q], the Schwinger functional of the 
full theory and its moments. Total functional deriva¬ 
tives (3.9) with arbitrary 5R°'^ scan the space of theo¬ 
ries given by W[J,R]. For computing W[J] it is suf¬ 
ficient to study one-parameter flows with regulators R 
depending on a parameter k € [A, 0] with R{k = 0) = 0 
and W[J, i?in], 0[J, i?(A)] well under control. These one- 
parameter flows derive from (3.9) as partial derivatives 
due to variations 


SR = dtdtR, (3.11) 

where t = ln{k/ko) is the logarithmic cut-off scale. The 
normalisation fcg is at our disposal, and a standard choice 
is ko = A leading to tin = 0. In the following we shall 
drop the normalisation. The flows with (3.11) lead to 
correlation functions Ok that connect a well-known initial 
condition at A with correlations functions O = Oq in the 


full theory. In most cases a well-defined initial condition 
is obtained for large regulator R^ oo. This is discussed 
in section III C 4. 

The most-studied one-parameter flow relates to a suc¬ 
cessive integration of momentum modes of the fields (p, 
that is fc is a momentum scale. More specifically, we dis¬ 
cuss regulators leading to an infrared regularisation with 
IR scale k of the theory under investigation, the scale 
k providing the parameter k G [A:in,0]. To that end we 
choose regulator terms AS'[(^] = R°-^ipaPb for a scalar 
theory with 


R = R{p^)S{p-p'), (3.12) 

with the properties 

(i) it has a non-vanishing infrared limit, p'^/k'^ 0, 

typically R ^ k^ for bosonic fields. 

(ii) it vanishes for momenta larger than the cut-off 

scale, for p^/k"^ ^ oo at least with ^ 

0 for bosonic fields. 

(ii) ’ (ii) implies that it vanishes in the limit fc —> 0. In 

this limit, any dependence on R drops out and all 
correlation functions Ok reduce to the correlation 
functions in the full theory O = Oq, in particu¬ 
lar the Schwinger functional Wk and the Legendre 
effective action L^. 

(iii) for fc ^ cx) (or k ^ A with A being some UV 
scale much larger than the relevant physical scales), 
R diverges. Thus, the saddle point approximation 
to the path integral becomes exact and correlation 
functions Ok tend towards their classical values, 
e.g. Tfc^A reduces to the classical action S. 

Property (i) guarantees an infrared regularisation of the 
theory at hand: for small momenta the regulator gen¬ 
erates a mass. Property (ii) guarantees the (ultraviolet) 
definiteness of W[J, i?]. The insertion /AS vanishes in the 
ultraviolet: no further ultraviolet renormalisation is re¬ 
quired, though it might be convenient. It is property (ii) 
that facilitates perturbative proofs or renormalisability. 
Properties (ii)’ and (iii) guarantee well defined initial con¬ 
ditions, and ensure that the full theory as the end-point 
of the flow. In most cases the regulator R = p^r{p'^/k/^) is 
a function of a; = p^/k"^, up to the prefactor carrying the 
dimension. For such regulators the condition (iii) follows 
already from (i). For regulators (3.12) with the proper¬ 
ties (i)-(iii) we can study flows from a well-known initial 
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condition, the classical theory or perturbation theory, to 
the full theory. Integrating the flow resolves the quantum 
theory. The properties (i),(ii) guarantee that the flow is 
local in momentum space leading to well-controlled limits 
a; ^ 0, oo. In turn, mass-like regulators violate condition 
(ii): additional renormalisation is required. Moreover, 
the flow spreads over all momenta which requires some 
care if taking the limits 0,oo, see e.g. [17]. 

General one-parameter flows are deduced from (3.1), 
(3.8) by inserting regulators R{k) where k G [A, 0]. The 
condition R{0) = 0 guarantees that the endpoint of such 
a flow is the full theory. For one-parameter flows, (3.1) 
reads 


with 


A^,[A] = A5[^,i?(fc)], 

and AS' is defined in (3.2). Similarly we rewrite (3.8) as 
Ik[J] = Ik[J, (3.14a) 

with 

W, i[J, . (3.14b) 

We also recall that (3.14) entails that /o['^] = I[J] and 

/[J]=0 ^ /fc[J] = 0 V/c, (3.15) 

that is a symmetry relation I[J] =0 is lifted to a rela¬ 
tion Ik[J] = 0 in the presence of the cut-off. The flow of 
/c-dependent quantities Ik, dtik with t = \nk at fixed cur¬ 
rent J allows us to compute I[J], if the initial condition 
Ia is under control. For momentum flows, this input is 
the high momentum part of / at some large initial scale 
A. Perturbation theory is applicable for large scales, and 
hence Ia[J] is well under control. The flow equation dtR 
can be evaluated with (3.6) for R{k). However, for later 
purpose it is more convenient to approach this question 
as follows. Let us study the operators 


fall into the class of functionals (2.10) and can be lifted 
to their i?-dependent analogues (3.8), and in particular 
to F’fc,7fc,A/fe as defined in (3.14). The full Schwinger 
functional W[J] = Wo[J] is independent of t, dtW = 0, 
and we derive from (3.6) that F = A/ and consequently 

Fk = Ah ■ (3.18) 

Moreover, the most interesting / are expectation values 
in the full theory and do not depend on t. For this class 
we have A/ = 0 leading to Fk = 0. Still, the consid¬ 
eration of more general Fk will also prove useful so we 
do not restrict ourselves to Fk = 0. The general Fk is 
derived from (3.14b) with help of 

[dt, A-"" (3.19) 

In (3.19) we have used that [dt , ^] = 0 as 9* = dt\j. 
The rhs of (3.19) commutes with AS'fc[^] and we con¬ 
clude that (dt + AS”)^ ^ R])exp—ASk = {exp—ASk)dt. 
Inserting F into (3.14b) and using (3.19) we are led to 
Fk with 

Fk = (9t + AS[^ , 7?]) Ik . (3.20) 

The expression in the parenthesis in (3.20) is an operator 
acting on everything to the right. Inserting (3.20) into 
(3.14a) we arrive at 

+ AS[A , 7?]) h = Ah , (3.21) 

valid for general Ik given by (3.14). Aik on the right 
hand side carries the explicit t-scaling of the operator I 
and vanishes for t-independent /. In order to get rid of 
the exponentials in (3.21) we use that + 

^1^). With this relation (3.21) turns into 

-l- Wk F AS'[^ -|- (j ), 7?]^ h = Aik, (3.22) 

where we have introduced the expectation value </> = {(j))j 
of the operator coupled to the current 


F[J,A]=dti[J,A]^ 

(3.16) 

AI = [dt, i]. 

(3.17) 


Here the t-derivative acts on everything to the right, 
i.e. 9t/G[J] = {dtI)G[J] + IdtG[J], and is taken at 
fixed J. The notation for partial derivatives is ex¬ 
plained in appendix B. The functionals I, F and AI 


MJ] ■■= Wk,4J]. (3.23) 

Eq. (3.22) involves the flow of the Schwinger functional, 
Wk, reflecting the normalisation of h- Independent 
flows of Ik are achieved by dividing out the flow of the 
Schwinger functional. The flow Wk is extracted from 
(3.22) for the choice /fc = 1 with Ah = 0, following from 
1=1 and AI = [dt , /] = 0. Then, (3.22) boils down to 
Wk + (AS'[^ -I- (f), 7?]) = 0, where both expressions are 



functionals and not operators. More explicitly it reads 

+ (3-24) 

n 

X• • • (^ + </-),„_, s^)Wk[J] = 0. 

Eq. (3.24) is the flow equation for the Schwinger func¬ 
tional. It links the flow of the Schwinger functional, 114, 
to a combination of connected Green functions II4^ai .a„- 
For quadratic regulators (3.3) we obtain the standard 
flow equation for the Schwinger functional, 

(dt + R‘^^^jj^+R^'^<Pe.^)Wk[J]=0. (3.25) 

We remark for comparison that the standard notation 
involves a factor 4 in the i?-terms. It has been shown 
in [47] that (3.25) is the most general form of a one loop 
equation. Eq. (3.24) makes this explicit in a more general 
setting as the one considered in [47]. Only flows depend¬ 
ing on 144_ai -a„ with n <2 contain one loop diagrams in 
the full propagator. Note in this context that J couples 
to a general operator not necessarily to the field. 

Eq. (3.24) is the statement that the flow operator 
AS'i[J, R] = Wk + -\- (j) ^ R\ with 

A5i[J, R] = A5[A +4>,R]- (A5[^ + <^, 7?]), (3.26) 

is given by all terms in AS'[^ -I- , 77] with at least one 

derivative jj acting to the right. For later use we also 
define AS'„[J, 77] as the part of AS' with at least n J- 
derivatives. Their dehnitions and properties are detailed 
in appendix C. The operator of interest here, ASi, can 
be written with an explicit J-derivative as 

ASi[J,7?] = AS’"[J,7?]^ . (3.27) 

The operator AS‘‘[J, 77] is defined in (C.l). Using (3.24) 
and the definition (3.27) in (3.22) we arrive at 

+ AS"[J, 7?]^^ h = Aik , (3.28) 

valid for general Ik, Aik given by (3.14). Aik carries 
the explicit t-scaling of I and is derived from (3.17). 
The partial 7-derivative is taken at fixed current J. The 
flow of a general functional 4 requires the knowledge 
of (j)a[J] = Wk,a.[J] and Aik- Only for those R that 
entail this information in a closed form, (j) = (j)[Ik\ and 
Aik = AIk[Ik], the flow equation (3.28) can be used 
without further input except that of / a - 


2. Flow of the Schwinger functional 

We proceed by describing the flow (3.28) for correlation 
functions (3.14) within basic examples. To begin with, we 
study the flow of the Schwinger functional M4. First we 
note that its flow (3.24) was derived from (3.22) with I = 
1. The final representation (3.28) was indeed achieved 
by dividing out (3.24). Nonetheless, the latter should 
follow from the general flow equation (3.28). Naively 
one would assume that R = 114 can be obtained from 
a 7-independent operator J, that is AI = 0. However, 
inserting the assumption R = Wk into the flow (3.28) 
and using (3.24) we are led to 

AR = AS^[J, 77]((.a - {AS[A + 7?]). (3.29) 

which does not vanish for all J, e.g. for quadratic regula¬ 
tors it reads Aik = 77‘‘'^^a'^b- Hence (3.29) proves that 
R = Wk implies AI 4 0. Indeed in general (3.29) can¬ 
not be deduced from a AI that is polynomial in the cur¬ 
rent and its derivatives. The above argument highlights 
the necessity of the restriction of (3.28) to functionals R 
constructed from (3.14). Still the flow equation for Wk 
can be extracted as follows. Let us study the flow of 
(7fe)a = W4,a = 4>a. which also is of interest as (j) is an 
input in the flow (3.28). R = (j> falls into the allowed 
class of R as 

/a=(4)a=^ ^ (4)a = lUfc.a = <()a . (3.30) 

Moreover, AR = 0. Consequently, the flow of the func¬ 
tional 4 introduced in (3.30) reads 

lUfe,a + A^[^ + , 7?](/)a - (A5[^ + , 7?])<^a = 0 . 

(3.31) 

With ^1=0 the second term on the left hand side can 
be rewritten as follows 

A5[A+^,77]((.a = AS[A+fg,R]{A+<f)^ 

= {^ + + cj,, 77]. (3.32) 

We emphasise that the first line in (3.32) is not an oper¬ 
ator identity. For the second line in (3.32) we have used 
the bosonic nature of the regulator term and the rep¬ 
resentation jj + (j) = This also entails that 

(/)a(A5'[A -I- 4]) = (AS'[^ -I- (j), 77])(()a. We have al¬ 

ready mentioned that dt4>a[J] yf 0 as the 7-derivative is 
taken at fixed J. For the same reason we can commute 7- 
derivatives with J-derivatives: dtWk,a[J] = {dtWk[J]),a- 
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We conclude that the flow of Wk,ii can be written as a 
total derivative 

dtWk + {AS[^ + <l>, R])] =0, (3.33) 

which upon integration yields 

dtWk + {AS[A + cj) ^ R]) = 0. (3.34) 

Eq. (3.34) agrees with (3.24) 

3. Standard flow 

For its importance within applications we also discuss 
the standard quadratic flow. In this case the flow (3.28) 
reduces to 

(a* + 4[J] = 0, (3.35) 

and (3.29) turns into Aik = R^^'t’hfpa which does not 
vanish for (j) ^ 0. That proves that there is no / leading 
to Ik = Wk- The flow of {Ik)a = Wk,a follows as 

{dtWk[J]),a = 

= - (Wfc,bc + 0b<^c)l . (3.36) 

J ,a 

Both sides in (3.36) are total derivatives w.r.t. J^. Inte¬ 
gration leads to 

Wk[J] = (VEfc.ab + <^a0b) , (3.37) 

where we have put the integration constant to zero. For 
the reordering in (3.37) we have used that the regulator 
is bosonic. Eq. (3.37) agrees with (3.25). It also 
follows straightforwardly from (3.34) for quadratic regu¬ 
lators. 

4- Flow of amputated correlation functions 

The results of the previous sections translate directly 
into similar ones for amputated correlation functions 
Ik [J {(j))] with the following fc-dependent choice of the cur¬ 
rent 

J-($) = [S + = {Pk)baJ'^ , (3.38) 


^ We have fixed the integration constant to precisely match (3.24). 


introducing the classical propagator Pk- With (3.38) the 
flow for general correlation functions Ofc[J($)] is com¬ 
puted as 

dtOklJm = IdtOklJ] + $a(9tA5)itoOfe.bl 

L JJ=J($) 

(3.39) 

in particular valid for Ok = Ik- The t-derivative on the 
Ihs of (3.39) is taken at fixed 4>: the first term on the 
rhs of (3.39) is the flow (3.28) at fixed J, and the second 
term stems from the fc-dependence of T($). For example, 
in the presence of a regulator the effective Lagrangian 
S'eff[4>] (2.17) turns into 

:=-Wfe[J(<I>)], (3.40) 

and hence has the flow (3.39) with (3.24). This flow 
further simplifies for quadratic regulators R^'°(j)a.(t>h- For 
this choice we arrive at 

= \{PkU - 5 ,^, S’X - . 

(3.41) 

Often (3.41) is rewritten in terms of the interaction part 
of the effective Lagrangian defined as S'lnt* = •S'eff;, -I- 
^ [S' -I- AS];^“q fl’afl’b- The flow of Smt^ follows as 

dtSintM = l{Pk)ab (S(„t - 5;:,, sQ , (3.42) 

where we dropped the <I>-independent term —{dt lu-Pkpa.- 
Flows for Seff;, and its A-point insertions can be found 
e.g. in [6, 9, 12, 13, 36]. They are closely related to 
Callan-Symanzik equations for A-point insertions for 
i? (x with a possible mass renormalisation, see also 
[181]. The flows (3.41), (3.42) can be extended to <I>- 
dependent Pk by using the general DS equations (2.12) 
in the presence of a regulator, see e.g. [40, 135]. Then it 
also nicely encodes reparameterisation invariance. 

We close this section with a remark on the structure 
of the flows (3.28),(3.39). They equate the scale deriva¬ 
tive of a correlation function to powers of held deriva¬ 
tives of the same correlation function. The latter are 
unbounded, and the boundedness of the flow must come 
from a cancellation between the different terms. Hence, 
within truncations the question of numerical stability of 
these flows arises, see [69]. 

C. Flows in terms of mean fields 

1. Derivation 

In most situations it is advantageous to work with the 
flow of IPI quantities like the effective action, formulated 
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as functionals of the mean field (ps, = W^g,. In other words, 
we would like to trade the dependence on the current J 
and its derivative in (3.28) for one on the expectation 
value (j) and its derivative Similarly to (2.18) we define 
the effective action F = i?] as 

T[(j), R] = sup (J®(/)a - W[J, R]) - AS'[(j), R]. (3.43) 

J 

where 

(3.44) 

n>2 

The exclusion of the linear regulator terms in A5'(, is nec¬ 
essary as they simply would remove the dependence on 
the linear regulator. F)^, R] is the Legendre transform of 
W[J,R], where the cut-off term has been subtracted for 
convenience. For R ^ 0 (3.43) reduces to (2.18). The 
definition (3.43) constrains the possible choices of the op¬ 
erators coupled to J to those which at least locally admit 
a Legendre transform of W[J,R\. Eq. (3.43) implies 

7"bJ‘’ = (F + A5')^", <^a = lF,a, (3.45) 

as well as 

Gac(r-i-A^)’‘^‘’ = 7\, (3.46) 

with 


Gac = LF^ac . (3.47) 

Here 7 ’^a leads to the minus sign in fermionic loops, see 
appendix A. For quadratic regulators (3.3) the above 
relations read 

7\jb =F’=" + 2i?'^‘’0b, (3.48) 

and 


Gac(r’"‘’ + 2R^'^) = 7^. (3.49) 


For (3.48),(3.49) we have used (3.5) and the bosonic na¬ 
ture of The operator in (3.46) is the full field 
dependent propagator. With (3.46) we are able to relate 
derivatives w.r.t. J to those w.r.t (j) via 


A = r A 

S(l)b ’ 


(3.50) 


where we have used that 0b,a = lF_ab = Gab- As in the 
case of the Schwinger functional we are not only inter¬ 
ested in the flow of F but in that of general correlation 
functions I as functions of 0. This is achieved by defining 
I[J, R] as a functional of J[0]: 


/[0,i?]=/[J(0),i?]. (3.51) 


We emphasise that / is not necessarily IPI, it only is for¬ 
mulated in terms of such quantities. Still, all IPI quan¬ 
tities can be constructed from the class of I. 

One-parameter flows for / are derived by using trajec¬ 
tories R{k). We extend the notation introduced in the 
last section for flows of I with 

/fc[0] = /[0,i?(fc)], (3.52) 

and 

Ffc[0] = F[0,i?(fc)]. (3.53) 

For reformulating (3.28) in terms of R we need the re¬ 
lation between dtR = dt\(j)Ik and dtik = dt\jlk, see also 
appendix B. With (3.50) we rewrite Ik,a = Gab-ffe’*^, and 
it follows from (3.51) that 

0t/fc[0] = dtlklJ] + (a* J"[0]) Gab/fc’^’, (3.54) 

with dtJ^^lcj)] = Now we insert the flow for R, 

(3.28), in (3.54). With (3.50) the operator AS'i[J, A] = 
AS'’^[J, 1?]^ is rewritten in terms of Gab^. As it is 
more convenient to use an expansion in plain derivatives 
^ we also employ the identity A^i [J, i?] = AAi [0, 7?], 
the terms that contain at least one derivative w.r.t. J 
are equivalent to those containing at least one derivative 
w.r.t. 0. Note that this fails to be true for higher deriva¬ 
tive terms, ASn with n > 1. Together with (3.50), (3.54) 
the above considerations lead to the flow (3.28) as an 
equation for R 

(a* - ((atJ")Gab + A5b[0,i?]) 4 ^) 4[0] = A4[0], 

(3.55) 

where AS'b is defined in (C.lb). It can be easily com¬ 
puted for general regulators. However, the higher the 
order of derivatives is in the regulator term, the more 
loop terms are contained in AS'i. For further illustration 
we have detailed the simplest case of the standard flow 
in appendix D. We proceed by evaluating (3.55) for a 
specific simple R: we use 4[0] = 0 already introduced 
via 4 = jja in (3.30). For this choice we have A4 = 0 
and dtR = 0, and the flow (3.55) reads 

4J"Gab-(AAb[0,i?])=O. (3.56) 

Here (AS'b[0,7?]) is the linear expansion coefficient of 
ASi in a power expansion in derivatives w.r.t. 0, see 
also (C.l). Note that (3.56) already comprises the flow 
equation for Tk'. it follows from the definition of the 
current in (3.45) that J'^7®b = (r -I- A5")’®. Moreover 
4(rfe’‘^) = (dtTk)'’^ as the partial t-derivative is taken at 
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fixed 4>. Then (3.56) contracted with (T^ + com¬ 

prises 9t(rfc b) and is a total derivative w.r.t. (j) which 
can be trivially integrated. This can be best seen for 
quadratic regulators (3.3) for which (3.56) boils down to 

J" + (Gi?G)bcrfc’"‘’‘' 7"d - = 0, (3.57) 

see also (D.l) in appendix D. We also remark that an 
alternative derivation of the identity (3.56) solely makes 
use of structural considerations which prove useful for 
general flows: for IPI Ik the related term in (3.55) is 
not IPI, whereas the other terms are. Accordingly these 
terms have to vanish separately which implies that 
the expression in the parenthesis has to vanish leading to 
(3.56). With dfJ^Gah = — (AAb)^, 7?]) the coefficient of 
^ in (3.55) takes the form 

AAb[((-,i?] - (AAb[((-,i?]) = A^ab[<(-,i?]47 , (3.58) 

where AS'abjl^ is the part of the operator AS'b contain¬ 
ing at least one (()-derivative. AS'ab follows from (3.58), 
see also (C.lb). With (3.58) the operator in the flow 
(3.55) is 

A52[0,i?] = A^ab[<^,i?]x74- , (3.59) 

OfPaOfph 

that part of AS'[G^ -I- <f), i?] containing at least two (jy 
derivatives, and we arrive at 

(dt + A52[<^, 7?]) hm = Aik , (3.60) 

for general functionals Ik as defined with (3.14) and 
(3.51). The functional Aik originates in the explicit t- 
scaling of I. The partial t-derivative on the left hand 
side of (3.60) is taken at fixed 0, and the operator AS2, 
(3.59), accounts for inserting the regulator R into the 
Green functions contained in correlation functions Ik- 
We also provide a representation of AS2 [(p, R] that only 
makes direct use of AS'[GA -j- cp ^ i?], 

AS2[<P,R] = AS[GA + ^^r] (3.61) 

-([A5[GA +^^R], - {AS[GA +<P,R]), 

where (G^)i 3 = Gbcj^- The relatively simple insertion 
operator AS '2 in terms of derivatives w.r.t. <p is related to 
the structural dependence of Ik on (p and R that is fixed 


^ Strictly speaking, one also has to use that the span of IPI I 
generates all IPI quantities. 


by the definitions (3.14),(3.51). In turn, changing the 
definition of Ik, Ik leads to different flows. The construc¬ 
tion oi Ik, Ik is a natural one as it includes general Green 
functions (0") as building blocks. Still, it might be worth 
exploring the flows of different correlation functions for 
specific problems, whose setting admit more natural vari¬ 
ables than the Ik- 

Let us now come back to the remark on numerical sta¬ 
bility at the end of section IIIB. In contradistinction 
to the flows (3.28),(3.39) the flow (3.60) relates the scale 
derivative of a correlation function to a polynomial of the 
full propagator, field derivatives of the effective action 
and the correlation function itself. In most cases both 
sides of the flow (3.60) are bounded, ensuring numerical 
stability and hence better convergence towards physics 
[69]. A notable exception is the case where the Legendre 
transform from Wk to Tfc -l- AS'j. is singular. This either 
hints at a badly chosen truncation, or it relates to phys¬ 
ical singularities that show up in the propagator G, see 
also [49]. In the scale-regime where such a singularity 
occurs one might switch back to the flow of Wk or SeS^ 
[84]. In the vicinity of S'eff.ab ~ 0 the flows (3.41),(3.42) 
are bounded. 


2- Flow of the effective action 


As in the case of the flow equation for Ik we describe 
the content and the restrictions of (3.60) within basic 
examples. From its definition (3.43) it follows that its 
flow is closely related to that of Wk, 

dWkicp] - dtWk[J] - AS'[(P, R] = 0, (3.62) 

where we have used (3.45) for J^{<pa — Wk,a[J]) = 0. 
Inserting the flow (3.34) for the Schwinger functional we 
are led to 

dtTk[4>] - {AS[GA -F(P,R]) + AS'[P),R] =0. (3.63) 
More explicitly it reads 

atrfc[<^] - ^ ■■■"- 

n>2 


^ (G^ + ^)ai • • • (G^ + <(')a„_i 


= 0 . 


(3.64) 


The explicit form of the flow (3.64) allows us to read off 
the one particle irreducibility of T^)^] as a consequence 
of that of the classical action £'[(()]: the flow preserves 
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irreducibility and hence it follows recursively from that 
of 

As for the Schwinger functional there is no / with A/ = 
0 leading to /^ = Ffc. The related consistency equation 
reads 

Ah[^] = AS2[(t>,R]rk 

+iAS[G^ + <l>,R])-AS'[<l>,R]. (3.65) 

The right hand side of (3.65) does not vanish for all (p 
implying AR ^ 0. Moreover, in general (3.65) cannot 
be deduced from a A1 polynomial in the current J and 
its derivatives. Again this highlights the necessity of re¬ 
stricting R to those constructed from (3.14) and (3.51). 

Similarly to the derivation of the flow of 114 we 
can derive the flow (3.63) from that of its derivative, 
We use I® = The corresponding Ik de¬ 
rived from (3.14b) as I® = i?]. The 

second operator AS'’®^ originates from the commutator 
term [AS', J*^]. The commutator gives the right (p- 
derivative of AS[^ + p, i?] at fixed J, see appendix C. 
Contracted with y^b we arrive at the left derivative, 
where we have also used the bosonic nature of AS. The 
corresponding R reads with (3.45) 

ik = Tfe’" + AS'’"[0, R] - (AS-"[G^ + P,R]), (3.66) 

Moreover, AR = 0. The choice (3.66) boils down to 
R = Tfc’® in the standard case. For general flows the 
last term on the right hand side of (3.66) is non-trivial 
by itself. Indeed, its flow can be separately studied and 
follows from / = AS’"^)^ , R] and F = ftAS’'^)^ , R]. 
This leads to AI = AS’‘^[^ , i?] and AR = (AS’‘‘[G^-I- 
p, S]). Inserting this into the flow (3.60) we are led to 

(St + AS2[</.,!?]) (AS'"[G^ + <^, R]) 

= iAS’‘^[G±+p,R]). (3.67) 

The above equation describes the flow of the functional 
(AS’®[G^ + p, i?]) at fixed second argument R. Using 
(3.67) within the flow of R of (3.66) it reads 

dRk'^ = -AS2[P,R]{Tk + AS'y’^ 

-AS'-"[</), R] + (AS’"[G^ + P,R]). (3.68) 

Eq. (3.68) looks rather complicated. However, note that 
AS 2 acts on the current as (Ffc -f AS'(,)’‘‘ = see 

(3.45). Hence the evaluation of (3.68) is simplified if rep¬ 
resenting AS '2 [p, i?] in terms of J-derivatives as all higher 


J-derivatives vanish. To that end we use that the sum of 
all derivative terms in either p or J coincide as in both 
cases it is given by the operator AS — {AS). The latter 
can be written as the sum of all terms with two and more 
derivatives, AS '2 and the linear derivative terms, ASa^|- 
and AS‘‘[J, A] jjiT respectively. This leads us to 

AS2[P,R] = AS2\J,R] 

-(ASa[<^, A])4;^ + (ASV,^])^- (3.69) 

The validity of (3.69) follows from the above considera¬ 
tions, but also can be directly proven by inserting (3.50) 
in the first term on the right hand side. Using the rep¬ 
resentation (3.69) of AS 2 [^, A] in (3.68), only the terms 
in the second line of (3.69) survive as (AS 2 [J, A] J) = 0. 
Furthermore AS‘^[J, i?]^y‘^b4'^ = (AS’‘‘[G^ + P, 7?]), 
and (3.68) reduces to 

dtVk^^ = {ASb[P, R]) (Ffe + ASO’^*" - AS'^^[P,R] . 


(3.70) 


Both terms on the right hand side of (3.70) are total 
derivatives w.r.t. pg,. For the first term this follows with 
(3.46) and it reduces to (AS'[G^ -I- p, 7?])’‘^. With this 
observation we arrive at 


dtTk’^ = 


{AS[G^ + P, 7?])-A5'[().,i?] ’ , (3.71) 


which upon integration yields (3.63). 


3. Standard flow 

The standard flow relates to regularisations 
quadratic in the fields pg. We also restrict ourselves to 
bosonic i?’s, that is no mixing of fermionic and bosonic 
fields in the regulator. Then, the flow of R can be di¬ 
rectly read off from (3.60) 

dJkiP] + {GRG)bJF^'^[P] =0. (3.72) 

The flow equation for F^ is extracted from Ik^ = F^’®. 
This R can be constructed from = ^'°'^Jb'. we get 
ik^ — 7*’“Tb—jyb- Inserting this operator into (3.14) 
we arrive at R^ = ^^’^Jh — R^'°Ph = Its flow is read 
off from (3.60) as 

atFfe’- = -(GAG)bcFfe’'^‘’"= [Gbci?^’"]’", (3.73) 

where we again have used (3.46). The flow (3.73) matches 
(3.56) and can be trivially integrated in p^ 

Ffc = Gbci?^’" , (3.74) 


13 



where we have put the integration constant to zero. 
Eq. (3.74) is the standard flow equation of as derived 
in [10] (up to the normalisation | absorbed in R). It 
matches the flow of Wk, (3.34), when using (3.46) and 
the definition of in (3.53) 

OtTkl^] = -dtWklJ] - Gbc . (3.75) 

In (3.75) we have used that — Wk,a) = 0, see 

(3.23). Note that we could have used (3.75) instead of 
evaluating I]^ = (j) for deriving (3.57) with help of 

= at(rfe'") + 7?-‘’0b 

= (a^Ffc)’- + 7?-‘>0b . (3.76) 

The derivatives in (3.76) commute as the partial t- 
derivative is taken at fixed (j). Indeed, it is the flow of 
the Schwinger functional Wk which is at the root of both 
derivations. The flow of Wk equals that of the effec¬ 
tive Lagrangian SeSk generating amputated connected 
Green functions. The relation between the flows StT k and 
in particular the (in-)equivalence within trunca¬ 
tions, has been explored in [12, 13, 39, 69, 88, 89], see 
also the reviews [16, 17, 19-22]. The numerical stability 
of the flows has been compared in [69]. 

Finally let us study the consistency condition (3.65) 
in the present case. It reads dWk + {G RG)\^W’k'^ = 0, 
which does not match (3.74). Hence there is no Ik leading 
to Ik = Ffc and AJj, = 0. Again this underlines the im¬ 
portance of (3.14) for devising flows: first one constructs 
an Ik or R from (3.14). Their flow is given by (3.28) and 

(3.60) respectively. 

4- Initial condition for general flows 

For IPI correlation functions R the Ihs of (3.60) con¬ 
sists of IPI graphs in the full propagator G. Furthermore, 

(3.60) is only one loop exact if d'l/S.Sk[x\ does not depend 
on X, that is for n < 2, see also [47]. For n = 2 the flow 

(3.60) boils down to the flow (3.72), whereas dtR = 0 for 
n = 1. For n > 2 we have higher loop terms in (3.60). 
Appropriately chosen ^11 loops finite. In 

the class of R that provide momentum cut-offs, these 
loops can be localised about the cut-off scale. Then the 
flows (3.60) are finite and numerically tractable, shar¬ 
ing most of the advantages with the standard flow (3.72) 
with n = 2. Indeed, for specific physical problems, in 
particular theories with non-linear symmetries, the gen¬ 
eral choice in (3.60) can pay-off. However, we emphasise 


that for general flows the limit R ^ oo has to be studied 
carefully. Here, it is understood that R ^ oo entails a 
specific limit procedure characterised by some parame¬ 
ter, i.e. the standard momentum regularisation Rk —> oo 
for k ^ A. For practical purposes an accessible limit of 
the effective action F^ is required as it usually serves as 
the initial condition for the flow. In particular regulator 
terms AS'fc[(/f)] that, after appropriate held rescaling, tend 
towards finite expressions which are more than quadratic 
in the fields require some care. The general case can be 
classified as follows. For a regularisation AS' of a theory 
with classical action S[:^] and a given limit procedure 
R oo we can find field transformations </> ^ f{R)4> 
with f{R cx)) = 0 that render S -I- AS finite: 

hm (S + AS)[/(i?)0] =S[0]. (3.77) 

R^oo 

For R that diverge for all (^-modes S only depends on AS. 
In the standard case with (j) = (p and S[(/>] = 
with held-independent S®*'’, the effective action F^ tends 
towards the classical action S of the theory In general, 
the corresponding effective action F k tends towards 

F[</), R] ^ SMcf)] + det |^(<^) + m , (3.78a) 

where F is given by 

f[(/)] = - In , (3.78b) 

and [d(f] is the flat ())-measure including renormalisation 
effects. The term in the exponent comprises the Taylor 
expansion of S[()) -I- </>] about (p leaving out the first two 
terms, 

. (3.79) 

n! 

n>2 

The representation (3.78b) relates F to a Wilsonian ef¬ 
fective action. We emphasise that {(j>) ^ p, the mean 
field computed from the path integral (3.78b) is not the 
original mean held. Indeed we compute 

^"[<^] = (<('a-(<^b))S'‘’". (3.80) 

Eq. (3.80) also entails that F has no classical part due to 
the classical action 5[^] in the exponent in (3.78b). Only 


^ More precisely all power-counting irrelevant couplings tend to 
zero. 
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those limits {ASk,S) admitting the computation of the 
effective action F in (3.78b) provide suitable initial condi¬ 
tions for the flow (3.60). They lead to consistent flows as 
defined in [47]. The standard limit {ASk, leads 

to a (()-independent F: the explicit integration of (3.78b) 
gives ^ In det up to renormalisation terms stemming 

from [d(jji\. Such a flow was coined complete flow in [47] 
as it connects the classical action with the full effective 
action. The flow provides for the complete integration 
of quantum effects, and the theory is determined by the 
parameters in the classical action S. The requirement 
of convexity of the effective action constrains the set of 
parameters in S, which can be evaluated with help of the 
regulator dependence, see [49]. The general case (3.78) 
with non-trivial, but accessible F was coined consistent 
flow. Eq. (3.78) also covers the interesting class of proper¬ 
time flows [50-54], where F comprises a full non-trivial 
quantum theory [45-48]. A detailed discussion of the 
general situation will be given elsewhere. 


D. General variations 


In the previous sections we have studied one-parameter 
flows (3.11). These flows can be used to compute ob¬ 
servables in the full theory starting from simple initial 
conditions like the classical or perturbation theory. For 
the question of stability of the flow or its dependence on 
background fields present in the regulator we are also in¬ 
terested in general variations (3.9) of the regulator. In 
particular functional optimisation as introduced in sec¬ 
tion V is based on studying general variations w.r.t. R. 
These variations are also useful for the investigation of 
physical instabilities [49]. They can be straightforwardly 
derived with the generalisation of (3.16): 


F[J, A] = ^ (3 3^) 


with 


A/V,^] = 






J[J,; 


■■■a„ ’ SJ 


. (3.82) 


The corresponding F[J, ^,i?] follows with the commu¬ 
tator 




6 


obi • -bn 6 
^ SJ^i 



= SR’^^ 


•an S 
6J^i 


5 


(3.83) 


as 


F[J,A^r] (3.84) 

With (3.81) and (3.84) the derivation of one-parameter 
flows in the previous sections directly carries over to the 
present case. Therefore we read off the response of R 
and Ik to general variations from (3.28) and (3.60) re¬ 
spectively: 

+^S,[J,SR]^ I[J,R] 

= AI[J, 6R] , (3.85) 

and 

+ A52[</>, SR]^ i[ct>, R] 

= Ai[<j3,SR], (3.86) 

with Ai[4), SR] = AI[J{(j)), (5i?]. For the choice R = R{k) 
and 6R = dtR the flows (3.85), (3.86) reduce to the one- 
parameter flows (3.28), (3.60). 


IV. RENORMALISATION GROUP FLOWS 
A. RG flows of general correlation functions 

The flows (3.28) and (3.60) comprise the successive 
integrating-out of degrees of freedom in a general quan¬ 
tum theory. The standard example is the integration of 
momentum modes, but the formalism allows for more 
general definition of modes. The current J and the regu¬ 
lator R couple to which is not necessarily the fun¬ 

damental field (j) = (p. In any case, with R we have in¬ 
troduced a further scale fc, thus modifying the RG prop¬ 
erties of the theory. Moreover, at any infinitesimal flow 
step k ^ k — Ak there is a natural fc-dependent repa- 
rameterisation of the degrees of freedom. Taking this 
reparameterisation into account should improve numer¬ 
ical stability. Hence the appropriate choice of Ja at the 
initial scale A is affected by the proper book keeping 
of the anomalous scaling, which becomes crucial in the 
presence of fine-tuning problems. It also is relevant for 
studying fixed point solutions of the flow. Hence the 
representation of RG rescalings in the presence of a reg¬ 
ulator is a much-studied subject, e.g. [34-44, 114, 115]. 
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From the formal point of view canonical transformations 
on the functional space govern both RG rescalings and 
general flows presented here. This point of view shall 
be developed elsewhere. In most practical applications 
an appropriate fc-dependent RG rescaling is simply in¬ 
corporated by hand, see reviews [15-23]. We emphasise 
that contrary to claims in the literature the incorpora¬ 
tion of RG rescalings is not a matter of consistency but 
rather one of numerical stability and optimisation. We 
will come back to this issue later in chapter V. 

The formalism introduced in the previous chapter al¬ 
lows us to derive RG equations in the presence of the reg¬ 
ulator. In general we deal with theories that depend on a 
number of fundamental couplings g, which also includes 
mass parameters. We are interested in the response of 
the theory to an infinitesimal total scale change of some 
scale s, e.g. s = k, the flow parameter, or s = /x, where 
/i is the renormalisation group scale of the full theory. 
The couplings and the currents may depend on this scale, 
g = g{s), = J^{s). An infinitesimal variation is intro¬ 

duced by the operator sjf. Here we consider a general 
linear operator Ds with 

Ds = sds + Ig'jgtdg^ + , (4.1) 


with 


DsW = 0, 

where the partial s-derivative is taken at fixed J (and 
g), see appendix B, and the anomalous dimensions 7 j 
do not mix fermionic and bosonic currents. With J- 
independent matrices 7 we only consider linear depen¬ 
dences of the currents. More general relations are eas¬ 
ily introduced but should be studied separately in the 
specific situation that requires such a setting. Still we 
remark that non-linear relations can be reduced to lin¬ 
ear ones by coupling additional composite operators to 
the currents. A relevant non-trivial example for (4.1) is 
sds = with renormalisation scale (or cut-off scale) 
g of IF[J, 0] and Jgjjj the corresponding anomalous di¬ 
mensions of couplings and fields respectively. We also 
could use sds = gdg + dt ■ We emphasise that the oper¬ 
ator Ds accounts for more than multiplicative renormal¬ 
isation. The matrices Jgjjj are not necessarily diagonal 
and the multi-index a possibly includes composite oper¬ 
ators. Hence (4.1) naturally includes the renormalisation 
of composite operators (e.g. in A^PI flows) or effects due 
to additive renormalisation. The operator Ds does not 
commute with derivatives w.r.t. J. Still DsW = 0 can 


be easily lifted to identities for general A^-point functions 
with 


= (T>slF).ai .a„ 

N 

~ 7j aikk.ai -ai-i bai+i - a^f 5 (4-2) 

i=l 

where we have used the commutator 


[Ds , = -T/ajjrr • (4.3) 

The derivation of Dg-Q-Ows for functionals Ik is done along 
the same lines as that of the t-flow in section HI. First 
we define an operator F similarly to (3.16) with 

F = Dsi and A/ = [Dg , /]. (4.4) 


With DgWk = 0 it follows that Fk = Aik which does 
not vanish in general. We shall use that still Aik = 0 for 
/ = 1. The only further input needed is the commutator 
of the regulator term AS' with the differential operator 
Dg defined in (4.1). For its determination we compute 

r a 7 b 5 pai - an S S 1 

[7j sja , -ft ^jai • • • 5jan J 

= (4.5) 


where we have used the symmetry properties (3.4) of 
i?. Eq. (4.5) enables us to compute the commutator 
[Hg, AS]. For the sake of brevity we introduce a short 
hand notation for the symmetrised contraction of 7 with 
R, 


n 

( 7 ^ T)'"!-’"" = , (4.6) 

i=l 


for a given n. The commutator of AS with the differen¬ 
tial operator Dg takes the simple form 


AS[^,R] 


AS[^ , {Dg - 7 ,)i?]. (4.7) 


With the above preparations the derivation of the RG 
flow boils down to simply replacing R in the commutator 
(3.19) with {Dg — 7 j)i? and allowing for a non-zero Fk = 
Aik- We finally arrive at 


{Dg + ASi[A , {Dg - 7 ,)S]) Ik = Aik , (4.8) 

where AI = [Dg,I]. The term Aik contains the s-scaling 
inflicted by the operator I, and ASiIk contains the addi¬ 
tional scaling inflicted by the operator AS. In summary 
(4.8) comprises general scalings in the presence of the 
regulator, and reduces to the flow (3.28) for s = fc, up 
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to an additional fc-dependent RG rescaling. We also em¬ 
phasise that for the derivation of (4.8) only the linearity 
of the operator Ds has been used. 

An explicit example for the content of (4.8) is pro¬ 
vided by the RG equation of A^-point functions = 
as defined in (2.11). Then Dg = imple¬ 
menting RG rescalings in the full theory. Furthermore 
we assume that the operator AS' does not spoil the RG 
invariance of the theory, i.e. the commutator (4.7) van¬ 
ishes. The requirements on the regulator R leading to a 
vanishing commutator are further evaluated in the next 
section IV B. The RG equation for Ik is read off from 
(4.8) as {D^ + Af7j'^ai)(4^\a2 - ajv = 0i where A4 
produces the explicit scaling Ay of the A-point function. 
This is the usual RG equation for A-point functions as 
expected. For connected A-point functions it is put down 
in (4.2). 

The general equation (4.8) simplifies in the case of 
quadratic regulators, 

[{Ds - 7.) sw) 4 = A4, (4.9) 

where 

[{Ds - = DsR^'^ - 27 /c R^^ . (4.10) 

In the last equality in (4.10) we have used R^^^ = 
(_l)ab 7 ^ba and the fact that jj does not mix fermionic 
and bosonic currents. The general s-scaling of the 
Schwinger functional for quadratic regulator is derived 
similarly to the flow (3.37): we use / = jA which leads 
to Aik = —7j'^abFfc^b- Moreover we have DsWk,a. + 
7j'^abFfe,b = {DsWk).ei- Inserting this into (4.8) we arrive 
at 

[DsWk + i(Gbc + <(>b4) [{Ds - 7.,)^]^’"] = 0. (4.11) 

Upon integration we are led to 

DsWk = -i(Gab + <^a<)>b) [{Ds - 7j)Rr'^. (4.12) 

Eq. (4.12) entails the response of the theory to a general 
scaling including the flow (3.28) as well as RG rescalings. 
For s = /X (4.12) expresses the modification of the RG 
equation Df^W[J, 0] = 0 in the presence of the regulator. 

B. RG flows in terms of mean fields 

We proceed by turning (4.8) into an equation formu¬ 
lated in terms of IPI quantities and fields. This is done 


by repeating the steps in the derivation of (3.60), and 
hence we shorten the details. First we lift (3.54) to op¬ 
erators Ds- This requires the definition of the action of 
Ds on functionals F[<j3\ as provided in appendix B: 

Ds = {sds + lg''39^^g^ + 7<^'"a0b ■ (4.13) 

With Ik = Ik[J{4’)] we rewrite Dglk in terms of Ik as 

Dsik[^] = DsIk[J] + {{Ds - 7.) J[<('])"Gab , 

(4.14) 

where ( 7 jJ)® = In (4.14) we have used that 

Ds\jlk = Dsh - {ljJ[(l)\)'^h,a- In case Dg stands for 
a total derivative w.r.t. s, the second term on the right 
hand side of (4.14) has to vanish, [Dg — 7 j)d = 0. Then, 
keeping track of dependences on </> or J is irrelevant. 
With (4.14) we get 

{Ds + A52[((-, {Ds - 7,)R]) ik[ct>] = Aik - Ab 4’^’ , 

(4.15a) 

where 

Ab = (D,J)"Gab + (A^b[<(-, {Ds-^,)R]). (4.15b) 

We emphasise that in Dglk is at our disposal. Now, 
as in the case of the t-flow for 4 , we simplify the above 
equation by solving it for R = (j) following from R = jj- 
Then, DgR = and AR = —^j4>. This leads to 

Ab =-(70 + 7j)4(/>a. (4.16) 

Inserting this into (4.15) the As-fiow equation for R 
reads 

{Ds + Aa^|^)4 + AS' 2 [^ , {Ds — ^j)R\ik = Aik , 

(4.17) 

where 

Dg + Aaj^ = sdg + 7 g j9i^gj ~ 7j a ‘^b ■ (4-18) 

The dependence on y^ has completely dropped out. Its 
role has been taken over by — 7 j. In other words, however 
we choose the fields </> to scale under Ds^ the RG flow 
(4.15) shows its natural RG scaling induced by DgW = 0 
and DgJ = y^ J. For the t-flows studied in section III this 
translates into dt(j) = 0 , corresponding to the natural 
choice 70 = 0. As y^ is at our disposal we take the 
natural choice 

'Y'P = -Ij > ( 4 - 19 ) 
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for which Ab = 0. With the choice (4.19) we arrive at 

{D, + AS 2 [(/-, {D, + 7^)i?]) Ik [4>] = Aik M , (4.20) 

where AJ derived from (3.14) with AI = [Dg,!]- 
Eq. (4.20) is the ^-based representation of (4.8), and 
hence comprises general explicit and implicit scalings in 
the presence of the regulator. A special case are those s- 
scalings with [Dg + "frj,)R = 0 leading to AS' 2 /fe = 0. For 
these choices of the pairs (i?, Dg) the s-scaling of the reg¬ 
ularised theory remains unchanged in the presence of the 
regulator. If Dg stands for a scale-symmetry of the full 
theory such as the RG invariance with s = /i, regulators 
with {Dg + = 0 preserve the RG properties of the 

full theory, see [42, 43]. We shall discuss this interesting 
point later in section VIIIB. 

The above equations (4.8),(4.20) can be straightfor¬ 
wardly lifted to include general variations (3.85),(3.86) 
by 

Dg^DR=5R^^-^-j^^^^ +SsDg (4.21) 

with variations 6R{k) about R{k) and 6s{R,6R),s{R). 
The operator Dr stands for the total derivative w.r.t. R, 
hence using Dr in (3.86) simply amounts to rewriting a 
total derivative w.r.t R in terms of partial derivatives. 
These general variations are important if it comes to sta¬ 
bility considerations of the flow as well as discussing fixed 
point properties. 

We close this section by illustrating the content of the 
RG flow (4.20) within some examples. First we note that 
by following the lines of the derivation for the t-flow of 
Tfe, (3.63), we can derive the RG flow of the effective 
action. It is given with the substitutions dt —>■ Dg and 
R {Dg -I- 7^)1? in (3.63). For quadratic regulators 
(4.20) reduces to 

Dgik + i(G [{Dg + j4,)R] G)ab/fc’"‘’ = Aik , (4.22) 
where 

[{Dg - j^)Rr^ = DgR^^^ + 27,^"e . (4.23) 

The Hg-flow of the effective action Ffc is derived with the 
choice i’^ = This leads to = Tfe’'^ and A/® = 

dj’^bTfe’*’. By also using the commutator [Dg,-^] = 
we are led to 

- iGbc[(i?« + ,a = 0 • (4-24) 

This is trivially integrated and we arrive at 

DgVk = ^Gb4{Dg + 70)i?]‘’" , (4.25) 


where we have set the integration constant to zero. The 
Ihs of (4.25) can be projected onto the anomalous di¬ 
mensions 7 with appropriate derivatives w.r.t. fields and 
momenta. Then the rhs is some linear combination of 
7 ’s. These relations can be solved for the 7 ’s, see e.g. 
[36, 42, 43]. With the choice s = fj, and (I?^ -I- 70 )i? = 0 
we are led to the equation Df^Fk = 0 , the regularised 
effective action satisfies the RG equation of the full the¬ 
ory. This interesting case is further discussed in sec¬ 
tion VIIIB. 


V. OPTIMISATION 

An important aspect concerns the optimisation of 
truncated flows. Optimised flows should lead to results 
as close as possible to the full theory within each order 
of a given systematic truncation scheme. This is inti¬ 
mately linked to numerical stability and the convergence 
of results towards physics as already mentioned in the 
context of RG rescalings in the last section. By now a 
large number of conceptual advances have been accumu¬ 
lated [60-71], and are detailed in sections VB, VG. In 
particular [64] offers a structural approach towards op¬ 
timisation which allows for a construction of optimised 
regulators within general truncation schemes. Still a fully 
satisfactory set-up requires further work. In the present 
section we take a functional approach, which allows us 
to introduce a general setting in which optimisation can 
accessed. This is used to derive a functional optimisation 
criterion, which admits the construction of optimised reg¬ 
ulators as well as providing a basis for further advances. 


A. Setting 

The present derivation of flows is based on the exis¬ 
tence of a finite Schwinger functional W and finite corre¬ 
lation functions 0\4>\ for the full theory. These quantities 
are modified by the action of an R-dependent operator, 
G[^, R] with 0[<j)[ = G[()), 0], see section III A. 
One-parameter flows (3.86) connect initial conditions, 
that are well under control, with the full theory. For 
most theories these flows can only be solved within ap¬ 
proximations. Typically truncated results for correlation 
functions 0[(j), 0 ] show some dependence on the chosen 
flow trajectory R{k) not present for full flows by defini¬ 
tion. Naturally the question arises whether we can single 
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out regulators R{k) that minimise this non-physical reg¬ 
ulator dependence. 

Consider a general systematic truncation scheme: at 
each order of this systematic expansion we include addi¬ 
tional independent operators to our theory, thus succes¬ 
sively increasing the number of independent correlation 
functions. At each expansion step these correlation func¬ 
tions take a range of regulator-dependent values. This 
regulator dependence should be rather small if the trun¬ 
cation scheme is well adapted to the physics under inves¬ 
tigation. In extremal cases the truncation scheme may 
only work for a sub-set of well-adapted regulators but fail 
for others. An optimisation of the truncation scheme is 
achieved if at each successive expansion step and for the 
set of correlation functions included in this step we ar¬ 
rive at values that are as close as possible to the physical 
ones of the full theory. In all cases such an optimisation 
of the truncation scheme is wished for as it increases the 
reliability and accuracy of the results, in the extremal 
case discussed above it even is mandatory. 

General correlation functions 0[<j)] are either given di¬ 
rectly by /[(/)] or can be constructed from them as the 
/ include all moments of the Schwinger functional, , 
see (2.11). From now on we restrict ourselves to /[</<]. 
Most relations directly generalise to correlation functions 
O [(/)], in particular to physical observables, except those 
whose derivation exploits the flows of I. The constraint 
of quickest convergence can be cast into the form of an 
equation on the single iteration steps within a given trun¬ 
cation scheme. We expand a correlation function R] 
in orders of the truncation 

/«[<^, R] = R] + R ], (5.1) 

where adds the contribution of the zth order. With 

adding the subscript k and keeping the variable R we wish 
to make explicit the two qualitatively different aspects 
of the i?-dependence of I^^^[(j),R]. Firstly, the I^'''>[(f),R] 
depend on the functional form of R{k) that singles out 
a path in theory space. Secondly, k is specifying that 
point on the path belonging to the value k of the cut¬ 
off scale ranging from k/A€ [0,1]. If we could endow 
the space of theories with a metric, optimisation could 
be discussed locally as a stationary constraint at each k. 
The resulting flows are geodesic flows, and k turns into a 
geodesic parameter. For now we put aside the problem of 
defining a natural metric or norm on the space of theories, 
but we shall come back to this important point later. 

The full correlation function in the physical theory is 
given by /[^] = /g°°^ [()), i?] and shows no i?-dependence 


except for a possible i?-dependent renormalisation group 
reparameterisation, not present for RG invariant quan¬ 
tities. Therefore, optimisation of a correlation function 
/ at a given order i of an expansion scheme is simply 
minimising the difference 

OO 

min||/[())]=min|| ^ A(”)io|| , (5.2) 

^ ^ ' n=i+l 

on the space of one-parameter flows R{k). An optimal 
trajectory i?opt(fc) is one where the minimum (5.2) is 
achieved. As already mentioned in the last paragraph, 
for the general discussion we leave aside the subtlety of 
specifying the norm ||.||. The constraint (5.2) also fixes 
the freedom of RG rescalings for a given I with fixed RG 
scheme in the full theory. 

How can such an optimisation (5.2) be achieved? A 
priori we cannot estimate how close to physics the re¬ 
sults are, that were obtained with a specific regulator 
and truncation step. If we could, we knew the physical 
results in the first place and there would be no need for 
any computation. Hence an optimisation of the ith order 
within a general truncation scheme has to be based ei¬ 
ther on structural aspects of the flow or on an evaluation 
of successive truncation steps; both procedures allow to 
evaluate (5.2) within the given ith order. For correlation 
functions I with 

OO OO 

II ^ A(")/||= ^ ||A(")/||, (5.3) 

n=i-\-l n=i+l 

we can reduce (5.2) to a constraint on at a given order 
i. The minimum in (5.2) is approached for regulators 
minimising each term ||A("Hj| separately. In this case 
optimised regulators i?opt(fc) are those with 

||A«/o[</),i?opt(fc)]|| =min||A«/o[0,R(A:)]||, (5.4) 

R{k) 

for almost all i, (j)- Eq- (5-4) is the wished for relation 
applicable at each order of the truncation. Note that 
(5.4) also eliminates the freedom of a /c-dependent RG 
scaling of general correlation functions. It picks out that 
implicit RG scaling which minimises the norm of A*^*Ho. 
One could argue that an optimisation with (5.4) possi¬ 
bly gives close to optimal convergence even if (5.3) is not 
strictly valid: in the vicinity of optimal regulators sub¬ 
leading orders I — are small in comparison to the 

leading rest term A^®^/ and a partial cancellation between 
them should not have a big impact on the optimisation. 
Still it is dangerous to rely on such a scenario. For its 
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importance we discuss the general situation more explic¬ 
itly: assume that we deal with Wmax observables 
TO = 1,TOmax, built off some set of I[(j), i?]’s. Examples 
are critical exponents, physical masses, particle widths 
etc.. Within the fth order of a given truncation scheme 
and a flow trajectory R{k) we get Am [i?] taking values in 
an interval . By construction the extrem- 

isation picks out either or This procedure 

entails an optimisation if A^^^" ^ (sub¬ 

ject to the correct choice of the closest extremum). In 
turn, if AP(*y® G A““^A] g, procedure picking out 

the boundary points decouples from optimisation, only 
by chance it provides close to optimal results. Indeed this 
scenario is likely to be the standard situation at higher 
order of the truncation scheme. An indication for this 
case is the failure of finding coinciding extrema for all 
observables, in particular if these extrema are far apart. 
The resolution of this problem calls for an observable- 
independent optimisation based on (5.2). 

The evaluation of the optimisation (5.2) is more conve¬ 
nient in a differential form. This equation can be directly 
derived from (5.2). However, there exists an alternative 
point of view which might also be fruitful: truncated 
flows may be amended with functional relations valid in 
the full theory. The hope is to carry over some additional 
information from the full theory that is not present in the 
truncation of the flow. This is the idea behind the use 
of symmetry relations such as STIs together with flows. 
In the context of optimisation the key relation is the reg¬ 
ulator independence of the full theory. By construction 
the end-points of one-parameter flows = Io[4>i di] are 
correlation functions in the full theory, being trivially 
independent of the path R{k) in regulator space: k is 
a further variable of R and any local variation of such a 
path about a regulator -anfe change R]- 

Moreover, in section IV we have seen that there is the 
freedom of /c-dependent RG scalings of the full theory, 
and the apparent independence of /[(()] = Io[4>i -R] on the 
path R(k) for full flows is expressed in the relation 

= (5-5) 

for all The variation on the Ihs of (5.5) stands for 
the total derivative w.r.t. g^jg^ including possi¬ 

ble i?-dependent RG scalings as in Dji, (4.21). The rhs 
of (5.5) accounts for a possible integrated i?-dependence 
of the renormalisation scheme at fc = 0: 6fj,{R, SR), fJ,{R). 
For RG invariant /[</>] the rhs of (5.5) vanishes. For RG 
variant /[</>] the rhs can always be absorbed in an ap¬ 


propriate redefinition of the variation w.r.t. R, though 
technically this might be difficult. The relation of (5.5) 
to the optimisation (5.2) is provided by enforcing (5.5) 
already for the ith order of the truncation scheme and 
absorbing the RG scaling on the rhs in an appropriate 
redefinition of the R-variation. Also assuming (5.3) we 
are led to 


SR’^^- 


^, S\\io[<P,R]-ll,^\<P,R]\\ 

(5Rai-a„fc 


= 0 , 


(5.6) 


which is the differential form of (5.2). Eq. (5.5) is an inte- 
grability condition for the flow. Its relation to reparame- 
terisations of the flow and the initial condition R(A)] 
become more evident by using 

/•O flh 

/o [</), R\ = I A [cS, R] + —dth [0, R \. (5.7) 

Inserting (5.7) in (5.5) leads to 




R(A) 



DRi[(k, R] 


R(k) 


= 6 (In^) D^iQ[4),R ], 


(5.8) 


with Dr defined in (4.21). The integrand in (5.8) is a 
total derivative, and with using that SR\r^o = Sfj, the 
Ihs in (5.8) equals the rhs. A variation of the initial 
regulator R(A) in general entails that /[</>, R(A)] cannot 
be kept fixed by adjusting an appropriate RG scaling. 
For example, a different momentum dependence of R(A) 
leads to different composite operators coupled to the the¬ 
ory via AS, and hence physically different theories. For 
sufficiently large regulators these differences are usually 
sub-leading. Neglecting this subtlety we conclude that 
in general a change of regulator with a vanishing rhs and 
fixed initial conditions I[(j), R(A)] entails a /c-dependent 
RG scaling of the flow. 


B. Principle of Minimum Sensitivity 


For the sake of simplicity we only discuss couplings A’s 
and not general functionals I or O. Eq. (5.6), evaluated 
for one or several observables A^, to = l,...,TOi„ax) at 
some order i of a given truncation scheme can be viewed 
as a constraint for truncated flows. This implies the 
search for local extrema of observables Am in regulator 
space. However, not knowing Aphys we have to resort to 
(5.5), most conveniently written as 


JRai...a„fe 


SXm 

(5Rai -a„fe 


= 0 . 


(5.9) 
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Eq. (5.9) can be seen as a symmetry constraint as sug¬ 
gested in the last section or as an optimisation with the 
assumption (5.3). As a constraint, (5.9) can have sev¬ 
eral solutions or none (the extremum could be a point on 
the boundary in regulator space). Eq. (5.9) in its integral 
form, only allowing global changes along the full flow tra¬ 
jectory, is related to the principle of minimum sensitivity 
(PMS) [59], which has been introduced to the functional 
RG in [60], for further applications see [61-63]. Its lim¬ 
itations have been discussed in [66]. Practically such a 
PMS extremum has been evaluated by computing observ¬ 
ables Ai,..., for a class of regulators R{ai,aj) 

labelled with ai,...,aj. Strictly speaking, mmax should 
increase with the order i of the truncation, as the num¬ 
ber of observables increase with the order i of the trun¬ 
cation scheme. The functional derivatives w.r.t. R turn 
into ordinary ones and we are left with the problem of 
finding a coinciding extremum for these A. As already 
mentioned before, even if they exist at all, these extrema 
need not coincide. There are several options of how to 
proceed in such a situation. We can constrain the set of 
regulators by fixing the value of some Ai,..., A^ to their 
physical value to all orders of the truncation, thereby 
sacrificing a part of the predictive power. Such a pro¬ 
cedure resolves (if r is big enough) the above mentioned 
problem and the optimisation is done for the other ob¬ 
servables Ar+i,..., Am„,ax this smaller set of regulators, 
see [60]. One also could argue that optimised values for 
each of these variables are obtained at their extrema. A 
regulator that optimises the flow of Ai is not necessarily 
optimising that for other Am. This idea has been used 
in [63] and in general requires the use of supplementary 
constraints. Both procedures have to be used with care 
as already discussed in general in the last section V A. 
Within the present explicit procedure this analysis hints 
at several short-comings: firstly, fixing the values of r ob¬ 
servables does not necessarily lead to small flow operators 
AS' 2 , and possibly constrains the values for Ar+i,...,Am 
to regions that are far from their physical values. Sec¬ 
ondly, non-coinciding optimal regulators also could hint 
at a badly working truncation scheme, or badly chosen 
Am- We emphasise again that searching for a solution of 
(5.9) for some variable equals an optimisation (5.2) 
only as long as the physical value is not included 

in the range of possible values of It is mainly for 

this reason that an observable-independent optimisation 
is wished for. 


C. Stability criterion 

The above mentioned problems are also directly re¬ 
lated to the fact that the preceeding use of (5.5),(5.9) is 
not a constructive one; it does not allow us to devise an 
optimal regulator that limits the contribution of higher 
orders of the truncation by construction. Moreover, an 
optimisation as in section VB always involves consider¬ 
able numerical effort. A constructive optimisation cri¬ 
terion, directly based on the fundamental optimisation 
condition (5.2) and on the structure of the functional RG, 
has first been suggested in [64]. The construction there 
also emphasises the link between optimisation, optimal 
convergence and global stability of the flows. We shall 
show later in section V D that the criterion developed in 
[64-68] relates to the local use of (5.5). 

The key point in [64] is the observation that optimi¬ 
sation of any systematic expansion implies quickest con¬ 
vergence of the expansion towards physics. Gonsequently 
we can turn the question of optimisation into that of 
quickest convergence. The latter allows to devise con¬ 
structive optimisation conditions. In [64] it was pointed 
out that for the standard flow (3.73) any such expansion 
includes an expansion in powers of the propagator G = 
1/ (r^ ^[(j)] + R). Hence minimising the norm of the propa¬ 
gator G relates to stability and fastest convergence. Gon- 
sider regulators introducing an IR cut-off with R = R{p^) 
as discussed at the end of section III A. The norm implic¬ 
itly used in [64] is the operator norm on L2' llG[())o, i?] ]] = 
suP||^II^^=i{11G[((.o,R]V’IIlJ, where [[V’lli^ = (/ 
is the L 2 -norm. The norm ||G[())o, i?] [12,2 is directly re¬ 
lated to the biggest spectral value of G aX <j)o, and hence 
is sensitive on the growth of the maximum of G" for 
n —>■ 00. A canonical choice for (jjQ is a field maximising 
llG[(;i', i?]ll on the space of fields (p. Within a truncation 
scheme that uses an expansion in powers of the field a 
natural choice for (j)o is the expansion point. Reformu¬ 
lating the optimisation criterion of [64] in the present 
setting leads to 

{Rstab} = {r with ||G[0 o,R]||l2 ^ \\G[(I)o,R']\\l 2 

V R' and R\klff) = R(fceff) = cfceff} ■ (5.10) 

The normalisation constant c is at our disposal. The 
condition R'{k'^g) = ck^g is required for identifying a 
parameter k'(kes) at which the norm of the propagator 
is taken. Eq. (5.10) allows to construct optimised regu¬ 
lators for general truncations schemes, even though the 
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key demand of stability might necessitate supplementary 
constraints, see e.g. section VIIIE. At a given order it 
singles out a set of stability inducing regulators as (5.10) 
does not restrict the shape of i?stab- An optimisation 
with (5.10) entails in the limit of large truncation order 
the PMS condition (5.5), if the latter admits a solution 
[66]. If the PMS condition has several solutions, by con¬ 
struction (5.10) is likely to pick out that closer to the 
physical value. 

The criterion (5.10) has very successfully been applied 
to the derivative expansion [67, 68], where also the above 
statements have been checked. In its leading order, the 
local potential approximation (LPA), a particularly sim¬ 
ple optimised regulator is provided by 

Roptip^) = {k"^ - - P^), (5.11) 

where 6 is the Heaviside step function. By now (5.11) is 
the standard choice in the field. It is a solution of (5.10) 
with fcgfj = and c = 1. As a solution of (5.10) in 
LPA it only is optimised for the LPA but not beyond, 
as has been already remarked in [65]. Beyond LPA a so¬ 
lution to (5.10) has to meet the necessary condition of 
differentiability to the given order. The related supple¬ 
mentary constraint is provided in (8.42). Solutions to 
(5.11) with (8.42) exist, being simple enhancements of 
(5.11) [71]. We add that (5.11) works within truncation 
schemes where the full momentum dependence of corre¬ 
lation functions is included from the onset. 


D. Functional optimisation 

In summary much has been achieved for our under¬ 
standing as well as the applicability of optimisation pro¬ 
cedures within the functional RG. Still, the situation is 
not fully satisfactory, in particular given its key impor¬ 
tance for the reliability of functional RG methods. In the 
present section we exploit the functional equation (5.5) to 
devise an optimisation criterion based on stability as well 
as discussing in more detail the link between stability- 
related criteria and the PMS condition. We also aim at 
the presentation of fundamental relations and concepts 
that are possibly helpful for making further progress in 
this area. 


1. Local optimisation 

So far we have only discussed the implications of (5.5) 
in its integrated form as done within the PMS optimi¬ 


sation in section VB. Such a procedure always requires 
the integration of the flow and hence involves consider¬ 
able numerical effort. On the practical side, the classes of 
regulators usually used for the PMS are not sufficiently 
dense for resolving the local structure: for the standard 
choice of a momentum regulator we parameterise quite 
generally R{p^) = p^r{x) with x = p^ /k"^. Then, a varia¬ 
tion of i? is a variation of r and as such an integral condi¬ 
tion as it implies a variation at all scales k. Consequently 
a resolution of the local (in k and ai) information of (5.5) 
is only obtained for regulator classes {i?} which include 
as differences R 1 —R 2 smeared out versions of the delta 
function in k: (i?i — i? 2 )“^ oc <5e(fc — fceff)Ai?®i '‘^". 
It is convenient to include these variations functionally: 
evaluating (5.5) for variations local in k we turn (5.5) 
into a local condition on I[(j),R]. As such it is the local 
form of the integrability condition (5.5) and can be read 
off from (5.7) and (5.8), 

6i[(j),R]=0, (5.12) 


the integral in (5.12) describing a small closed curve in 
the space of regulators. Within truncations, (5.12) is a 
non-trivial, physically relevant constraint. For example, 
gradient flows cease to be gradient flows within trunca¬ 
tions that violate (5.12). In turn, this property is kept 
intact if satisfying (5.12). A consequence of (5.5) and its 
local form (5.12) is 




si[<i>,Rm 

(5i?ai-a„fc' 


DRi[(j),R{k)], 


(5.13) 


for all I[(f>,R{k)] and variations SR that vanish at A. 
The right hand side in (5.13) accounts for a total scale 
variation of the end-point R{k) with Dg as defined in 
(4.21). We emphasise again that (5.12) and (5.13) are 
non-trivial constraints within truncations. Moreover, at 
finite fc yf 0, 00 the rhs in general does not agree with 
S (In/i) D^I[(j>, R{k)] even for full flows, as already men¬ 
tioned in section III A: firstly, a general variation w.r.t. 
R leads to the flow (4.20) with (4.21), a special case 
being the one parameter flow (3.60) with SR = dkdkR 
and Dg = dt- Secondly, in the presence of two differ¬ 
ent regulator functions R, R' at some fixed scales k, k' 
the two regularised theories cannot completely agree as 
they differ by their coupling to different composite oper¬ 
ators AS[(j),R] and AS[(j),R']. Still it might be possible 
to identify hyper-surfaces of regularised theories at the 
same physical cut-off scale kes- So far k was just a pa¬ 
rameter labelling one-parameter flows, only its end-point 
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k = Q (and to some extend R = oo) defining a specific 
theory. For k ^ 0 this is a priori not clear, the trivial 
example being two momentum regularisations R{p^) and 
R'{p^) = R{c^p^)/c^. Obviously k cannot be the physical 
cut-off scale in both cases. In this trivial case it is sim¬ 
ple to identify the relative effective cut-off scale for R, R' 
with k = kes and k'{kes) = ckes- In general the natural 
relation k'{k) is less obvious, apart from not being unique 
anyway. Nevertheless let us assume for the moment that 
we have overcome this subtlety. Then we can define a 
variation of R on hyper-surfaces = {i?(A:(A:eff))} 

regularising the theory under investigation at the same 
physical cut-off scale kes- Stability of the flow is achieved 
by minimising its action on the set and (5.13) 

translates into 




di[4), R] 

^^ai-a„fc' 


—-^stab 


51n/r£>^/[((),i?stab], 


(5.14) 


lifting (5.5) to non-vanishing regulators. Eq. (5.14) is a 
non-trivial constraint already for full flows. Subject to a 
given foliation of the space of theories with {i?±} for all 
cut-off scales kgS, (5-14) entails maximal (in-)stability of 
the flow at its solutions i?stab- With (5.13) we rewrite 

(5.14) as 


Dr±I[4’j -R] 


= 0 , 

R —-Rstab 


(5.15) 


where we have absorbed the RG rescaling on the rhs of 

(5.14) in = Dr{5R = 5R±I). A solution i?stab(fc) 
of (5.15) is achieved by varying the flows of variables I 
in regulator space. In its form (5.15) it cannot be used 
to construct regulators i?stab- To that end we have to 
rewrite (5.15) as a criterion on the flow operator AS' 2 . 
This is done as follows: if a one-parameter flow = 
I[(j),R{k)] obeys the constraint (5.15) for all k, so must 
dtik- Varying dtik with SR± it follows with (3.60) and 

(5.15) that 


(^Dr^as2[<i>,r]) i[<i>,R] 


^—-Rstab 


= 0, (5.16) 


where we have used that dt and commute up to RG 
scalings. For most practical purposes the RG scaling will 
be neglected and (5.16) boils down to 






i[<P,R] 


= 0. (5.17) 


^—-Rstab 


Finding a globally stable one-parameter flow i?stab(fc) 
amounts to demanding the validity of (5.16) for all I 


and k. This implies that the variation of AS '2 in the di¬ 
rections 6R± has to vanish at all scales k and all index 
values ai • • • an, that is pointwise zero. Clearly there is 
the danger of overconstraining the regulator. In practical 
applications we limit ourselves to a restricted set of / for 
which we solve (5.16). As any truncation scheme is based 
on the assumption of dominance of certain degrees of free¬ 
dom the related {/rei} should be taken. Then the choices 
Rstab(A) lead to extrema of the action of AS 2 [^, i?stab] 
on {/rei} for all scales k. Such a flow, if it exists, is either 
most stable (minimal AS' 2 ) or most unstable (maximal 
AS 2 ). Eq. (5.16) implements the PMS condition (5.5) 
on {/rei}, as the /-flow vanishes identically at fc = 0 and 
integrating (5.16) over all scales still is zero. We also em¬ 
phasise that (5.16) defines local (in-)stability. We could 
have global extrema at the boundary of the hyper-surface 
{i?_L} defined with fceff- 


2. Optimisation and effective cut-off scale 

So far we have not fixed the hyper-surfaces {/?j_} which 
amounts to the definition of a metric on the space of reg¬ 
ularised theories. Before embarking on a discussion of 
natural definitions of such metrics we would like to elu¬ 
cidate the subtleties within a simple example: assume 
we restrict ourselves to the set of regulators given by a 
specific flow Rhase{k) and possibly momentum dependent 
RG rescaling of /?base(/)- Then the definition of a nat¬ 
ural (relative) physical cut-off scale is uniquely possible; 
the set of regulators {R±}k is defined by those regula¬ 
tors with correlation functions I[(j), R] that only differ by 
RG rescalings (fixed physics) from I[(j), Rbase]- Note in 
this context that the RG scalings also change the field 
4>. The {R±}k cover the restricted space of regulators 
we started with, and by definition (5.16) is satisfied for 
all R G {R±}k- This should be the case as their physi¬ 
cal content is indistinguishable. In turn, if we had cho¬ 
sen another foliation the result would have been differ¬ 
ent. Then, necessarily R{k),R{ck) G {R±}k for at least 
one regulator R and (5.16) differentiates between them 
even though the one-parameter flows R{ck) and R{k) are 
the same. Suitable foliations are those where the hyper¬ 
surfaces {R±} do not contain such pathologies. 

So far k is only a parameter that provides a scale order¬ 
ing without identifying physical scales (except for fc = 0). 
Gonsequently we have to answer the question of how to 
define the distance d of two points R and R' in theory 
space given by their set of correlation functions I[(j),R], 
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I[<f), R'], or more generally 0[(j), i?], 0[(j), R']. To that end 
we define 

do[R,R'] = sup{\\0[<l>,R]-0[^,R']\\}, (5.18) 

where the supremum is taken in an appropriate space 
of fields S, and we have to specify an appropriate 
norm ||.||. A natural choice for S is the configura¬ 
tion space of the theory under investigation. How¬ 
ever, the definition (5.18) only is useful if do is finite 
for almost all R, R'. This can be achieved by turn¬ 
ing O f{0) in an operator or functional that has a 
spectrum that is bounded from below and above, e.g. 
O \/{C + |Op) with positive constant C. Alterna¬ 
tively, one can restrict the space of fields (j), e.g. with 
<!> e Sc = W\\ \m,R]\\ , \\0[cP,R']\\ < C}. Here, the 
constant C < oo is introduced to get rid of singular fields 
with 0[(j),R] = oo that possibly would render the dis¬ 
tance c? = oo for all R, R'. Obviously allowing for these 
fields would spoil the construction. We could also evalu¬ 
ate the norm in (5.18) for a specific configuration <f) = (j)o 
with S = {(('o}- This is an appropriate choice if (j)o could 
be singled out by the truncation scheme, e.g. as the ex¬ 
pansion point in an expansion in powers of the field. 

As general flows (4.20) for Ik depend on T^"^ via AS 2 
which is the crucial input for the optimisation, a natural 
choice for O is the effective action 0[(j), R] = r[^, R] — 
r[0, i?] or its second derivative T^^^. Of course, any 
correlation function R (or set of correlation functions) 
that entails the full information about the theory and 
has no explicit regulator dependence is as good as the 
above suggestion. From now on we drop the subscript o, 
keeping it only if discussing a specific choice for O. The 
distance d between two regularisation paths R{k), R'{k') 
of a theory at the effective cut-off scale k = kes is given 

by 

d[R, R']{k) = min d[R{k), R'{k')] , (5.19) 

k' 

which implicitly defines the relative effective cut-off scale 
k'{k) as that k' for which the minimum (5.19) is obtained 

d[R{k),R!{k'{k))] = d[R,R!]{k). (5.20) 

In general d[R, R']{k) = d[R', R\{k'(k)) yf d[R',R\{k). A 
priori, k'{k) is not necessarily continuous. Indeed one 


® In (3.74) we have put an integration constant to zero, here we 
choose it to be —r[0, it]. At finite temperature the effective ac¬ 
tion r cannot be renormalised that way as r[0, iJ] is related to 
the thermal pressure. 


can even construct pathological regulators that lead to 
discontinuities in k'{k). In most theories such subtleties 
are avoided by using regularity restrictions on the reg¬ 
ulators R{k) such as monotony in k: R{k) < R{k') for 
k < k'. 

The basic building block of the flow operator AS '2 is 
the full propagator G = l/(r(^i -h R), and it would seem 
natural to use da- However, dclR^R'] does not qualify 
directly for measuring the distance: for physically close 
regularisations i?, R' the distance c?r( 2 ) [i?, R'] is necessar¬ 
ily small ®. Then, da[R, R'] is determined by the differ¬ 
ence (R — R') evaluated in the regularised regime which 
has no physical implication. Still, dc can be turned into 
a simple relation for the effective cut-off scale keS with 

dG.sup[i?,^] = \\G[R]\U^ = ^ktr , (5.21) 

with 

||G[i?]||,,p = sup{||G[((.,i?]|LJ, (5.22) 

where the supremum is taken in configuration space. The 
norm ||.||^2 i® tbe operator norm on L 2 already used for 
the criterion (5.10). In (5.21) dime is the momentum 
dimension of G, e.g. dime = —2 for bosons and dime = 
— 1 for fermions. is the wave function renormalisation 
of the field and makes the definition of fceff invariant 
under RG rescalings. In most cases the norm (5.22) will 
be evaluated in momentum space where it reads explicitly 

||G[i?]|L„p = sup \( [\G[cj,,R]i;\\p)) ^ ] . (5.23) 

,^,||V-||2 = 1 [\jp ) ] 

Note that the use of Z^ is not necessary as long as one 
uniquely fixes the endpoint of the flows, the theory at 
vanishing regulator. If one allows for simultaneous RG 
rescalings of the flow trajectories the prefactor in (5.21) 
arranges for an RG invariant fceff. For including rela¬ 
tive RG rescalings of trajectories the supremum in (5.21) 
also has to be taken over RG transformations. For most 
practical purposes these more general scenarios are not 
of interest. 

The expression k'^^^ relates to the biggest spectral 
value the propagator G[(f>,R] can achieve for all fields (p. 
Therefore fceff is the smallest relevant scale and hence is 


More precisely this applies to the distance dy^p( 2 ) j [i?,where 
the function |/(fr)| is bounded from above. 


24 



the effective cut-off. In the limit fc ^ 0 the effective cut¬ 
off scale fceff tends towards the smallest mass scale in the 
theory ^ . 

As an example we study a scalar theory with R 
in the leading order derivative expansion: rfc[(/)] = 
f + Vk[(j)]) ■ For regulators providing a momen¬ 

tum cut-off we can adjust fc as a physical cut-off scale by 
taking as a reference regulator the sharp cut-off 

i?sharp(/) = - k^) - 1) . (5.24) 

For i?sharp it is guaranteed that is the momentum scale 
below which (j) modes do not propagate. Inserting (5.24) 
in (5.21) with = 1, the effective cut-off scale is 

, (5.25) 

( 2 ) . . . ( 2 ) 
where '■ is the minimal value for Vt ■ Hence, in 

theories with a mass gap the effective cut-off scale fceff 

does not tend to zero but settles at the physical mass 

scale of the theory. In the present example Ae"ff(fc = 0) = 
( 2 ) 

^0 mill! the minimum of the second derivative of the full 
effective potential. Note that the full effective potential 
is convex and hence Iq ^in — 

3. Optimisation criterion 

The analysis of the previous two sections allows to put 
forward a general optimisation criterion in a closed form: 

DrJ[(^,R] _ =0, (5.26a) 

—^stab 

with 

{R^} = |i? with ||G[i?]||,„p = , (5.26b) 

where /[(/), i?] are correlation functions in the given order 
of the truncation. The norm ||.l|sup and the effective cut¬ 
off fceff have been introduced in (5.21). For the sake of 
completeness of the definition (5.26) we recall its prop¬ 
erties here: dime is the momentum dimension of G, and 
the effective cut-off kes is related to the biggest spectral 


^ In a regime with anomalous momentum scaling G oc 
one should rather define ||G[t?]||sup = ^ with di¬ 

mensionful Z^. 


value of the propagator k'^^^/Z^ . The norm in (5.26) 
is the supremum of the L 2 operator norm, 

||G[i?]|L„p = sup{||G[<^,i?]|LJ, (5.27) 

4 > 

see also (5.23). If the theory or the truncation scheme 
admits a natural expansion point (j)Q, the supremum in 
(5.27) might be substituted by evaluating the propagator 
at (j)Q, e.g. a configuration (j)Q for which the minimum of 
the effective potential is achieved. 

As shown in section VD 1, the constraint in (5.26) can 
be rewritten as the constraint of minimal action of AS' 2 , 
(5.16): 

(DR^^S2[(t>,R]) i[(t>,R] _ =0. (5.28) 

The criterion (5.26) is not bound to specific truncation 
schemes. The trivial starting point at i? = 00 is eval¬ 
uated for fceff(i? = 00 ) = 00 (assuming dg < 0), the 
end-point at i? = 0 represents the mass gap of the the¬ 
ory, kesiR = 0) = (lll/F^^^llsup)^/'^™'^. The monotone 
parameter kgS defines the effective cut-off scale and in¬ 
terpolates between the classical theory at keS = 00 and 
the full theory at fceff(O). If the theory undergoes a phase 
transition, in particular if it is first order, the monotony 
of kes{k) within truncations is at stake. If this hap¬ 
pens it hints at a truncation scheme that is not well- 
adapted. Nonetheless it can be dealt with in (5.26), it 
simply demands a more careful comparison of regulators 
at an effective cut-off scale defined by (5.27). Indeed, 
such pathologies can be avoided if restricting the space 
of regulators to those with monotony in k, R{k) < R{k') 
for k < k' which entails that regulators implement a true 
mode (scale) ordering. There are further secondary regu¬ 
larity constraints, but we do not want to overburden the 
criterion (5.26) with technicalities. 

The general form of the optimisation criterion (5.26) is 
achieved by substituting ||G||sup by a general norm do as 
defined in (5.19). For example, an interesting option can 
be found in [75]. In most cases the norm (5.27) applied 
to G supposedly is the natural choice: the propagator G 
is the key input in AS' 2 , any iterative truncation scheme 
involves powers of G and hence the importance of its 
supremum is enhanced within each iteration step 

Even in its form (5.28) the optimisation criterion (5.26) 
shows some dependence on the correlation function un- 


First investigations within LPA reveal the suggested equivalence 
of different choices for d(p, see also [70]. 
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der investigation. Bearing in mind the discussion about 
observable-independent optimisation we apply this idea 
to (5.28). First of all, AS 2 depends on F^^^ (and pos¬ 
sibly higher derivatives of F). The functional optimisa¬ 
tion implies (5.26) for these correlation functions which 
maximises the physics content of AS' 2 . Consequently the 
derivative in (5.28) is taken at 

=0. (5.29a) 

-^stab 

Eq. (5.29a) facilitates the evaluation of (5.28) as it only 
requires the evaluation of derivatives w.r.t. the explicit in¬ 
dependence. An optimisation for almost all relevant cor¬ 
relation functions I within a given truncation order im¬ 
plies the vanishing of the operator AS 2 [</>, i?] on the 
span of these I. Assuming that we can embed this span 
in a normed vector space V/ we arrive at a correlator- 
independent optimisation 

Dr^AS2[(^,R] _ =0, (5.29b) 

with (5.29a) with the operator norm ||.|| on V/. The 
optimisation (5.29) minimises the action of AS '2 on cor¬ 
relation functions I within a given truncation order. The 
representation (5.29) allows for a clear understanding of 
the result of the optimisation with the example of the 
two-point function. Eq. (5.29a) entails that for optimal 
regulators i?stab the spectrum of F^^^ at the effective cut¬ 
off scale kes is as close as possible (for the set of regulators 
R±{kes) ) to that of the full two-point function at fc = 0: 
the physics content of F*^^^ is optimised. It also implies a 
monotone evolution of the spectral values of F^^^ for opti¬ 
mal regulators. In case F^^^ has negative spectral values 
at the initial scale, e.g. a non-convex potential, the above 
investigations lead to one fc-independent spectral value, 
up to RG rescalings. 

The criterion (5.26), (5.29) can be rewritten as a sim¬ 
ple criterion on the full propagator and the full vertices. 
For its importance and for the sake of simplicity we con¬ 
centrate on the standard flow (3.72) with 

A52 = (Gi?G)bc4^4^ 

= i9t|r(2)(G - Go)bc^5|^ , (5.30) 

where Gq is an appropriate i?-independent normalisation, 
that leads to well-defined insertions for correlation func¬ 
tions I if applying (G — Gq)-^. In the presence of a 
mass gap a possible choice is e.g. Gq = G[<j), i? = 0]. The 
partial t-derivative at fixed F^^^ commutes with Dr^ at 


.Rstab- There, Dr^ = Dr^ |p( 2 ). Now we use that the sec¬ 
ond functional derivative /54)^54)'° does not vanish on 
almost all /. Therefore a vanishing norm (5.29b) implies 

||9t|r(=)7?Kjr(^)(G-Go)|U=«^^^^ =0. (5.31) 

The norm in (5.31) derives from the operator norm on 
Vj, and hence is related to the truncation scheme. A so¬ 
lution of \\Dr^ |p( 2 ) (G — Go) = 0 for all k implies 

a solution of (5.31). Consequently we search for extrema 
on the spectrum of the positive operator G. Now we 
use that the positive operator G vanishes identically for 
R = 00 and tends towards the full propagator G[4>,R = 0] 
with positive spectrum at vanishing regulator. Then with 
(5.29b) and (5.31) we conclude that optimal flows max¬ 
imise G at a given keS for all spectral values, with the 
constraint that —dtG > 0 is a positive operator. The lat¬ 
ter constraint guarantees that the maximisation is glob¬ 
ally valid for all k. We conclude that optimal flows are 
those where G[4i, R] is already as close as possible to the 
full propagator for a given cut-off scale keS- This crite¬ 
rion can be cast into the form 

de^(G)[Rsted,,0] = min(ie,^(G)[i?±,0], (5.32a) 

R± 

for all A € R"'’ with {i?±} as defined in (5.26), and 9x is 
defined via its action on eigenvectors IV’Ag) of G 

0x{GMxa) = [a + (Ag - A) 0{X - Ac)] IV’Ag) , (5.32b) 

with Heaviside step function 0{x)^. The operator 0x used 
in (5.32a) resolves the full spectral information of G. The 
criterion (5.32a) entails the constraint that G[4>,Rstah] 
takes the closest spectral values (according to the norm) 
to the full propagator G[((), 0] for all R G {R_l}, starting 
from the boundary condition G[4>, 00 ] = 0, or alterna¬ 
tively at G[4>, 0]. This implies a minimisation of the flow, 
as well as monotony of the spectral values of G in fc: 
G[^, 0] > G[^, i?]. These considerations enable us to re¬ 
formulate (5.32a) without relying on the full propagator 
G[4>, 0]. We are led to 

||0A(r(2)[R,tab] +Rstab)|| = min||0A(r(2)[Rj_] +R^)|j 

R± 

(5.32c) 

for all A G R"*". We remark that in deducing (5.32c) from 
(5.32a) we have again used —dtG > 0 and F^^^ [R = 0] > 


^ Ox is required to be a bounded operator. Hence for general norms 
used in d[R, R'] (5.32b) has to be modified, see e.g. section VIIIE. 
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0. If the distance d is defined with the L 2 -norm in the 
given order of the truncation, (5.32c) is also conveniently 
written as 

<^ex(G)[^stab,oo] = maxde,^(G)[i?±,oo]. (5.33) 

R± 

Note that in general (5.32c) can be written as (5.33) 
and some supplementary constraints depending on the 
norm used in (5.32c), see e.g. section VIIIE. For each 
norm these supplementary constraints are straightfor¬ 
wardly derived from (5.32a). 

Eq. (5.32) is a simple optimisation procedure indepen¬ 
dent of the correlation functions I under investigation. It 
already works without computations of full flow trajec¬ 
tories. In its form the criterion (5.26) has already been 
successfully applied to Landau gauge QCD [128, 129], 
see also section VIIIC. We emphasise again that the 
appropriate norm relates to the truncation used. The 
above analysis extends to general regulators. There, one 
also has to take into account the evolution of higher ver¬ 
tices r^”^. Their properties under -variations at i?stab 
derive from (5.29a) by taking field-derivatives. Spectral 
considerations are more involved but it can be shown that 
an optimisation for general regulators implies (5.32). 

We close the section with some comments concerning 
the generality of (5.26), the existence of solutions, and 
its connection to the criterion (5.10) 

the definition of the set R± in (5.26b) guarantees the 
existence of i?stab for a general expansion scheme: within 
any given truncation scheme the set of {i?±} is bounded 
by possibly smooth modifications of the sharp cut-off and 
the optimal cut-off (5.11) as functions on the spectrum 
of and for spectral values A(r(^^) < To¬ 

gether with positivity and monotony of the regulators R 
this proves the existence of a stable solution of (5.26), if 
neglecting the i?j_-variation of T), L Indeed such a proce¬ 
dure defines a further truncation scheme on top of that 
at hand. Note also that possibly one has to introduce a 
A-ordering: we search for a solution to (5.32a), (5.32c) 
for a given A on the sub-space of solutions to (5.32a), 
(5.32c) for A' < A. 

The argument above fails for generalisations of reg¬ 
ulator functions where the demand of positivity and 
monotony of the regulator are dropped. Still, for rea¬ 
sonable choices the set R± sweeps out basically the area 


For its connection to the PMS condition (5.9) we refer the reader 
to the discussion below (5.17). 


bounded by, possibly smooth modification, of the sharp 
cut-off and the optimal cut-off (5.11). However, it is not 
guaranteed anymore that the boundary curves are them¬ 
selves in R±. Therefore, a strict extremisation for all 
momenta (spectral values) as demanded in (5.26) might 
fail for generalisations of (5.26). More details will be 
provided elsewhere. 

Both criteria, (5.10) and (5.26), are based on the same 
key idea of global stability. In (5.10) the set of regulators 
is defined by normalising the regulators at some 
momentum. Then the inverse gap !|G[.R, </>o]||l 2 of the 
full propagator is minimised. In (5.26) the set of regu¬ 
lators {Ra_} is defined as those with the same maximal 
spectral value (inverse gap) ||G[i?]||sup s-od the action of 
the flow operator AS '2 is minimised. With (5.10) one is 
comparing regulators with different effective cut-off scales 
but, roughly speaking, close physics content. Then, op¬ 
timal regulators are those where this physics content is 
achieved for the biggest effective cut-off scale. In turn, 
with (5.26) we compare regulators leading to the same 
effective cut-off scale and single out those that lead to 
correlation functions as close as possible to those in the 
full theory. 


VI. APPLICATIONS TO FUNCTIONAL 
METHODS 

In this chapter we discuss immediate structural conse¬ 
quences of the setting developed so far. First of all this 
concerns the interrelation of functional methods like the 
general flows studied here, Dyson-Schwinger equations 
[I49-I57], stochastic quantisation [158-160], and the use 
of iVPI effective actions [161-177]. All these methods 
have met impressive success in the last decade, in partic¬ 
ular if it comes to physics where a perturbative treatment 
inherently fails. Here, we discuss structural similarities 
as well as functional relations between these approaches 
that open a path towards a combined use as well as non¬ 
trivial consistency checks of respective results. We also 
highlight the important aspect of practical renormalisa¬ 
tion schemes that can be derived from general flows for 
either DS equations or iVPI methods. However, given 
the scope of the present work we only outline the rele¬ 
vant points, leaving a more detailed analysis to future 
work. 
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A. Functional RG and DS equations 


as well as 


1. DSEs as integrated flows 


9i(5’“) + A52[^,i?]rfc’“ =0. (6.4b) 


Formally Dyson-Schwinger equations (2.14) are inte¬ 
grated flows. They constitute finite functional relations 
between renormalised Green functions as well as bare ver¬ 
tices. They have been successfully used for the descrip¬ 
tion of the infrared sector of QCD formulated in Landau 
gauge, initiated in [149, 150], for a review see [151]. This 
approach is also tightly linked to a similar analysis in 
stochastic quantisation [158-160]. 

More recently, these investigations have been extended 
to finite temperature QCD, e.g. [157] and the review 
[156]. The formal finiteness of the DS equations is more 
intricate if solving them within truncations [149-157]. 
Here, we discuss Dyson-Schwinger equations and their 
flow in the presence of a standard regulator coupled to 
the fundamental fields. This allows us to construct a gen¬ 
eral consistent BPHZ-type renormalisation of DS equa¬ 
tions from integrated flows being valid beyond perturba¬ 
tion theory. The extension of the results to the general 
setting is straightforward. 

Recall the DS operator / given in (2.14) with (f = (p, 
the source J coupled to the fundamental fields: /dse = 
J — |S. Inserting this into (3.14) leads to 

/SsE[<^,i?] = r’“[((.,R] - (5’“[0]) ^ 0 (6.1a) 

with 

^SsE[^,^,i?] = (6.1b) 

Note that in (6.1a) has to be read as a function 

of The flow of /dse is given by (3.60) and reads 

(dt + AS2[^,R])iusE = 0. ( 6 . 2 ) 

The first term in the DSE (6.1) Tfe/ already satisfies 
(6.2), see (3.71). This leaves us with the separate flow 

[dt + AS2[AR]^{S^-)=Q. (6.3) 


Eq. (6.4) highlights the aspect of the functional RG 
as a differential DSE. The use of the above identities 
(6.1),(6.2) and (6.3) is twofold. Firstly they allow us to 
relate DSEs and flow equations in similar truncations, 
hence providing non-trivial consistency checks for both 
approaches. Secondly they open a path towards a com¬ 
bined use of functional RGs and DSEs dwelling on the 
advantageous features of both. For example, an infrared 
analysis within both functional approaches usually pro¬ 
vides a set of possible solutions whose intersection is pos¬ 
sibly unique. In QCD this can be directly achieved by a 
fixed point analysis of (6.4a) along the lines in [128, 129]. 


2. Renormalisation 


Furthermore the flow equation in its integrated form 
can be used to set up an explicit renormalisation proce¬ 
dure within general truncation schemes. Such a renor¬ 
malisation is not necessarily multiplicative but gener¬ 
alises the BPHZ renormalisation of perturbation theory 
to general expansions. As it relies on a functional equa¬ 
tion for the effective action its consistency is guaranteed 
by construction. Hence it is possible to derive consistent 
subtraction schemes for Dyson-Schwinger equations from 
the integrated flow in a given truncation. 

We illustrate the above statements within the standard 
flow (3.75) for the effective action. Assume that we have 
solved the theory within the zth order of a given general 
truncation scheme, leading to Generally the flow 

can be written as 


dtV = k^'^Gab = (In G)aa " k^fCab ■ (6.5) 


In its integrated form this leads to 


rfc=rA-l (InG)aalA 



dttfG 


ab ’ 


( 6 . 6 ) 


Eq. (6.3) also follows directly from considering I = 
S'’“[^]. By construction the corresponding correlation 
function I satisfies the flow equation (3.60) and is given 
by I[(j),R] = From the above identities we also 

relate t-derivatives of Ffc’“ and i.e. 

atrfc^“ + A52[^,i?](^’“)=0, (6.4a) 


The integrated flow (6.6) represents an integral equa¬ 
tion for the effective action F^ with the boundary con¬ 
dition Fa. Note that its solution for a given k re¬ 
quires its solution for k' € [fc,A]. As such it consti¬ 
tutes a Dyson-Schwinger equation. It provides an ex¬ 
plicit (re)normalisation procedure involving two different 
aspects. Firstly the choice of a finite boundary condition 
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Fa implicitly renormalises the theory: it ensures finite¬ 
ness. The renormalisation conditions for the full effective 
action, i.e. fixing the relevant operators (of Fq) at some 
renormalisation scale /i translate to similar conditions for 
Ffe for all k. In particular its choice at A:' = 0 relates to 
an appropriate normalisation at fc = A. As can be seen 
from the representation of the integrated flow in ( 6 . 6 ) the 
renormalisation is done in a BPHZ-type way with sub¬ 
tractions — 5 lnG(fc' = A) -I- (Fa — S'), the t-integral also 
comprises some sub-leading subtractions. 

With ( 6 . 6 ) we have resolved the notorious consistency 
problem for explicit renormalisation procedures within 
Dyson-Schwinger equations. Practically it can be solved 
within an iteration of F^ about some zeroth iteration 
step Ffe^o for k G {A, 0}, e.g. F^^o = 5'ci, the classical 
action. This works for paths R{k), for which the initial 
condition Fa is sufficiently close to the classical action, an 
example being regulators R implementing a momentum 
regularisation with A setting a high momentum scale. 

An interesting option are non-trivial Ffc^o that already 
incorporate some non-trivial physics content of the the¬ 
ory under investigation. If the zeroth iteration step is al¬ 
ready close to the full solution the numerical effort is min¬ 
imised. Accordingly such a procedure benefits from any 
information already collected by other means about the 
physics content. In comparison to the standard (numer¬ 
ical) solution of DS-equations involving momentum inte¬ 
grations one has to perform an additional t-integration. 
In general this is bound to increase the numerical costs. 
However, this additional integral comes with the bene¬ 
fit that now the integrand is localised in momenta and 
t which stabilises the numerics. Indeed, the above ideas 
have been used for resolving the infrared sector of QCD 
within the Landau gauge thus furthering the evidence for 
the Kugo-Ojima/Gribov-Zwanziger confinement scenario 
in this gauge [128, 129], and providing a general con¬ 
sistent renormalisation procedure for related DS-studies 
[151, 152]. This aspect will be further discussed in sec¬ 
tion VIII. We also remark that the present analysis 
can be extended to the stochastic quantisation [158-160]. 
There it helps that we do not rely on an explicit path in¬ 
tegral representation. This shall be detailed elsewhere. 

Still the question arises whether ( 6 . 6 ) can be used more 
directly for setting up a renormalisation procedure for 
functional equations in the full theory at fc = 0 , solved 
iteratively within a given general truncation scheme 

F«[<^,i?] =Fi*-')[<^,i?] + AfoFfc[((.,i?], (6.7) 

as introduced in (5.1) for general R- Assume we have 


managed to construct regulators R that lead to a sup¬ 
pression of modes in the path integral related to orders 
i > ik of our truncation scheme. As an example we 
take the derivative expansion. Here we can use regulators 
that suppress at fc = A all momentum-dependent fields, 
lA = 0. By decreasing k we add more and more deriva¬ 
tives, ik —>■ 00 with k 0, either continuously switching 
on their effects or adding more and more derivatives in 
discrete steps. 

If R implements the truncation in discrete steps the 
flow only is non-zero at the discrete set of R. Integrating 
the flow from ki < ki < ki+i and ki+i < < ki +2 we 

arrive at 

rb+i) = Ffo-l((lnG)i),+i)-(lnG)fo) 

-i / 'dtf-“'’Ga6. (6.8) 

Jki 

Eq. ( 6 . 8 ) recursively implements the renormalisation at 
a given order i -|- 1 of the truncations by subtraction of 
appropriate terms of the order i. Naively the integral in 
( 6 . 8 ) can be performed as F’®^ only is non-zero at 
However, this has to be done carefully for similar rea¬ 
sons to those that do not allow for a naive integration 
of sharp-cut-off flows: at fcj+i, the flow F’^*' is singular 
and G jumps. Nonetheless, as in the case of the sharp 
cut-off ( 6 . 8 ) can be easily integrated within explicit iter¬ 
ation schemes. For example, perturbation theory within 
BPHZ-renormalisation can be reproduced with ( 6 . 8 ) but 
it extends to general schemes as well as general functional 
relations and correlation functions I of the theory that 
require explicit renormalisation if it comes to truncations. 


B. Composite operators and NPI flows 

The analysis of the last section extends naturally to 
flows in the presence of composite operators, in particular 
to flows of iVPI effective actions [161-163]. Flows with 
the coupling to composite operators have been considered 
in e.g. [21, 41, 77, 79-82]. Flows for the 2PI effective 
action have been studied in [77, 79, 82]. 

In the presence of sources for composite operators the 
renormalisation of these operators has to be taken into 
account. In particular, the construction of practical con¬ 
sistent renormalisation schemes within truncations poses 
a challenge, see e.g. [164-172]. Such a renormalisation 
has to respect the symmetry and symmetry breaking 
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pattern of the theory under investigation. We discuss 
the use of general flows for the construction of consistent 
subtraction schemes in general truncations by extending 
the renormalisation ideas of the last section. We also 
discuss the direct relation between flows in the presence 
of composite operators and A^PI effective actions, relying 
on the interpretation of the regulator i? as a source for a 
composite operator. 


1. Linear flows 


The structure of the flows (3.28),(3.60) always allows 
us to reduce the order of derivatives in ASk at the ex¬ 
pense of introducing further tensorial currents. In general 
we have 


5J“1' 
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'e-’ 

• • br 

‘^by ■ bn ^ 


A) 




xe 

jbl ■ t 

*n + i 



—Y 

. + m J 

7bl 

•bn 

0bi-.- 


xe 

jbi ■ 1 

’n + 1 


••bn^b^_|_j^--.b^_|_^ 


(6.9) 


with Sin+j = a'. Eq. (6.9) is valid for all i G IN. We also 
could have substituted only a part of the derivatives, ob¬ 
viously the relation is not unique. In case the source 
term ■'^"+”‘ 0 bi - b„ 0 b„+i - b„+„ was not present in 
the Schwinger functional W[J] it has to be added. Note 
that the derivatives w.r.t. t are taken at fixed arguments 
J and 0 respectively. Hence the reduction to lower pow¬ 
ers of derivatives is accompanied by holding the corre¬ 
sponding Green functions fixed. With (6.9) a part of 
the regulator term (3.2) with nth order derivatives, is re¬ 
duced to order n — m -I-1 by adding a further source term 
to W[J] 

J®(^a ^ = J® <^a' , (6.10) 


with 0 ai • ■ • ^a„ ■ The 

above relation is not unique, and we could have further 
reduced the order of derivatives by identifying additional 
products 4>eLi---sLm = 4>eLi ' ’' for n — m>m. By recur¬ 
sively using ( 6 . 10 ),( 6 . 12 ) with general m we can substi¬ 
tute AS” by an expression with only quadratic derivative 
terms, and the flow reduces to the standard form of the 
flow equation (3.72). Reducing AS one step further we 
arrive at first order derivatives w.r.t. J and (3.60) boils 
down to 


dtik[<l^]=0. (6.13) 

It seems that (6.13) is rather trivial but it should be read 
as a fixed point equation for the flow. When evaluating 
= 7 ’^b(T*’ — i?*^) = resulting from = y'^bT'^ 
the flow (6.13) reads 

=i?-, (6.14) 

where the partial t-derivatives is taken at fixed fields (jja. 
Eq. (6.14) yields upon integration 

atrfc[<^] = R>a, (6.15) 

which also can be read off from (3.64). If ■an _ q 
for n 7 ^ 2, (6.15) resembles the standard flow equation 
with G fo particular for (f) = (p. However, even 

for general n its integration is trivial: we exploit that for 
k = Q the regulator vanishes, i? = 0 and get 

rfc[<^] = ro[<^] + = r[<^] + A 5 fe[<^]. (6.16) 

Eq. (6.16) can directly be obtained by evaluating the 
Legendre transformation (3.53) for the present scenario. 
For regulator terms linear in (j), ASk[4>\ = R’^4’a, there 
is a simple relation between the Schwinger functional 
of the full theory and that of the regularised theory: 
lEfe[J] = Wo[J — R]. Moreover = 0. With these 

observations we can rewrite (3.53) for linear ASk as 


where 

J-''-' , ( 6 . 11 ) 


with enlarged multi-indices a' = a, ai • • • a^ and 7 = 
(y'^'^). Eq. (6.10) implies = 0ai---^a„,- With 

( 6 . 10 ) we are led to 





m + 1 
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gja'n-m + l 




'((3.12) 


Efc = sup(J"<^a - VE[J - i?]) 

J 

= sup((j - - W[J - R]) + ASk 

J 

= r-bAS'fc. (6.17) 

In (6.17) we have used that the supremum over the space 
of functions J is the same as that over the space of func¬ 
tions J — R. Strictly speaking, the last equality in (6.17) 
is only valid for the subset of regulators R that can be 
absorbed in currents J. 
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From the above definitions and the flow (6.15) we can 
step by step resolve the composite operators by us¬ 
ing the related equations of motion. Here we show how 
such a procedure can be used to finally recover the reg¬ 
ularised effective action rfc[(/)o] in (3.53) and the general 
flows (3.60). The equations of motion for (/)oi...o„. for 
rii >2 read 

Vni>2. (6.18) 

d(Pai ■ ■ -Onj 

Using the solution ^{(pa) = {(pa, of 

(6.18) in (6.16), we end up with the effective action 
(3.53). As AS'). = 0 for linear regulators we have 

^k[(pa] = ^k[^] - ASk[(p{(pa)], (6.19) 

where ASk[^{(pa)] = Z)* - Due to the 

linearity of the t-derivative the flow (6.15) holds true also 
for the effective action rfc[())a]. This statement reads more 
explicitly 

dtTk[(Pa] = dt\^rS] + r^[^]dtM^a] 

= dt\^rk[4>]. ( 6 . 20 ) 

The second term on the rhs of the first line in (6.20) van¬ 
ishes due to the equations of motion (6.18) for m > 2 and 
due to dt(pa[(pa] = 0 for the fundamental field (pa '■= (pa, 
that is not a solution to the related equations of motion 
but a general field. Hence the flow equation for the IPI 
effective action reads 

dtrk[(pa] = ( 6 . 21 ) 

The equations of motion (6.18) relate the fields (pa[(pa] to 
a combination of Green functions 

^a[0a] = (0a[<^a])Ja = (J<»,0) • (6.22) 

The relations (6.22) can be written in terms of functional 
(((-derivatives as 

^a= (0a[Gafc4 +((.,]) . (6.23) 

As an example we use (6.23) for the two-point function 
(^aias = 0 ai(^a2 and {(pa_) = {(pa, (pa^a.^) ■ It followS 


(3.60) starting at the trivial flow in (6.13), dtik = 0. The 
flow for ik[(pa] = Ik[H(pa)] reads 

dtik[(pa]-h"[4']dt^a = 0, (6.25) 

similarly to (6.20). In (6.25) we have used (6.13), there is 
no explicit t-dependence. In contradistinction to (6.20) 
the remaining term on the rhs of (6.25) does not vanish as 
general correlation functions do not satisfy the equations 
of motion (6.18). Note also that the fields (p trivially 
satisfy the flows (6.25). The fields (p{(pa) belong to the 
correlation functions Ik and hence they obey the flow 
equation 

dt<Pa[(pa] + AS2[(pa, I?] ^a[(pa] = 0 . (6.26) 

Inserting (6.26) into (6.25) we arrive at the flow 

dJklP^a] + {AS2^a)i’k[^] = 0 , (6.27) 

which implies (3.60). The latter statement follows only 
after some algebra from (6.27). For its proof one has 
to consider that AS '2 acts linearly on Ik which it does 
not on general correlation functions Ok However, 
it is more convenient to work with the flow (3.28) for 
and with the definition = (J“, 0 )]. 

By using the equivalence of J-derivatives (6.9) valid for 
the Ik, the flow for derives from that of 

as {dt + AS'i[J“,i?])/fc[J“] = 0, implying the flow (3.60) 
for I [(pa]- It is worth noting that truncated flows derived 
from either the representation (3.60) or (6.27) differ. This 
fact can be used for consistency checks of truncations as 
well as an improvement in case one of the representations 
is better suited within a given truncation. 

Accordingly there is a close link between API formula¬ 
tions of the effective action and general flows. Moreover, 
it is possible to switch back and forth between these for¬ 
mulations, thereby combining their specific advantages. 

2. 2PI flows 

As an explicit example we study the standard flow re¬ 
lated to the quadratic regulator term 


(p - (^{G-i^+(p)a^{G-^+(p)aJ^^ [ A] = ^ ^ , (6.28) 

= {(pa , Gaia2 + <pai(pa2) ■ (6.24) _ 


Inserting (6.23) into the flow (6.21) we recover the flow 
(3.63). The relation (6.23) also leads to the general flows 
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The proof can be worked out for A^-point functions (6.23) from 
where it extends straightforwardly. 



which can be linearised in terms of 2PI quantities 

4’aia2 ~ 0ai<^a2 ) (6.29) 

where (j)a is not necessarily a fundamental held. For (jia^a^ 
as dehned in (6.29) the relation (6.9) reads 

(5J“i (5J“2 

= , (6.30) 

^jaiaa ^ > 

Using (6.30) we reduce (6.28) to a linear regulator at 
the expense of also keeping the corresponding 2 -point 
functions hxed, 


dt4’aia2 — dt{G + 4>ai4’a2) — ^tG — 0 . (6.31) 

We substitute AS'^ in (3.13),(3.14) with 

r r r 

(5J“i (5J“2 (5 J“i“2 ^ ^ 


(6.32) 


and are lead to (6.13), = 0. The effective action 

and its flow are functions of the held (pa and the two-point 
function (pab'- 

rfe[<^a]=r[((.a]+i?“Va 6 , (6.33) 

with {(p^) = ipai.fpaia^) and 


9trfe[().a] = . 


(6.34) 


The how (6.34) resembles the standard how equation 
(3.60) and follows directly from the dehnition of F^ in 
(6.33). It also follows by integration w.r.t. p from (6.14) 


with -^Al = Ra’b' and = 0. The equation of 


motion in pab according to (6.18) is given by 

Ppab 0=0 

Its solution (6.23) reads 0a = [PaiPaia^) with 
pab — Gab T PaPb • 


(6.35) 


(6.36) 


The above relations lead to the standard how equation 
for the IPI effective action rfc[0o] = Tk[pa,pab]—R^‘^pbpc 
dehned in (6.19). With (6.35) it follows that [77, 79, 82] 


dtTkiPiP)] = dt\^n[p]+G^[p]dtP.[p] 


= dt\^rk[p,p]. (6.37) 

Using the how (6.34) in (6.37) we arrive at 

dtrk[Pa]=R'’"Gbc, (6.38) 

the standard how (3.74). Hence linear hows of 2PI quan¬ 
tities and its hxed point equations rehect the standard 
how equation and offer the possibility of using 2PI ex¬ 
pansions as well as results in standard hows. 


3. Renormalisation 

The setting in the present work hinges on the bootstrap 
idea that the path integral, more precisely the Schwinger 
functional 1U[J, i?], is hnite and uniquely dehned. Re¬ 
sorting to Weinberg’s idea of non-perturbative renormal- 
isability [30] this simply implies the existence of a hnite 
number of relevant operators in the theory. If not only 
the fundamental helds p = p are coupled to the path 
integral but also general composite operators p’^ some 
care is needed. As an example let us consider 0 ‘^-theory 
in d = 4 dimensions in the presence of a source for 
p^{x). More generally we deal with a Schwinger func¬ 
tional 1 U[J, i?] with J^pa = ' Pa&- 

The composite p^{x) operator is coupled with the choice 
jai - ae^^^ • • • ^ae = Ae J^p^(x}. However, at face value 
we have changed the theory to a ( 0 ®(a;)-theory with cou¬ 
pling Ae that is not perturbatively renormalisable in 
d = 4. Still, within functional RG methods one can 
address the question whether such the theory is con¬ 
sistent. In particular if the theory admits a non-trivial 
ultraviolet fixed point the problem of perturbative non- 
renormalisability is cured. Leaving aside the problem of 
its UV-completion the flow equation can be used to gen¬ 
erate the IR-effective action from some finite initial con¬ 
dition. Then, the flow equation introduces a consistent 
BPHZ-type renormalisation. 

In turn, as long as the composite operator p’^ is 
renormalisable we deal with the standard renormalisa¬ 
tion of composite operators [182]. Moreover, functional 
RG flows can be used to actually define finite generat¬ 
ing functionals in the presence of composite operators 
as well as practical iterative renormalisation procedures 
[41, 77]. The general case is covered by the RG equations 
(4. 8 ),(4.20) and the full flows (4.20). In particular we deal 
with a matrix 70*^0 of anomalous dimensions, and the 
corresponding renormalisation conditions, for the general 
perturbative setting see e.g. [182]. We resort again to the 
above example of in d = 4 but coupled to the 

2 -point function: J^pa = J°'pa + We have 

extended the number of (independent) relevant opera¬ 
tors {p‘^{x)), {{dp)^{x)) and {p‘^{x)) with {p{x,x)) and 
{p{x,x)p'^{x)) and {p‘^{x,x)), where p{x,y) = p{x)p{y). 
The anomalous dimensions of these operators are related 
by the matrix 70 and coincide naturally on the equations 
of motions. 

Apart from these more formal questions there is the 
important issue of practical renormalisation, i.e. consis- 


32 



tently renormalising the theory order by order within a 
given truncation scheme. The general flows (3.60) to¬ 
gether with the considerations of this section allow to 
construct such a renormalisation. Again we outline the 
setting within the 2PI effective action with a = a, 0102 
and 0a = ( 0 o) 0 ai 0 a 2 )- As distinguished from the last 
section VIB 2 we couple a quadratic regulator to the 
fields, 

A,S[0,i?] =i?«‘’0a0b, (6.39) 

where we also allow for insertions of the operators 0 a 0 hi &2 
and 001020f)ib2- The regulator (6.39) leads to the stan¬ 
dard flow (3.72) for general correlation functions, for the 
effective action it is given by (3.74). In the present case 
it reads 

dtTk[<l>] = 7?“*'Gob + i?“'’^'’=Gabib2 +^“^“^"^01026 

■ (6.40) 

In the first term on the rhs of (6.40) we could also identify 
Gab = 4>ab — 4'a4'b, See (6.24). 2PI expansions relate to 
loop (coupling) expansions in the field 0a and hence, via 
the equations of motion, to resummations of classes of 
diagrams. For general expansion schemes we refer to the 
results of section VIA 2 that straightforwardly translate 
to the present multi-index situation. 

We proceed by discussing an iterative loop-wise resolu¬ 
tion of the flow (6.40) that leads to a BPHZ-type renor¬ 
malisation of diagrams as in the standard case. This anal¬ 
ysis is not bound to the 2PI example considered above 
as the index a could comprise higher 7V-point functions. 
From now on we consider the general case. Still we keep 
the simple quadratic regulator (6.39). Assume that we 
have resolved the theory at ith loop order leading to a fi- 
nite f-loop contribution F),', the full effective action being 
Ffc = Then, the i-|-lst order reads in differential 

form 

, (6.41) 

and is finite. At one loop, i = 1, its integration results in 

ri'^[</']=-|(lnG)\|^ + rW, (6.42) 

where the A-dependent terms arrange for a BPHZ-type 
renormalisation procedure and, in a slight abuse of nota¬ 
tion, G stands for the classical propagators of the fields 
0a. The superscript indicates the one loop order, not 
the one point function. The subtraction at A makes the 


rhs finite. Fa ensures the A-independence as well as in¬ 
troducing a finite (re)-normalisation. For z = 2 we have 
to feed F*^^^ [0] and its derivatives into the rhs of the flow 

(6.41) . Again the t-integration can be performed as the 
rhs is a total derivative w.r.t. t. It is the same recur¬ 
sive structure which reproduces renormalised perturba¬ 
tion theory from a loop-wise integration of the IPI flow. 
At two loop the flow (6.41) reads 

^ r(l)’"‘lGdb , (6.43) 

assuming no coupling dependence of R. The two-point 
function at one loop, is the second derivative of 

(6.42) w.r.t. the field 0a, and (6.43) turns into a total t- 
derivative. Finally we arrive at the two-loop contribution 

p(2) ^ Ip,a,a2a3a.(^ _ G|A)a,a2(G - G|A)a3a. 

_^p,a3a2a3r.a3a3a2(G _ G|A)a3a2(G - G|A)a3a4 
x((G-G|A) + 3G|A)a3ae 

+ ^ri2ia2(G - G|A)a3a2 + . (6.44) 

Higher orders follow similarly. Such a procedure allows 
for a constructive renormalisation of the theory under 
investigation, and also facilitates formal considerations 
concerning the renormalisation of general truncations 
schemes. The first two terms in (6.44) are already fi¬ 
nite due to the subtractions. The terms proportional to 
3G in the third line of (6.44) and in the 4th line con¬ 
stitute finite (re-) normalisations. Eq. (6.44) stays finite 
if the vertices and propagators are taken to be full ver¬ 
tices and propagators in the sense of an RG improve¬ 
ment. Within the 2PI example considered in (6.40) the 
integrated flow (6.44) is the consistently renormalised re¬ 
sult for the 2PI effective action at two loop. It translates 
into a resummed renormalised IPI effective action by us¬ 
ing the equation of motion (6.18) for the composite field 
0ah. However, the above result also applies to VPI effec¬ 
tive actions or more general composite operators coupled 
to the theory: the integrated flow (6.44) constitutes a 
finite BPHZ-type renormalised perturbative expansion. 
Moreover, the above method straightforwardly extends 
to general expansion schemes: in general the integrated 
flow constitutes a finite BPHZ-type renormalised expan¬ 
sion. The consistency of the renormalisation procedure 
is guaranteed by construction. 

The renormalisation conditions for the full theory are 
set implicitly with the choice of the effective action at 
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the initial cut-ofF scale A. We emphasise that any RG 
scheme that derives from a functional truncation to the 
flow (3.60), and in the particular the loop expansion 
(6.41), is consistent with the truncation. Moreover, the 
iterative structure displayed in (6.41), (6.42) and (6.44) 
allows us to discuss general renormalisation conditions in 
the present setting. By adding the operator (j)ab we have 
extended the number of relevant vertices in the effective 
action and hence the number of renormalisation condi¬ 
tions. In case 4>a. includes only marginal and irrelevant 
operators the renormalisation proof can be mapped to 
that of the IPI case. 

The basic example is provided by {(pa) = {(paiPaia^)^ 
where the field (paia^ with = PaiPa-i) counts like 

paipa 2 - RG conditions for e.g. the 2-point function and 
the 4-point function 


(5r (5r 

dpaSpb' Spat-'-Spai" 


(6.45) 


trigger additional RG conditions for 


(5r (5r 

^Pab ^Paia2^Pa2a4 


ST 

Spa\a2^Pa2^Pa4 


.(6.46) 


Using the relation (6.9) between derivatives w.r.t. pa and 
pab we are left with the same number of independent RG 
conditions as in the IPI case. In other words, the matrix 
is highly symmetric. This symmetry can be imposed 
on the level of Pa and evolves with the flow as its rhs only 
depends on (derivatives of) Pfc. We observe that formally 
any choice of Pa independently fixes these RG conditions 
at all scales (via the flow) but violates the relation (6.9). 
A priori there is nothing wrong with such a procedure 
that simply relates to an additional additive renormali¬ 
sation (at IPI level) and can be absorbed in a possibly 
/c-dependent rescaling of the 2PI fields. The above dis¬ 
cussion extends to the general case with fields pa- We 
shall detail these observations and structures elsewhere 
and close with the remark that for general truncation 
schemes that do not admit a direct resolution of the flow 
as in perturbation theory, the costs relate to an addi¬ 
tional t-integration as already discussed in the IPI case 
of section VI A. 


theories [98-134] and gravity [I42-I47] with the stan¬ 
dard quadratic regulator one has to deal with modified 
Slavnov-Taylor identities [98-115]. These identities tend 
towards the Slavnov-Taylor identities of the full theory in 
the limit of vanishing regulator. It is crucial to guarantee 
this limit towards physical gauge invariance. 

The subtlety of modified Slavnov-Taylor identities can 
be avoided for thermal flows. This is achieved by ei¬ 
ther modifying the thermal distribution [121, 122], or by 
constructing the thermal flow as a difference of Gallan- 
Symanzik flows at zero and finite temperature in an axial- 
type gauge [17]. The resulting thermal flows are gauge 
invariant. We remark that Gallan-Symanzik flows in ax¬ 
ial gauges at zero temperature [116-119] are formally 
gauge invariant, but the approach towards the full the¬ 
ory at vanishing regulator has severe consistency prob¬ 
lems. This problem is related to the missing locality in 
momentum space combined with the incomplete gauge 
fixing [112]. One expects a better convergence for Gallan- 
Symanzik flows within covariant or Abelian gauges [120] 

Alternatively one can resort to gauge-invariant degrees 
of freedom [140, 141], gauge-covariant degrees of freedom 
[135-139], or higher order regulator terms with regulators 
jjai -an n > 2. Then, V-point functions directly 
relate to observables and allow for the construction of 
gauge-invariant flows. In general such a parameterisation 
is payed for with non-localities, in particular in theories 
with a non-Abelian gauge symmetry. 

In this chapter we discuss the structural aspects of the 
above formulations. In particular we deal with the ques¬ 
tion of convenient representations of symmetry identi¬ 
ties that facilitates their implementation during the flow. 
Moreover we discuss the related question of adjusted pa- 
rameterisations of gauge theories, and evaluate the fate 
of symmetry constraints in gauge-invariant formulations. 


A. Parameterisation 


VII. APPLICATIONS TO GAUGE THEORIES 

The generality of the present approach fully pays off 
in gauge theories, and the present work was mainly trig¬ 
gered by related investigations. In flow studies for gauge 


In gauge fixed formulations of gauge theories, and in 
particular in strongly interacting regimes, the propaga¬ 
tors and general Green functions are only indirectly re¬ 
lated to physical observables. Firstly, only combinations 
of them are gauge invariant and secondly, the relevant 
degrees of freedom in the strongly interacting regime are 
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not the perturbative ones Good choices are observ¬ 
ables that serve as order parameters; e.g. the Polyakov 
loop^^ 

f0 

P(a;) = TrT^exp / Ao(x)dT , (7-1) 

Jo 

and its two-point function {P{x)P^ {y}) in the case of the 
confinement-deconfinement phase transition. These ob¬ 
servables fall into the class of Ik defined in (3.14). For 
the Polyakov loop variable (7.1) the corresponding op¬ 
erator is / = P{x)[Aq = which implies Ik = I, 
see (3.14b). Hence their flow can still be described in 
terms of field propagators and vertices of the fundamen¬ 
tal fields via (3.14),(3.51). It amounts to the following 
procedure: compute the flow of propagators and vertices, 
even though partially decoupling in the phase transition. 
Then, the flow of relevant observables / is computed 
from this input with the flow (3.60), i.e. the heavy quark 
potential from the flow of the Wilson loop or Polyakov 
loop. Such a procedure allows for a direct computation 
of physical quantities from the propagators and vertices 
of the theory in a given parameterisation, and it applies 
to gauge fixed as well as gauge invariant formulations. It 
also emphasises the key role played by the propagators 
of the theory, and matches their key importance within 
the functional optimisation developed in section V. 

One also can use appropriate fields (j) coupled to 
the theory. In the above example of the confinement- 
deconfinement phase transition a natural choice is pro¬ 
vided by the gauge invariant field 4>(x) = P{x) with 
(7.1). Such a choice has to be completed by additional 
<jp that cover the remaining field degrees of freedom. Al¬ 
ternatively one can integrate out the remaining degrees 
of freedom and only keep that of interest. Another in¬ 
teresting option are gauge covariant degrees of freedom, 
e.g. (j)^v{x) = F^i, or = Ff_cu, that is the dual field 

strength [123, 124]. Both choices can be used to derive 
(partially) gauge invariant effective actions, and aim at 
a description of gauge theories in terms of physical vari¬ 
ables. 

We emphasise that the above suggestions usually gen¬ 
erate non-local and non-polynomial effective actions even 
at the initial scale. We have to keep in mind that gauge 


Basically by definition; the relevant degrees of freedom should 
only weakly interact. 

The definition (7.1) only applies in the case of periodic boundary 
conditions for the gauge field. 


theories are formulated as path integrals over the gauge 
field supplemented with a polynomial and local classi¬ 
cal action. Gauge fixing is nothing but the necessity to 
deal with a non-trivial Jacobian that arises from the de¬ 
coupling of redundant degrees of freedom, and Slavnov- 
Taylor identities (STIs) carry the information of this 
reparameterisation. If coupling gauge invariant or gauge 
covariant degrees of freedom to the theory the necessity 
of decoupling the redundant degrees of freedom remains, 
and hence the symmetry constraints are still present. In a 
gauge invariant setting the corresponding STIs turn into 
a subset of DSEs. Their relevance might be hidden by the 
fact of manifest gauge invariance, but still they carry the 
information about locality. In other words, approxima¬ 
tions to gauge invariant effective actions or general corre¬ 
lation functions still can be in conflict with the Slavnov- 
Taylor identities and hence violate physical gauge invari¬ 
ance. Indeed it is helpful to explicitly gauge fix the theory 
within a choice that simplifies the relation (j) = 4){A) for 
gauge-fixed fields A as it makes locality more evident in 
the variables 4>. For example, in case of the confinement- 
deconfinement phase transition we choose 4>{x) = P{x) 
defined in (7.1), and use the Polyakov gauge or variations 
thereof, e.g. [178-180]. 

In summary we conclude that it is vital to study 
the fate of symmetry constraints such as the Slavnov- 
Taylor identities for general flows, be they gauge invari¬ 
ant or gauge variant. This is done in the next three 
sections VIIB, VBG,VBD. 

B. Modified Slavnov-Taylor identities 

The propagators and vertices of a gauge theory are con¬ 
strained by gauge invariance of the theory. A non-trivial 
symmetry /^ = 0 is maintained during general flows 
(4.8), (4.20): if /^ = 0 is satisfied at the starting scale, its 
flow vanishes as it is proportional to Ik ■ In particular this 
is valid for Dg = dt ■ The corresponding flows include that 
of modified Ward-Takahashi or Slavnov-Taylor identities 
for the effective action [104, 110, 112], and that of Nielsen 
identities [141] for gauge invariant flows [140, 141]. 

The above statements imply that the generator of the 
flow, Ug, commutes with the generator of the modified 
symmetry Ik- Within truncations this property does not 
hold, and it is not sufficient to guarantee the symmetry 
at the starting scale. Consequently a symmetry relation 
Ik = 0 should be read as a fine-tuning condition which 
has to be solved at each scale. This is technically rather 
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involved, and any simplification is helpful. Here we aim 
at a discussion of different representations of symmetry 
constraints and their flows. 


Matter fields and a Higgs sector can be straightforwardly 
added. The action (7.5) is invariant under the BRST 
transformations 


1. STI 

First we concentrate on a pure non-Abelian gauge the¬ 
ory with general gauge fixing iF[A]. For its chief impor¬ 
tance we shall explain the structure with sources coupled 
to the fundamental fields and a standard quadratic 
regulator term R°-^ipa(pb- We keep the condensed nota¬ 
tion and refer the reader to [181] for some more details. 
The Schwinger functional is given by 

^W[J,Q] ^ j _ (7 2) 

In (7.2) we have also included source terms for 

the symmetry variations of the fields as introduced in 
section H. Here s generates BRST transformations de¬ 
fined below in (7.8). The fields (l)a\0\ depend on the fun¬ 
damental fields ipa given by 

{<Pa) = {A,, , C„), (7.3) 

where we have dropped the hats on the component fields. 
The component fields in (7.3) read more explicitly Ai = 
A^{x), the gauge field, and Ca = Ca{x), Ca = Ca{x), 
the ghost fields. A more explicit form of the source term 
in the case of </> = (^ reads 


{5ip) = {D^C^ , , A„), (7.8) 

and s acts trivially on A: sAq- = 0. The operator s can be 
represented as a functional differential operator on the 
fields (p, A with 

S = {B(pa) , (7.9) 

making the anti-commuting (Grassmann) property of s 
explicit. The invariance of the action, sS'[(,J] = 0 can be 
proven straightforwardly by insertion. Moreover, s is a 
differential with 5^Lp = 0 allowing for a simple form of 
the symmetry constraint. The only BRST-variant term 
is the source term The related Slavnov-Taylor 

identity (STI) is cast into an algebraic form with help of 
the source terms for the BRST variations (7.8) included 
in (7.2). For (j> = p> this source terms reads 

+ Q“A„, (7.10) 

where Q" sCq, = Q°‘ Xa could also be considered as a 
standard source term for the auxiliary field A. The gen¬ 
eral BRST source term reads 

Q"(4)a = Q"(5^)a4“[<^], (7.11) 


a — Ji-^i H” JocCqc Jqc^oc 


(7.4) 


= ^ (J“(x)A“(x) + J“(x)C“(x) + C“(x)J“(x)) . 

The action S in the path integral (7.2) is given by 
A] = S'ym[i^] 


-w(A) + A„.A„(A) - , (7.5) 

with 

(A), Lo{X) = |a„A„ , (7.6) 

the latter equation for to leading to the standard gauge 
fixing term upon integration over A. Then, in a 

less condensed notation, (7.5) turns into 




ca 


dAt 


(7.7) 


following with (7.9). The Slavnov-Taylor identity follows 
from 

J s [dipdX exp{-S'[(^] -|- J®(^a + g“(s<^)a}) = 0 . (7.12) 

Eq. (7.12) is of the form (2.12). It follows with (7.9) after 
a partial functional integration and = Dffia + 

faapCp = 0 (for compact Lie groups). Except for the 
source term J‘^0a all terms in (7.12) are BRST-invariant: 
sd(p = 0, sS'[(p] = 0, s(g‘^S(^a) = 0. The operator s 
commutes due to its Grassmannian nature with bosonic 
currents J and anti-commutes with fermionic ones. For 
example, for the fundamental fields and currents this en¬ 
tails that s commutes with J* but anti-commutes with 
J“, J“ and 5J°'(pa = J^7“b(si^a). Using all these proper¬ 
ties in (7.12) leads us to the Slavnov-Taylor identity 

J dip dXJ^^‘^b{5ipa) exp{-S'[l^] -I- J®i^a + Q^^Sipa} 

= = 0 . (7.13) 
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Eq. (7.13) is of the form I[J,Q] = 0 leading to (3.7) 
with / defined in (2.10) for 

= ( 7 . 14 ) 

The operator 1^ generates BRST transformations on the 
Schwinger functional W. Accordingly the STI (7.13) can 
be written as 

/,1E[J,Q] = 0, (7.15) 

that is the Schwinger functional is invariant under BRST 
transformations. The STI (7.15) can be generalised to 
that for correlation functions I. To that end we use that 
(7.15) can be multiplied by any operator / from the left. 
We are led to 

= with Ws,i=ih, (7.16a) 

where W/ is derived from W/ with (2.10). The symmetry 
relation (7.16) is a direct consequence of (7.13), which 
is reproduced for 7=1. We can write the correlation 
function W/ in terms of I as 

W,,i[J, Q] = h I[J, Q] + Shj[J, Q], (7.16b) 

with 

(7.16c) 

For the derivation of (7.16b) we have used that I Is = 
is I + 61 as well repeatedly using [7^ , VF] = 0, which is 
the STI (7.13). ^ 

For Q-independent 7 the commutator SI substitutes 
one of the J-derivatives in 7 by one w.r.t. Q. Applied on 
this generates a (quantum) BRST transformation on 
7. Consequently we write 

TlsAJ,^]e'^ = - {miiJA]) (7.17) 

which we evaluate a.t (f> = A- Accordingly, for BRST- 
invariant 7[J, </>] the second term on the rhs of (7.16b) 
disappears. Hence, if 7 is the expectation value of a 
BRST-invariant 7[J, </>], the second term on the rhs of 
(7.16b) vanishes and 7 is BRST-invariant, 7^ 7 = 0. 

We remark that (7.16), as the flow (3.28), does not di¬ 
rectly encode the STI for the Schwinger functional. This 
comes about since we have divided out the STI for W, 
(7.15) in its form [is , W] in the derivation of (7.16). In 
turn, it has to be trivially satisfied. Indeed, for either 
7 = I or 7 = W[J,Q], leading to 7 = 1 and 7 = W, the 
STI (7.16) is trivially satisfied. The situation is similar 


to that of the flow (3.28) where the flow of the Schwinger 
functional has been divided out. Without resorting to 
the STI for W, (7.15), the STIs Wsj derived with VVsj 
in (7.16a) read 

WsAJ, Q] = [is - {isW)) 7[J, Q] +61, (7.18) 

and, for 7 = 1 or 7 = IF[J, Q] the STI for the Schwinger 
functional, (7.15) follows. Hence, we shall refer to the 
STI (7.15) as Ws,i = 0. Note also that its trivial resolu¬ 
tion does not imply that it is not encoded in the repre¬ 
sentation (7.16b). Similarly to the derivation of its flow 
from the general flow (3.28), the STI for the Schwinger 
functional derives from 7 = ^, inserted in (7.16). We 
are led to ^75lF[J, Q] = 0 which entails (7.15). 

g. mSTI 

So far we have adapted the analysis of the STI in its 
algebraic form to the present setting. Now we consider 
regularisations of the Schwinger functional W[J,Q, R] de¬ 
fined in (3.1), as well as general operators I[J,Q,R] de¬ 
fined in (3.8). The operator ^,^,7?] correspond¬ 
ing to Is[J, Q, R] is derived from (3.8b) as 

is = {J^-[AS,A]h\s^, (7.19) 

where the second term generates BRST transformations 
of the regulator term AS', and we have used that AS 
is bosonic. As an example we compute (7.19) for the 
standard flow, (j) = ip and a quadratic regulator term 
R°‘^(pa’Pb- This leads us to the symmetry operator 

is = {A - 2R^^A_A\j^ , (7.20) 

where we have used the symmetry properties of R in (3.5) 
for standard flows. The STI for the Schwinger functional 
(7.13) turns into [98-115] 

isW[J,Q,R] = {), (7.21) 

with is defined in (7.20). It entails that only the source 
terms J^<j)a. and the regulator term are BRST-variant. 
The relation (7.21) was coined modified Slavnov-Taylor 
identity (mSTI) as it encodes BRST invariance at 7? = 0, 
and shows its explicit breaking via the regulator term at 
7? 7^ 0. 

The general case with W/ leads to the same general 
STI (7.16) with all operators and correlation functions 
substituted by their 7?-dependent counterparts defined 
in (3.8), 

W,./[J,Q,77] = 0, (7.22a) 
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with 


with 


Q, R] = 4 I[J, Q, R] + SI. (7.22b) 

The correlation function SI[J,Q,R] is the i?-dependent 
counterpart derived from (7.16c) with (3.8); 

SI[R] = [i, h]n=o , (7.23) 

where AS = AS'[^,i?]. Hence the second term on 
the rhs of (7.22b) still vanishes for a BRST-invariant 
I[J, 0]. The modification of BRST invariance is solely 
encoded in the modification of the BRST operator Ig in 

(7.20) . The flow of (7.22) is governed by (3.28). 

The mSTI (7.21) for the Schwinger functional follows 
as yVs,i = 0 with the alternative representation (7.18). 
As in the case without regulator, it also can be derived 
from (7.22) from Ws^w,^- Inserting I = -A into (7.22b) 
leads to Ww,a = -^hW[J,Q,R\ = 0 and hence to 

(7.21) . 

As in the case of the flows we can turn the general 
mSTIs (7.22) into mSTIs for correlation functions / in 
terms of the variable (j). The definition of the effective 
action (3.43) extends to the case with external currents 
Q: 


T[<P, Q, R] = Q)<l)g, - W[J{<P, Q), Q, R] 

-AS'[<I>,R], (7.24) 


the source J now depends on the fields (j) and the source 
Q. Eq. (7.24) entails that 

^IT <f(r + AS’) 5T 

as AS' does not depend on Q and the Q-dependence of 
J cancels out. In (7.25) the Q-derivatives of W and T 
are taken at fixed arguments J and (j) respectively. The 
correlation functions I derive from (3.51) as 


i[^,Q,R]=I[Ji^,Q,R),Q,R]. (7.26) 


For the mSTIs Wj = 0 we have to rewrite Q"derivatives 
at fixed J in terms of Q"derivatives at fixed <j). This reads 


5 

SQg 



6 

SQ 


■ 7 b 






6 

6(j)c ’ 


(7.27) 


where we have used (3.45) and (3.50). With the above re¬ 
lations we arrive at the modified Slavnov-Taylor identity 


Wgj[<l>,Q,R]=0, (7.28a) 


Wgj[(l),Q,R]=igi[<l>,Q,R]+Shj[<l>,Q,R], (7.28b) 

where the operator Ig is defined in (7.20). In (7.28b) it 
acts on functionals of the variable (j). With (7.27) it can 
be written as 

(7.28c) 


The sum of the A5'-terms in (7.28c) give the part of 
A5'’®[G^ -I- </<, R] with at least one (/(-derivative acting 
to the right. The operator Ig defined in (7.28c) generates 
BRST transformations while keeping the regulator term 
fixed. Consequently the mSTI (7.28) entails that such a 
BRST transformation of I is given by the explicit BRST 
variation due to 51. The correlation function 51 is the 
expectation value of 51 defined in (7.23). Similarly to 
(7.17) we write 

Tlgj[J, R] = (s[(^] /[J, ^]) . (7.29) 

Eq. (7.29) entails that 51 vanishes for BRST-invariant 
correlation function I. In this case Igl = 0. Finally 
we remark that the representation (7.18) of VV translates 
into 


W../ = 


(4[^] + (4[^]r)) / + W, (7.30) 

where we have used 7(3^1 j] IE[J, i?] = r[(/), i?], 

following from (7.19), (7.25) and (7.27). 

We proceed with elucidating the general identity (7.28) 
with two examples. Firstly we discuss the standard regu- 
larisation with a quadratic regulator R’^^4’a4’h- Inserting 
this into (7.28) we are led to 

51 


mj = 


iL_9Rb«r — 
54>a 54>e 


5Q\ 


+ 51. (7.31) 


The second important example is provided by the mSTI 
for r. It can be read off from the alternative representa¬ 
tion for Wgj in (7.30) for / = 1 (/ = 1) leading to 


5r 

5(j)eL 


-AS^^'^[G^ + ^,R] Geb 


5r 

5 ( 1 )'^ J 


= 0. (7.32) 


We emphasise that the Q-derivative in (7.32) is that at 
fixed (/) and not at fixed J. It is also possible to derive 
it directly from (7.28) with I = (j). For the standard 
regulator the mSTI (7.32) reads [103] 


5T JF 


dF’'^ 

The terms proportional to derivatives of AS'' cancel in 
(7.33). 
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3. Flows and alternative representations 

The compatibility of (7.28) with the flow is ensured by 
the flow (3.60) for VV/, 


[dt + AS 2 ) Wi=0, (7.34) 


for the effective action and quadratic regulator see [104, 
110, 112, 113] 

Eq. (7.34) implies that a truncated solution to /sti = 0 
stays a solution during the flow if the flow is consistent 
with the truncation. Then it suffices to solve the mSTI 
for the initial condition /[^, Q, i?in]. How¬ 

ever, the search for consistent truncations is intricate as 
(7.28) involves loop terms. It is worth searching for alter¬ 
native representations of the mSTI (7.28) that facilitate 
the construction of such truncations. For the sake of sim¬ 
plicity we discuss this for the mSTI (7.33) for the effective 
action in the presence of quadratic regulator terms. The 
generalisation to correlation functions / and general AS' 
is straightforwardly done by substituting the correlation 
function T with / (leaving the T-dependence of L, un¬ 
changed) as well as the quadratic regulator R^'° with a 
general R. We can cast (7.33) into an algebraic form 
using the fact that —R serves as a current for G: 


W .,1 = 


<5r jr 

5Q^ SR'^^ ■ 


(7.35) 


The algebraic form of the STI (7.35) can be used to en¬ 
sure gauge invariance in a given non-trivial approxima¬ 
tion to r by successively adding explicitly i?-dependent 
terms. Such a procedure accounts for gauge invariance of 
classes of resummed diagrams. We add that in most cases 
it implicitly dwells on an ordering in the gauge coupling. 
We also remark that (7.35) seems to encode a preserved 
symmetry. This point of view becomes even more sug¬ 
gestive if introducing anti-fields [114, 115]. Note that 
in general the related symmetry transformation is inher¬ 
ently non-local. 

Eq. (7.35) constitutes an ordering in R. This can be 
made explicit by fully relying on the interpretation of 
i? as a current. There is a simple relation between Q- 
derivatives and J-derivatives: BRST variations of the 
fundamental fields (p are at most quadratic in the fields, 
see (7.8). Hence, the t^-order of the BRST transforma¬ 
tion of a composite held S()) is at most increased by one. 
Therefore, the source term sej)^ can be absorbed into 


a redefinition of 

+ Q" = (7.36) 

(j" + Q" (s<^c)^"^J k - . 

The tensors (s^c)’'^'^ are the structure constants of the 
gauge group as can be seen within the example of the 
fundamental fields (7.3) and their BRST variation (7.8). 
With (7.36) we can rewrite Q-derivatives of W and T 
in terms of J,i?-derivatives of W and i?-derivatives and 
fields (fi for T. The key relation is 

= + (7,37) 

where we also have to admit source terms with source 
—R for AiCa and CaCp. With (7.37) we can substitute 
the Q-derivatives in (7.35) and eliminate Q. Then the 
correlation function VVsp[^, i?] = >Vsp[0,0,i?] reads 

-2R>'> A; , (7,38) 

At i? = 0 the second line vanishes and we deal with 
the standard STI. The parameterisation (7.28) and (7.35) 
of the STI emphasise the gauge symmetry and are cer¬ 
tainly convenient within a coupling expansion. The pa¬ 
rameterisation (7.35) and (7.38) naturally relate to the 
’importance-sampling’ relevant in the flow equation. The 
latter, (7.38), requires no BRST source terms and hence 
reduces the number of auxiliary fields/terms. 

The derivation of (7.38) highlights the fact that (7.31) 
also constitutes the Slavnov-Taylor identity for the 2PI 
effective action, e.g. [175, 176]. To that end we restrict 
ourselves to a = a and quadratic regulators With 
the substitution we are led to the Slavnov- 

Taylor identity for More explicitly we 

have 

jab ^ _j^ab ^ (-7 gg) 

and 

+ r’’^a<l>b , (7.40) 

with the implicit definition = {4>a,4>bc = ‘t’bfkc)- We 
perform a second Legendre transformation with 

r 2 Pl[^i’a, (t>ab, Q] 

= sup {j-^^ab + r[().a, Q, 7?“*' = - , (7.41) 
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leading to and cj)ab = G. Note that 

r[(j)a,Q, R‘^^] already includes the standard subtraction 
—J°‘^(j>a(t>b- We arrive at 


<^r2Pi (5r2Pi i5r2Pi (5r2Pi’‘^ 

Ha 5Q°- 5(j)ab SQ^ 


(7.42) 


The last term on the Ihs of (7.42) accounts for the BRST 
variation of (pab that derives from the BRST variations 
of its field content papb- The BRST variation of papb 
can be added with a source term Q’^^sHapb) the path 
integral leading to r2Pi = T 2 Pi[ 4 )aHab,Qa,Qab]- Then 
we have 


We close this section with a short summary of the 
derivation of STIs without the use of BRST transfor¬ 
mations. To that end we integrate out the auxiliary held 
A and use the classical gauge-hxed action (7.7). In view 
of the auxiliary nature of the ghost helds we derive iden¬ 
tities that describe the response of general correlation 
functions to (inhnitesimal) gauge transformations of 
the physical helds, the gauge held and possible matter 
helds. Gauge-invariant correlation functions /, I are in¬ 
variant under these transformations. 

= i{Du)i, [u, C]„ , [a;, C]^). (7.47) 


i^r2Pi _ (5r2Pi’° 

“ SQ^ 



(7.43) 


Eq. (7.43) and the symmetry property pab = Hbpca lead 
to (7.42). Collecting the helds into a super-held pa = 
iPaHbc), and with r2Pi = T 2 Pi[pa,Q‘^], 

we get an appealing form of the STI (7.42) 


(5r2Pi (5r2Pi 

Spa SQ’^ 


(7.44) 


In its spirit (7.44) is close to the mSTI written as a mas¬ 
ter equation [114, 115]. As in (7.44) the master equation 
emphasises the algebraic structure of the mSTI but hides 
the symmetry-breaking nature of the identities. Nonethe¬ 
less algebraic identities are useful if constructing consis¬ 
tent truncations as well as discussing minimal symmetry 
breaking due to quantisation in the sense of Ginsparg- 
Wilson relations [29]. 

As in (7.38) we can absorb (^“-derivatives with help 
of (7.36), (7.37). As the source Q is a spectator of the 
Legendre transformation (7.41) we have and 

(7.37) reads for the 2PI ehective action 

- {,p-+ li,p-y^p,,) , (7.45) 


where we have used that = — J®** and hence 
—pab- Using (7.45) we arrive at 


(sr + (sPH^Hbc) 


jr2Pi (5r2Pi 

Spab <5(5“'’ 


0. (7.46) 


The BRST variation of pab involves pcdpe and (5“^ is a 
source for a specihc tensor structure pcdpe- Within 

regularisation of the 2PI effective action that regularises 
three point functions the source can be eliminated 
analogously to (7.38). This is an interesting option for 
A^PI regularisations of gauge theories, in particular in 
view of consistent approximations [174-177]. 


The linear operator g is bosonic as distinguished to s. It 
can be cast into the form (7.9) as a functional derivative 
operator g = {Q(pa)s^- The related generator Ig reads 

4 = (^"(fl0)a-(fl^M))^^^ , (7.48) 

leading to the STI (7.16) for Wsj. Restricting ourselves 
to J-independent Ps (7.16c) reads 

ri={Ql[P])[P=^]. (7.49) 

In the presence of a regulator term the generator of sym¬ 
metry transformations turns into 

4 = (j"(g0)a - (g(^ + AS))[P])^^^ , (7.50) 

leading to the mSTIs (7.22) and (7.28) for Wgj and Ws,/ 
respectively. We close this section with exemplifying the 
mSTI lUg.7 at / = 1 and the standard flow. Then, with 
the alternative representation (7.30) we are led to [113] 

gT[p,R] = + As)[G^ + P]^ 

-gAS'{P,R\ - [g{Cc.^D^,G 0 ){GA + p]J . (7.5I) 

The right hand side of (7.51) reproduces the gauge vari¬ 
ation of the classical action gA]^] as well as loop terms. 
The highest loop order (in the full propagator) is given 
by the highest order of the held p in the gauge fixing 
term and the ghost term in the classical action as well 
as the regulator term. For linear gauges and p = p the 
modified STI (7.51) involves one loop (gauge fixing, ghost 
term) and two loop terms (ghost term) apart from the 
regulator-dependent terms. Thus a purely algebraic form 
of the mSTI (7.51) can be achieved for regulator terms 
with R involving 7?“'’ and 
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C. Gauge-invariant flows 

An interesting option for flows in gauge theories is the 
construction of (partially) gauge-invariant flows. The 
gain of such formulations is twofold. Firstly they allow 
for a more direct computations of physical observables. 
Observables are gauge-invariant as opposed to Greens 
functions in gauge-fixed formulations. Secondly one can 
hope to avoid the subtleties of solving the symmetry re¬ 
lations in the presence of a regulator. However, gauge- 
invariant formulations come to a price that also has to 
be evaluated: if the corresponding flows are themselves 
far more complicated than the standard gauge-fixed flows 
the benefit of no additional symmetry relations is, at least 
partially, lost. Also, gauge invariance does not rule out 
the persistence of non-trivial symmetry relations, mostly 
formulated in the form of Nielsen identities or, alterna¬ 
tively, in the form of specific projections of the general 
Dyson-Schwinger equations valid within such a setting. 

In the present work we concentrate on gauge-invariant 
flows formulated in mean fields and the effective action 
Ffc. An alternative construction of gauge-invariant flows 
is based on the Wilsonian effective action SeS^, (3.40), 
formulated in Wilson lines and using gauge-covariant reg¬ 
ulators. For details we refer the reader to [135-139] and 
references therein. 

1. Background field flows 

The first and most-developed gauge-invariant flow 
originates in the use of the background field formalism. 
We couple the fundamental fields to the currents, (f> = ip 
with 

((. = (ai,C'„,C'„), (7.52) 

where the full gauge held is defined as 

A = A + a. (7.53) 

The gauge field A is split into a background held configu¬ 
ration A and a fluctuation field a coupled to the current. 
BRST transformations and gauge transformations are de¬ 
fined by (7.8) and (7.47) respectively at fixed background 
held A, sA = qA = 0. Note that the covariant deriva¬ 
tive reads D = D{a -\- A). Therefore, the mSTIs (7.28) 
for yVsj and Wg,/ persist. Within appropriate gauges, 
e.g. the background field gauge T = D{A)a, there is an 
additional symmetry: the action (7.7) is invariant un¬ 
der a combined gauge transformation of the background 


held A A + D{A)uj and the fluctuation field a [w, a]. 
This invariance follows by using that the fluctuation held 
(j) in (7.52) as well as the covariant derivatives D{A) and 
D[A transform as tensors under this combined transfor¬ 
mation. Defining background held transformations 

= {-DiA)iv ,0,0, D{A)u;) , (7.54) 

the transformation properties under the combined trans¬ 
formation are summarised in 

(0 + 0)<.^(7’, D{A), D{A)) = [w, ((p, D{A), D{A))\. 

(7.55) 

with g defined in (7.47). As the action S in (7.7) with 
T = D{A) or similar choices can be constructed from 
((/?, D{A ), D{A) ) this leads us to 

(0 + 0)5[0,A]=O, (7.56) 

Then, the corresponding effective action r[(/). A] is invari¬ 
ant under the above transformation, in particular we de¬ 
fine a gauge-invariant effective action with 

r[A] =r[(/) = 0,A]. (7.57) 

We have (g -I- g)r[(/), A] = 0, where 0,g act on <p = ip 
according to (7.47) and (7.54). This implies in particular 
gr[A] = 0. The gauge invariance of r[(/). A] persists in the 
presence of a regulator if AS'[a, i?(A)] is invariant under 
the combined transformation of a and A. This is achieved 
for regulators R that transform as tensors under gauge 
transformations A ^ A + D{A)uj. This amounts to the 
definition of a background held dependent R{A) with 

0i?(A) = [w , R{A)] . (7.58) 

For example, standard flows follow with the regularisa- 
tion AS = R^AA)aiaj + R°'^{A)CaCf 3 . The invariance 
property (g -I- g)A5' = 0 follows immediately from (7.55) 
and (7.58). The relation (7.58) is e.g. achieved for regula¬ 
tors in momentum space depending on covariant momen¬ 
tum D{A). Correlators I still satisfy the modified STI 
(7.28), but additionally there is a modified STI related to 
the background held gauge transformations (7.54). The 
related generator is 

4 = (•/"(00)a - (0(^ + AS)M) , (7.59) 

leading to mSTIs (7.22) and (7.28) for Wg,/. For the 
effective action (7 = 1) the mSTI reads 

emR] = + <(-]) 

/ _ f) 'P^ \ 
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Adding (7.51) and (7.60) we arrive at 

(0 + 0)r[0,A,i7] = O. (7.61) 


The derivation makes clear that, despite background 
gauge invariance (7.61), the effective action r[((), A] still 
carries the BRST symmetry (7.28) displayed in or 
Wg, where the background field is a spectator sA = 0. 
In other words, the non-trivial relations between Ap¬ 
point functions of the fluctuation field are still present. 
However, for fV-point functions in the background field 
they play no role which has been used for simplifications 
within loop computations. 

Therefore it is tempting to use these features for the 
construction of gauge-invariant flows. General flows 
within such a setting are still provided by (3.60). In 
particular with (3.63) we arrive at the flow of rfe[A] as 
[42, 43, 98-100, 106, 111-113, 132-134] 


f fe[A] = (A^[G^ + ct >, i?(A)])^^o - A^'[0, R{A )]. 

(7.62) 

It has already been discussed in [42, 43, 112] that the 
flow (7.62) is not closed as it depends on 

AM, ,,,3) 

the propagator of the fluctuation field, and possibly 
higher derivatives w.r.t. (j) evaluated at vanishing fluctua¬ 
tion field (p = 0. The Ihs of (7.62) cannot be used to com¬ 
pute this input as it only depends on A = A. Moreover, 
as has been stated above, these A^-point functions still 
satisfy the modified Slavnov-Taylor identities discussed 
in the last section. The differences between r(^^[A] and 
the fluctuation propagator (7.63) become important al¬ 
ready at two loop. The correct input (7.63) at one loop 
was used to compute the universal two loop /3-function 
which cannot be reproduced by using r^^)[A] [42, 43]. 
Still one can hope that qualitative features of the theory 
are maintained in such a truncation. Then, a measure 
for the quality of such a truncation is given by the dif¬ 
ference between a derivative w.r.t. A and one w.r.t. a of 
the effective action. This relation reads [42, 43, 112] 


and can be understood as a Nielsen identity. Eq. (7.64) 
relates Green functions of the background field with that 


of the fluctuation field. The latter satisfy mSTIs whereas 
the former transform as tensors under gauge transforma¬ 
tions reflecting gauge invariance. Hence, (7.64) encodes 
the mSTIs. Note also that the background field depen¬ 
dence stemming from the regulator should be understood 
as a parameter dependence and not as a field dependence 
An improvement of the current results in gauge theo¬ 
ries [42, 43, 98-100, 105, 106, 108, 109, 111-113, 132-134] 
requires an implementation of the Nielsen identity (7.64) 
beyond perturbation theory. 

It is possible to enhance background field flows to fully 
gauge-invariant flows with standard STIs by identifying 
the background field with a dynamical field. There are 
two natural choices: A = A and A = (A) = A. The 
latter leads to the definition of the effective action as a 
higher order Legendre transform. Then we have addi¬ 
tional terms to those (3.60) as 

r’“ = . (7.65) 

With (7.65) we get additional terms in the relations 
between ^-derivatives of T and J-derivatives of W. 
Eq. (7.65) is actually implementing the Nielsen identity 
(7.64) on the level of the Legendre transformation. This 
entails that in particular the basic relations (3.45) and 
(3.46) receive modifications originating in (7.65). As an 
example we study the standard flow for the effective ac¬ 
tion which reads 

f[A,(/.] = R‘^^W,ab 

= k^^Gab + -terms^ , (7.66) 

where the propagator G is defined with G = l/(r*^^^ -|- 
AS”*^^^). The propagator G of the dynamical field trans¬ 
forms as a tensor under gauge transformations reflecting 
gauge invariance. However, it can be shown in a pertur¬ 
bative loop expansion that effectively the flow equation 
can be rewritten as that in the background field formal¬ 
ism: the effective propagator W'A) {wA)Y behaves 
as that of the fluctuation field in the background field 
formulation. This is already indicated in (7.65). The 
correction terms involve the same correlation functions 
already relevant in the Nielsen identity (7.64). So still 


For infrared diverging regulator R(A) even the computation of 
the one loop /3-function requires a subtraction of the field depen¬ 
dence of R(A) [42, 43, 112], 

This choice can be only used in the regulator. 
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we deal with non-trivial symmetry identities. Nonethe¬ 
less the above formulation furthers the knowledge about 
truncation schemes that expand about F’® = J“, or al¬ 
ternatively about A, i?] = 0. Details shall 

be provided elsewhere. 

The other suggestion A = A relates to the use of a 
regulator term AS'[A, Such a regulator term can 

be written as AS'[A, R] = A5'[A, where _ 

AS"*! [0, i?(0)]/(n!) is the nth expansion coefficient in 

a Taylor expansion of AS'[A, i?(yl)] in the gauge field A. 
This flow is covered by the general flow (3.60) and in¬ 
volves all loop orders in the full propagator. Again this 
effectively reduces to the background field flow and comes 
at the expense of an infinite series of loop terms in the 
flow. In this context we remark that the latter set-back 
is avoided within the Polchinski equation. This follows 
in the present setting with (3.40) and the flow (3.28) for 
the Schwinger functional. 


2. Geometrical effective action 


We have seen in the last section that the flow of the 
gauge-invariant effective action within the background 
field formulation is not closed. In the process of curing 
this problem we encounter the persistence of non-trivial 
symmetry relations, conveniently summarised in (7.64). 
Both aspects originate in the fact that the sources are 
coupled to fields that do not transform trivially under 
gauge or BRST transformations. Hence the question 
arises whether one can do better. Within the framework 
of the geometrical or Vilkovisky-DeWitt effective action 
the fields cj) coupled to the sources are scalars under gauge 
transformations. 

Then, gauge-invariant flows can be formulated [140, 
141]. We do not want to go in the details of the general 
construction that can be found in [141]. The configura¬ 
tion space is provided with a connection Ty (Vilkovisky 
connection) which is constructed such that the disentan¬ 
glement between gauge fibre and base space is maximal. 
The gauge fields Ai are substituted by geodesic normal 
fields (pi that are tangent vectors at a base point (back¬ 
ground field) A. As a consequence the geodesic fields (pa 
tangential to the fibre drop out of the path integral, only 
the fields (pA tangential to the base space remain and 
are gauge-invariant. This construction is lifting up the 
relation between fluctuation field and background field 
(7.53). The linear background relation can be read as 
the limit in which the connection Ty is neglected. The 


full relation reads schematically 

(Pi = A,-Ai + Fv (Pj(Pk + Oi(Pf ), (7.67) 

with Q(p^ = 0 = Qtp^- This is used to construct a gauge- 
invariant effective action r[(p^,A,R] which is gauge- 
invariant under both sets of gauge transformations g and 
g [141]. Again a gauge-invariant effective action in one 
field can be defined as F[A, i?] = F[^ = 0,A, i?]. The 
flows of T[(p,A,R] and F[A, i?] are given by (3.63) and 
(7.62) respectively, both being gauge-invariant flows. We 
still have a Nielsen identity equivalently to (7.64). In the 
underlying theory without regulator term it reads 

r* + F,,(r.i) = 0, (7.68) 


where (p'^-i stands for the covariant derivative with the 
Vilkovisky connection Fy. The related symmetry opera¬ 
tor is provided by 


In 



(7.69) 


With (3.8b) this turns into 
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+ (j“ - A^'“*'[G^]G6e^) ^.;4GA + <(,], (7.70) 

in the presence of the regulator term. For standard flows 
the choice Wn,i in (7.30) reproduces the Nielsen identity 
derived in [141], 


Ffc,, = iG'*" Rba,i + (Ffe,a - Rab • (7.71) 


For more details and its use within truncation schemes 
we refer to [141]. The formalism discussed above provides 
gauge-invariant flows that are closely linked to the back¬ 
ground field formalism (in the Landau-DeWitt gauge) as 
well as to standard Landau gauge. This comes with the 
benefit that results obtained in the latter can be partially 
used within the present formalism. Indeed the present 
setting can be used to improve the gauge consistency of 
these results. We hope to report on results for infrared 
QCD as well as gravity in near future. 

To conclude, we have discussed the various possibility 
of defining gauge-invariant flows and their relations to 
gauge-fixed formulations. These relations come with the 
benefit that it allows to start an analysis in the gauge- 
invariant formulations on the basis of non-trivial results 
already achieved in gauge-fixed settings, one does not 
have to start from scratch. 


43 



D. Chiral symmetry and anomalies 

We want to close this chapter with a brief discussion of 
FRG flows in theories with symmetries that are flawed by 
anomalies on the quantum level, e.g. [106-108, 114], A 
more detailed account shall be given elsewhere. In partic¬ 
ular a discussion of the chiral symmetry breaking requires 
a careful investigation of chiral anomalies. The deforma¬ 
tion of the chiral symmetry from a general RG transfor¬ 
mation has already been considered in [29], and leads to 
the Ginsparg-Wilson relation This has been empha¬ 
sised in [114]. A discussion of chiral symmetry breaking 
requires a careful investigation of chiral anomalies. Inte¬ 
grated anomalies are tightly linked to topological degrees 
of freedom like instantons via the index theorem. FRG 
methods have been shown to be sensitive to topologi¬ 
cal degrees of freedom [105, 106], an interesting quantum 
mechanical example can be found in [148]. In the present 
section we consider the gauge field action (7.7) together 
with a Dirac action 

SD[<l>]=MP + mr’^i;b, (7.72) 

with a possible mass term and (j) = {A, C, C, ip, ip). The 
Dirac operator Ip reads 

lp'^'> = {^ + 7^4)“*' (7.73) 

with the free Dirac operator ^ and a coupling to the gauge 
field with a possible projection V either proportional to 
the identity "P = F, or projecting on right- or left-handed 
Weyl fermions V± = Here we consider 

P = P+ = i^, m = 0, (7.74) 

a theory with left-handed Weyl fermions coupled to a 
gauge field, and free right-handed Weyl fermions. The 
symmetry transformation that leaves the action (7.72) 
invariant is given by 

S+P = (0+^, Q+C , 0+G, uV-ip, -PpV+uj). (7.75) 

The transformations (7.75) cover both, BRST transfor¬ 
mations with = s with u) = C, and g+ = g with gauge 
transformation parameter w. Here we stick to g+ = g. 
The chiral anomaly comes into play since the fermionic 
path integral measure dip dip is not left invariant under 


The derivation in [29] makes no use of the lattice. 


the transformation (7.75). In other words, (7.75) is not 
unitary. We quote the result 

s{dip dip) = dip dip , (7.76) 

with infinitesimal variation a;“Aa. The non-Abelian 
anomaly A reads 

A“(a:) = {dpA„Fp„ - \A,,ApA„) .(7.77) 

Then, the generator of gauge transformations /g in (7.48) 
receives a further contribution and reads 

h = (•7"(0^)a - (0^^) - A[P]) . (7.78) 

and with (3.8b) in the presence of the regulator term we 
arrive at 

/g = (j‘^(g^)a-(g(^[(^] + A^))-AM). , .(7.79) 

Eq. (7.79) can also be read off from (7.50) since the 
anomaly term A[(p] commutes with AS”. 

We conclude with briefly discussing the {7/i(l)-anomaly 
relevant for anomalous chiral symmetry breaking. We 
restrict ourselves to standard flows with quadratic regu¬ 
lator. The Dirac action (7.72) with P = F is invariant 
under global axial {/^(Ij-transformations. The related 
Noether current is derived from the C/t( 1) transforma¬ 
tions of the fermions 

qaP = W75V’, 9aP = V’75W . (7.80) 

The rest of the fields transforms trivially with qaA = 
qaC = qaC = 0. The related anomaly reads 

A = ^Cp,ypcrtr Fpi/Fp,y . (7.81) 

The anomalous Ward identity for the effective action, 
VVg^4, in the presence of the regulator reads 

(0.4</>)ar’“ + {9a{Sd + A5))[G^ + <P] 

= (A[Gi^ + <P]) ■ (7.82) 

The space time integral of (7.82) produces the (expecta¬ 
tion value of the) topological charge on the rhs, as well 
as the analytical index of the modified Dirac operator on 
the Ihs. In [106] it has been shown that the number of 
zero modes stays the same for regulators with chiral sym¬ 
metry. The chiral anomaly has been investigated in [107]. 
In general the Ihs of (7.82) is computed directly from the 
effective action. Accordingly we can use (7.82) for testing 
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the potential of given truncations to the effective action 
for incorporating the important topological effects. Ad¬ 
ditionally its provides non-trivial relations between the 
couplings. For example, the leading order effective action 
derived in [106] satisfies (7.82) up to sub-leading terms 
(in 1/k). Eq. (7.82) can be used to determine coefficients 
and form of these sub-leading terms, in particular in view 
of CP-violating effects. 


VIII. TRUNCATION SCHEMES AND 
OPTIMISATION 

The reliability of results obtained within the functional 
RG rely on the appropriate choice of a truncation scheme 
for the physics under investigation, as well as an optimi¬ 
sation of the truncation with the methods introduced in 
section V. The truncation has to take into account all 
relevant operators or vertices. In theories with a com¬ 
plicated phase structure this might necessitate introduc¬ 
ing a large number of vertices to the effective action in 
terms of the fundamental fields. A way to avoid such a 
drawback is to reparameterise the theory in terms of the 
relevant degrees of freedom [41, 72-76, 78-82]. 

Fixed point quantities like critical exponents and gen¬ 
eral anomalous dimensions have very successfully been 
derived within the flow equation approach, mostly in the 
derivative expansion, see reviews [15-17, 19-22]. For the 
evaluation of these results in view of quantitative reliabil¬ 
ity one has to assess the problem of optimisation. To that 
end we evaluate the consequences of the relation between 
RG scaling and flow for an appropriate choice of classes of 
regulators. As an example for the optimisation criterion 
developed in section V, we discuss functional optimisa¬ 
tion within the zeroth order derivative expansion. The 
unique optimised regulator is derived and its extension to 
higher order of the truncation scheme is discussed. For 
explicit results we refer the reader to the literature, in 
particular [68]. 

A. Field reparameterisations 

The derivation of the flow in section III was based on 
a bootstrap approach in which the existence of a renor¬ 
malised Schwinger functional in terms of the possibly 
composite fields </> was assumed. This already took into 
account that the fundamental fields ip may not be suit¬ 
able degrees of freedom for all regimes of the theory under 


investigation. For example, we could consider fields 4'{'p) 
that tend towards the fundamental fields in the pertur¬ 
bative regime for large momenta, 

( 8 . 1 ) 

while being a non-trivial function of ip in the infrared. 
This includes the bosonisation of fermionic degrees of 
freedom [72, 73, 78, 80], e.g. in low-energy QGD, where 
the relevant degrees of freedom are mesons and baryons 
instead of quarks. More generally such a situation applies 
to all condensation effects. 

In such a case the Green functions of ip will show a 
highly non-trivial momentum dependence or even run 
into singularities. Moreover, physically sensible trunca¬ 
tions to the effective action in terms of ip could be rather 
complicated. These problems can be at least softened 
with an appropriate choice of (j) that mimics the relevant 
degrees of freedom in all regimes. Such a choice may 
be adjusted to the flow by implementing the transition 
from p to 4>{p) in a fc-dependent way [72, 73]. This can 
be either done by coupling the current and the regulator 
to a /c-dependent field (j)k, or by choosing a fc-dependent 
argument (j)k of the effective action F^: 

The former option leads to additional loop-terms in the 
flow. The relation (3.18) is modified as the full Schwinger 
functional W[J] couples to a fc-dependent field (f)k with 
dtW[J] = J‘^{dt4>ka), and the flow operator AS'2 changes 
as the regulator term has an additional /c-dependence via 
the field, AS'2 [(j), R] —> AS'2 [^, R'] where i?' is defined 
with 

AS[$k, R'] = AS[$k, R] + dt$ke. A5’"[(^fe, R ], (8.2) 

where dt(i>ka. = With these modifications the 

derivation of the flow can straightforwardly be redone. 

The latter option keeps the flow (3.60) as the partial 
derivative is taken at fixed argument (j): dtik = dt\^Ik- 
For integrating the flow the total derivative is required, 

= -AS2[(^k,R]ik[M+dt^kJ'k[<^k] ■ (8.3) 

We can also combine the above options. For the sake of 
simplicity we restrict ourselves to the flow of the effective 
action which reads in this general case 

= (A5[G^ + , R'\) - R'\ 

+ [dt(k.-{dtk))K- ( 8 - 4 ) 
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In (8.4) we dropped the subscript k with (j) = <j)k- The 
first term on the rhs is the expectation value of AS”);/), i?'] 
defined in (8.2). The second term originates in the defini¬ 
tion of Tfc in (3.43). The expectation value in the second 
line in (8.4) can we written as {dtfj)) = -f (j)]), 

and R' is defined in (8.2). We remark that (8.4) is finite 
for fc-dependences of (j) that are local in momentum space. 
General /c-dependences may require additional renormal¬ 
isation. The flow (8.4) can be used in several ways to 
improve truncations. 

A given truncation scheme can be further simplify in 
a controlled way by expanding the effective action about 
a stable solution (j) of the truncated equations of motion, 
= 0. Then the second line in (8.4) is sub-leading 
for (j) — (j) small and can be dropped if restricting the flow 
to the vicinity of (j). As this is an expansion about a 
minimum of the effective action, such a truncation has 
particular stability. 

The second line also vanishes for dt(t>— {dtcj)) = 0. Sub¬ 
ject to a given (j) we demand (j) to satisfy 

{dt$)=dt(l>. (8.5) 

With (8.5) the second line in (8.4) vanishes identically 
and the flow reduces to the first line. The construction of 
R' requires the knowledge dt4>{4>). Within given trunca¬ 
tions (8.5) turns into a set of loop constraints that accom¬ 
pany the flow. These constraints resolve the dependences 
of the flowing composite fields ipk on the microscopic de¬ 
grees of freedom. This is more information than required 
for solving the flow. Indeed, we also can use (8.5) to cir¬ 
cumvent the necessity of finding We write for the 

expectation value of the second term in (8.2) 

AS^^[^+cl),Rh\{dtk) 

= {AS’^[G^+<P,Rh\dt<l>), ( 8 . 6 ) 

where we have used (3.50) and (8.5). With (8.5) and 
(8.6) we can substitute all dependences on in the 

flow (8.4) by that on (l)k,dt(l)k- We are led to a closed 
flow for the effective action 

dtTk[<P] = {AS[GA+^, i?]) _ AS'[<I>,R'] 

+ (AA’‘’[GA + R]j\ dt(j>) - . (8.7) 

The first term in the second line keeps track of the k- 
dependence in (f>k necessary to satisfy (8.5). The last 
term carries the /c-dependence of (pk- For the standard 
quadratic regulator (8.7) reads 

dtTklcp] = + 2i?-‘’Gac<^b’" - k^k . (8.8) 


We illustrate the above considerations within simple ex¬ 
amples for quadratic regulator terms (8.8). Furthermore 
the examples are based on linear relations between dtp 
and p. Then (8.5) can be resolved explicitly and up to 
rescalings (8.8) simplifies to the standard case: we ab¬ 
sorb a t-dependent wave function renormalisation 

1 /2 

into the field: pk = Z^ pQ with dtpk = 'ftjypk with 
'y<p = dtlnZfj,. Eq. (8.5) is satisfied with pk = Z^ pQ. 
Then (8.8) reduces to 

(j)t + '^cj>Pa.-^f^ Ffc)^] = Gbc(3t -l- , (8.9) 

which also can be obtained by explicitly using pk = 
Z^ p. The flow (8.9) also makes explicit that the trans- 
formation p Z^ p is a RG rescaling. This procedure 
can be used to fix the flow of vertices. 

Another simple example is the expansion of the effec¬ 
tive action about its minimum at pminik), implying 
p ^ pk = P — Pminik). Such a reparameterisation guar¬ 
antees that the minimum is always achieved for pk = 0. 
The flow (3.60) only constitutes a partial t-derivative, as 
it is defined at fixed fields p. With pk = p — pmin{k) we 
satisfy (8.5) and we are led to (8.7) with dtpk = —dtpmin 
with dtph’'^ = 0. The flow (8.7) reduces to the standard 
flow, 

dtTkiP] = {AS[G^ +P,R])- ^S'iP, R] 

+Ff [</*] (at0mi„)a , (8.10) 

now describing a total t-derivative of the effective action 
Tfe. For quadratic regulators R’^'^ it reads 

dtTklP] = R'^’Gab + G^IP] {dtPmin). ■ (8.11) 

The flow of the minimum pmin can be resolved with help 
of ['^min]) = 0, and reads 

{dtpinin)a = ( -75^- ) ■ (8.12) 

V Tfc ['(’min] / ab 

As mentioned before, the examples used linear depen¬ 
dences of dtp on p. Then (8.8) can also be derived ex¬ 
plicitly as p is known. In the general case this is not 
possible, and (8.7) or (8.8) are the fundamental flows. 

B. RG scaling and optimisation 

The reliability of results obtained within functional 
RG flows hinges on an appropriately chosen truncation 
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scheme and a regulator choice that optimises the given 
truncation scheme. Without specifying the truncation 
scheme the following observation can be made: the renor¬ 
malisation group analysis in section IV relates the RG 
equation for the full theory with that in the presence of 
a regulator. In particular we deduce from (4.25) and by 
identifying s with the RG scale /i, that the RG equation 
for the regularised effective action reads 

■ (8-13) 

The right hand side of (8.13) entails the modification of 
the RG properties in the presence of the regulator. In 
(8.13) we have restricted ourselves to quadratic regula¬ 
tors. As explained in detail in the context of optimisation 
in chapter V, for full flows without truncations differ¬ 
ent choices of regulators, in particular those with differ¬ 
ent RG properties, lead to a RG rescaling of fields and 
coupling in the full effective action T. However, within 
truncations this modification usually leads to a physical 
change of the end-point of the flow. In turn, this problem 
is softened if restricting the class of regulators to those 
with [42, 43] 


(D^+j^)E = 0, (8.14) 

where (y^i?)®*^ = 270®cR°’^- The constraint (8.14) leads 
to 


to normalise It could be proportional to eval¬ 
uated at some momentum, e.g. at vanishing momentum. 
By construction (8.16) satisfies (8.14) as the two-point 
function does, (Dg + = 0. If evaluating the 

standard flow (3.74) for the effective action at the back¬ 
ground field (j), it takes the simple form 



r(2) 


be 


(8.17) 


where for the sake of simplicity we have taken = 

T’'^’^]^], that is no subtraction. The first term on the 
rhs of (8.17) can be integrated explicitly and contributes 
to the effective action only at perturbative one loop or¬ 
der. The second term gives non-trivial contributions if 
the spectral density changes. Eq. (8.17) is a spectrally 
adjusted flow. 

In most truncation schemes used in the literature 
(8.16) simply amounts to the multiplication of the wave 
function renormalisation Z^. Then the propagator fac¬ 
torises G[Z(j,] = Z^^G[1] which facilitates the computa¬ 
tions. It is for the latter reason that (8.14) is a standard 
choice for regulators and it is a fortunate fact that the 
simple structure of flows for the choice (8.14) goes hand 
in hand with better convergence towards physics. 


D^rk = 0. (8.15) 

For the class of regulators with (8.14) the regularised 
correlation functions satisfy the same RG equation as in 
the underlying full theory, in particular this holds for the 
effective action, (8.15). Apart from the general optimi¬ 
sation arguments made above this facilitates the identi¬ 
fication of anomalous dimensions and critical exponents. 
Indeed, the choice (8.14) with the additional identifica¬ 
tion t = In /i allows for the straightforward identification 
of t-running and RG running within fixed point solutions 
at all orders of the truncation scheme. 

An explicit example for a class of regulators in stan¬ 
dard flows that satisfy (8.14) is provided by [42, 43] 

(8.16a) 

with 

D^r = 0, (8.16b) 

where is evaluated at some background field 0, 
with a possible subtraction. The subtraction can be used 


C. Integrated flows and flxed points 

An optimisation with (5.29b) requires the minimisa¬ 
tion of the norm of the difference between the regularised 
propagator and the full propagator with the constraint 
of keeping a fixed gap, see (5.32a). This implies a fine- 
tuning of the regulator in dependence of the two-point 
function, T^ab- Here we outline a way of solving the flow 
equation which naturally incorporates such a task and 
hence minimises the additional numerical effort. First 
we turn the flow (3.60) into an integral equation 

Iq = Ia 7- f dt AS 2 lk , (8.18a) 

Ja 

where A is the initial cut-off scale and the integrated flow 
for the effective action derives from (3.63) as 

ro = rA +£dt(^{AS[G A+ ^,R]) 

-kAS"[^,i?]) . (8.18b) 
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Eq. (8.18) constitutes DSEs as already explained in sec¬ 
tion VI A. As distinguished to standard DSEs they only 
involve full vertices and propagators. Such a set of equa¬ 
tions can be solved within an iteration about an ansatz 
for the full flow trajectory R{k)]. The better such 

an ansatz fits the result, the less iterations are needed 
for convergence towards the full result R{k)]. A 

benefit of such an approach is that it facilitates an im¬ 
plementation of the optimisation criterion (5.26) in its 
form (5.32a). After each iteration step we can prepare 
our regulator according to (5.32a) for the next step. Such 
a preparation is in particular interesting for truncations 
with a non-trivial momentum dependence for propaga¬ 
tors and vertices. Furthermore the integral equations 
(8.18) are likely to be more stable in the vicinity of poles 
of the propagator. 

The integral form (8.18) also is of use for an analy¬ 
sis of asymptotic regimes and in particular fixed point 
solutions. In general functional RG methods have been 
very successfully used within computations of physics at 
a phase transition. In particular critical exponents can 
be accessed easily. 

At fc = 0 the flows (3.60) have a trivial fixed point, 
i9t/|fc=o = 0. In case the theory admits a mass-gap Agap, 
this can be used to resolve the theory below this scale 
hence getting access to the deep infrared behaviour. For 
the sake of simplicity we further assume dimensionless 
couplings. The dimensionful case will be discussed else¬ 
where. Then, in the regime 

e « A^^p , (8.19) 

the flow of correlation functions Ik is parametrically sup¬ 
pressed by powers of fc/Agap, 

dth = 0(A:/Agap). (8.20) 

Eq. (8.20) applies in particular to the effective action 
and its derivatives. It is convenient to parameterise the 
correlation functions R as 

Ik = ioil + 6ik). (8.21) 

Inserting this parameterisation into the integrated flow 
(8.18) we arrive at an integral equation for SR, 

SR = -J° dt' AS 2 (/o(l + Sik')) , (8.22) 

where AS'2 depends on T^®*^ (and its derivatives) that 
admit the same parameterisation (8.21). Assume for the 


moment that SR' and on the rhs of (8.22) only 

depend on dimensionless ratios 



where the pi are momenta of the correlation functions R, 
e.g. external momenta of n-point vertices. This assump¬ 
tion reads 

SR = SI[pi, ...,p„] + 0(k/Agap). (8.24) 

Inserting (8.24) into the rhs of the integrated flow (8.22) 
we deduce from a scale analysis that the resulting SR 
on the Ihs can only depend on dimensionless ratios R. 
A good starting point for the iteration is R = R with 
SI = 0. Such a choice trivially only depends on the ratios 
(8.23). Hence this holds true for each iteration step, and 
we have proven (8.24). 

Now we invoke the optimisation (5.26) with Dr^I = 
IoDr±SI, and we are led to the constraint 

f dt'AS 2 (loiDR'jik')) =0. (8.25) 

J ' -Rstab 

For positive definite SI a solution to (8.25) is given by 
AS 2 SI = 0. In this context we remark that SI is not a 
correlation function /, and the above resolution does not 
imply a vanishing flow of SI. An optimisation along these 
lines was put forward in the infrared regime of QCD, for 
details see [128, 129]. 

D. Optimisation in LPA 

We continue with a detailed analysis of the optimi¬ 
sation (5.26), (5.29) in the LPA a scalar theory with a 
single scalar field a = x. We shall show that within the 
LPA the regulator (5.11) follows as the unique solution to 
(5.28), see also the more explicit form without RG scal¬ 
ing, (5.17). For the sake of simplicity we use the standard 
flow (3.72) with AS 2 R = {G RG)bcik^ ■ In the LPA we 
have to evaluate (5.28) for constant fields. Moreover we 
consider correlation functions R that are functionals of 
(j) and not operators. For example, in the present trun¬ 
cation scheme relevant correlation functions are provided 

by 

j = (J , (8.26) 

and combinations thereof. In LPA all quantities are eval¬ 
uated for constant fields R On the rhs of the standard 
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flow (3.72) the second derivatives are required. In 
LPA they are parameterised as 

=^k{^,p'^)Sip- q). (8.27) 

We also need the full propagator G{p, q) = — 

which reads l/{T^k\$\ + R){p,q) = 
l/(p^ + V"[4i] + R{p^)) 5{p — q). Inserting these objects 
into (5.17) we arrive at 


rpaiaa S{G RG)bc 


( d^W] 

\ (5(/>)2 


cb 


R —^stab 


0, (8.28) 


which we recast in a more explicit form 


/ 


d'^p 


<5i?±(/) 


6 

SRip"^) 


Ik 


(8.29) 


X 


d‘^q _ Ikitj), q^) 

{2ttY + R{qY + V”[$\? 


dtR{q^) 


= 0 . 


^—-Rstab 


Now we use that a general regulator R can be writ- 
ten as R{qY = q^r{x) with x = q'^/k'^, if no fur¬ 
ther scale is present in R. This entails that dtR = 
q^dtr{x) = q^{—2x)dxr(x). Furthermore we can rewrite 
the integration over q as one over x: d'^q/{2'K"^Y = 
d^lddxx ^‘^/'^~'^'>With these identifications we get for 
the g-integral in (8.29) after partial integration 


^dlk{(t), 0)S2d + 



dxx‘^^‘^ ^Ik^{{d/2 — \) 


~\~xdx InT/j; 


r + V"lx 
1 -I- r -I- V"/x 


V" IX 

{l + r + V"/xf 


} . (8.30) 


Now we are in a position to discuss the extrema (8.28). 
Searching for minimal flows is equivalent to searching for 
r that minimise the absolute value of the integrand in 
(8.30) 


min 

r 


{{d/2 - 1) + xdx\nXk) 


r + V'/x 
l + r + V'/x 


V'/x 

{l + r + V"/xf ’ 


(8.31) 


where we have left out the overall factor x'^l‘^~‘^Tk- A 
simpler condition is achieved by neglecting the model- 
dependent second term proportional to V'/x leading to 


min 

r 


r + V'/x 
1 + r -h V"/x ■ 


(8.32) 


We proceed with the extremisation of the full integrand 
by taking the r-derivative at fixed Ik of the function in 


(8.31). We arrive at 


{{d/2 — 1) -I- xdx lnXfc)(l + r + V"/x) + V'/x 
(1 + r -I- V"/xY 


(8.33) 


We remark that subject to (((i/2 — 1)-|-a;i9a; InZ^) > 0 
and V"/x > 0 the r-derivative (8.32) is positive. Note 
also that r + V"/x > 0 cannot be obtained for all x 
and 0 if the potential V is not convex yet. This state¬ 
ment holds for all regulators'^. However, for optimised r 
the region V"/x < 0 for x should have small impact on 
(8.30). If (i > 4 we regain positivity for vanishing or pos¬ 
itive dx InZfc. Leaving aside this subtlety we solve (8.31) 
for positive regulators r. As its derivative is positive, 
(8.33), this amounts to minimising r 


rsta.b<r, Vr,x. (8.34) 


So far we have not used the definition of in (5.26). 
With its use we are straightforwardly led to (8.36). Still 
we would like to evaluate how unique or natural the 
choice R± is. If r was an arbitrary positive function of x, 
(8.34) leads to r{x) = 0. However, as r has been intro¬ 
duced as an IR-regularisation it is inevitably constrained: 
it entails an IR-regularisation in momentum space only 
with 


X -I- xr{x) > c (8.35) 

for some positive constant c. For a proper IR-regulari- 
sation the full propagator G has to display a maximum 
G < l/cmin with Cmin = c + V//^^ > 0, where is the 
minimal value of V", possibly negative. For momenta 
X > c the solution of (8.34) with (8.35) is r{x > c) = 0. 
For X < c we saturate the inequality (8.35) with r{x) = 
c/x — 1. This leads to a unique solution rstab of (8.34) 
for r € {rx} defined by (8.35). 

rstab{x) = {c-x)d{c - x), (8.36) 

which is equivalent to (5.11). Note that in between (8.35) 
and (8.36) we have implicitly introduced the set {R±} of 
(5.26) by keeping c fixed while minimising r. Still, such a 
procedure was naturally suggested by the computation. 

Above we have restricted ourselves to correlation func¬ 
tions Ik with {{d/2 — 1) -I- xdx InZfe) > 0. If we discuss 
optimisation on the set of f d'^x , (8.26) they lead 


All regulator functions have to decay with more than 1/x, the 
exception being the mass regulator with r = 1/x. 
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oc l/{q^ + R+V" Y. For large n the contributions 
of xdxR^lk will dominate the x-integral in (8.30). Min¬ 
imising the absolute value of the integral then amounts to 
solving (8.32), so we still have to minimise r. Note also 
that this does not extremise the flow of all correlation 
functions f d'^x I 

It is also interesting to speculate about the most insta¬ 
ble regulator. It is found by maximising the integrand 
in (8.30) in the regularised momentum regime. This is 
achieved for rinstab = oo. If we also demand that r is 
monotone and that the gap (8.35) is saturated at some 
momentum, this leads to 


nnstab(a;) = - c) - I, (8.37) 

the sharp cut-off. Note that this argument concentrates 
on instability of the low momentum region of the flow. 
For high momenta maximal instability is obtained for the 
regulator i?cs = the mass cut-off. The related flow 
equation is an un-renormalised Callan-Symanzik equa¬ 
tion. Indeed, the results for critical exponents for scalar 
models in LPA are worse for the mass regulator [68] than 
that for the sharp cut-off. 

The stable and instable regulators (8.36) and (8.37) 
have been derived from (5.26) by dropping correlator- 
dependent terms. The regulators (8.36) and (8.37) can 
also be derived from (5.32c) in a very simple manner. In 
the present truncation (5.32c) has to be evaluated on L 2 
and boils down to 

__ ('o qo'j 

X + xrstab(x) + V" ~ X + xr±(x) + V" ' 

which can be converted into (8.34). This nicely shows 
the advantage of a simple functional criterion. 

Beyond LPA we are led to integrals as in (8.29) that 
also contain derivatives w.r.t. q. Then r also has to be 
differentiable to the given order. Such regulators exist, 
they are simply differentiable enhancements of (8.36). 


E. Optimisation in general truncation schemes 

In a general truncation and higher truncation order the 
correlation functions Ik resolve more structure of the flow 
operator AS' 2 . Roughly speaking, a solution to the func¬ 
tional optimisation criterion (5.26) minimises the expan¬ 
sion coefficients of AS '2 for a given truncation scheme. 
For example, in higher order derivative expansion the 
flow AS 2 /fc is projected on the part that contains higher 
order space-time derivatives. In momentum space and 


resorting to the representation (5.32) of the functional 
optimisation criterion (5.26), this amounts to differen¬ 
tiability of Giplp) w.r.t. momentum at the given order. 
Consequently the norm has to be taken in the space of 
differentiable functions with 




E 

|a|<n 


(n 


n! 

|Q;|)!ai!---ad! 


dp^i... dpo^d 


2 

L 2 


(8.39) 

where a G IN'^ and jaj = Eq. (8.39) defines the 

norm on Sobolev-spaces 77" with n G IN. Applied to 
the functional optimisation criterion, and leaving aside 
the intricacies discussed in section V D 3 we arrive at the 
following optimisation in nth order derivative expansion: 

||0A(G[((.o,Rstab])-0A(G[0o,O])L 


= min \\0x{G[(j)o, Rj_]) - 6 »a(G[^o, 0])L > (8-40) 

K_l 

for all A G R'*’. Here (f>o is either defined by the minimum 
of the potential or it maximises the propagator. 9\ has 
to meet the requirement of boundedness w.r.t the norm 
||.||„, as already discussed below (5.32b). This is achieved 
by using a nth-order differentiable version of (5.32b). We 
emphasise that the form of 9x is of no importance for the 
present purpose. The optimisation with (8.40) seems to 
depend on the full two-point function R = 0]. 

Now we proceed with the specific norm ||.||„ as indicated 
in section VD3 below (5.32a). The constraint (8.40) en¬ 
tails that the spectral values of G[(()o, i?stab] are as close 
as possible to that of the full propagator G[(j)o, Oj. More¬ 
over it entails maximal smoothness. Hence (5.32a) can 
be reformulated as 

|i0A(r(2)[0O,i?stab]+i?stab])L 

= inin ||0A(r(2)[0o, i?±] + R±])L , (8.41) 

H± 

for all A G R^. A solution of (8.41) provides a propaga¬ 
tor G[^ 0 j.Rstab] which is as close as possible to the full 
propagator G[(j)o, 0] as well as having minimal derivatives 
of order i <n. Eq. (8.41) also leads to the supplementary 
constraint for the stability criterion (5.10). The maximi¬ 
sation of the gap has to be supplemented by the minimi¬ 
sation of 

||r(2)[<^0,Rstab](p^) + Rstab(p^)]L , (8.42) 

within the class of i?stab singled out by (5.10). Here po is 
the momentum at which the propagator takes its maxi¬ 
mum. For an implementation of (8.42) see [71]. 
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In truncation schemes that carry a non-trivial momen¬ 
tum and field dependence [96, 97, 99, 100, 112, 125, 126, 
128, 129, 131-133], functional optimisation suggests the 
use of background field dependent regulators, or even reg¬ 
ulators with a non-trivial dependence on the full field. 
Evidently in the latter case structural truncations of the 
flows are inevitable, see also [47, 48]. In case momen¬ 
tum and field dependence are intertwined, as happens in 
the interesting truncation scheme put forward in [96, 97], 
functional optimisation directly implies the use of (back¬ 
ground) field dependent regulators. 

We continue with a brief discussion of a peculiar case 
relevant for the optimisation of QCD-fiows in Landau 
gauge QCD as initiated in [128, 129]. In case the spec¬ 
tral values A(p^) of the full propagator are not mono¬ 
tone in momentum, an optimised regulator does not re¬ 
solve the theory successively in momentum. This hap¬ 
pens for the gluon propagator in Landau gauge QCD 
[128, 131, 149, 151]. A propagator that is monotone in 
momentum violates the condition —dtG > 0 for some in¬ 
terval in t and some spectral values. This implies that 
the flow trajectory is not minimised for these spectral 
values. In turn, an optimised gluonic regulator can be 
constructed from 

~ {Z4,kls - T^^\p^))6[Z^kls - 
+ (r(")(p2) - - Z^kls\ , (8.43) 

where rg^^(p^) is the full two-point function at vanish¬ 
ing regulator, and a possibly smoothened step-function 9. 
Note that (8.43) boils down to the regulator (8.35) within 

LPA. The practical use of the suggestion (8.43) calls for 

( 2 ) 

an iterative solution of the flow about a suggestion Tq 
as described in section VIIIC. Apart from guaranteeing 
the mSTI, it also necessitates an appropriate choice of 
the renormalisation conditions. The latter ensures UV 
finiteness of such a flow. We also remark that within 
this approach further terms are required on the rhs of 
(8.43) in order to guarantee that the regulator vanishes 
if the cut-off scale tends to zero. The gluonic regulator 
(8.43) has to be accompanied with appropriate choices 
for ghost and quark regulators Rc and Rq respectively. 
A combined optimisation in {Ra, Rci Rq) may lead to 
a successive integrating out of fields as found already in 
the IR-optimisation in [128, 129] More details will be 
provided elsewhere [130]. 


The proof of an extremum being global is intricate 


In the light of the above results we add a further brief 
comment on the physical interpretation of optimisation 
as introduced in chapter V . The optimisation crite¬ 
rion is constructed from stability considerations. Sta¬ 
bility implies minimal integrated flows and hence quick¬ 
est convergence towards physics. At each order of the 
given truncation scheme the optimised propagators and 
correlation functions are as close as possible to the full 
propagator and correlation functions respectively. This 
minimises regulator artefacts, and triggers a most rapid 
approach towards the full theory. Moreover, optimised 
flows preserve the RG properties of the full theory within 
the regularisation as well as gradient flows, see (5.12). 
The above arguments emphasise the close structural re¬ 
lation of the optimisation criterion to the construction of 
both improved and perfect actions in lattice theory 
We emphasise that the optimisation can be implemented 
within an iterative procedure which leads to small addi¬ 
tional computational costs. 


IX. CONCLUSIONS 

The present work provides some structural results in 
the functional RG which may prove useful in further 
applications, in particular in gauge theories. We have 
derived flows (3.86) and their one-parameter reductions 
(3.28),(3.60) and (4.20) valid for a general class of correla¬ 
tion functions Ik defined in (3.14) with (3.51). This class 
of correlation functions Ik includes V-point functions as 
well as Dyson-Schwinger equations, symmetry relations 
such as Slavnov-Taylor identities, and flows in the pres¬ 
ence of composite operators, e.g. IV-particle irreducible 
flows. The present formulation also allows us to directly 
compute the evolution of observables in gauge theories. 
For example, the flows (3.60), (4.20) hold for the Wil¬ 
son loop and correlation functions of the Polyakov loop, 
see section VII A. This is a very promising approach 
to the direct computations of observables in the non- 
perturbative regime of QGD, e.g. the order parameter 
of the confinement-deconfinement phase transition. In 
section VIIIA we derived closed flows in the presence of 
general scale-dependent reparameterisations of the the¬ 
ory. This extends the options for scale-adapted param- 


An adaptation of the criterion (5.26) for lattice regularisations 
leads to improved actions and operators at lowest order of an 
expansion scheme based on the lattice spacing. 
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eterisations of the theory, and is particularly relevant in 
the context of rebosonisation. 

The functional framework developed here was used to 
systematically address the important issue of optimisa¬ 
tion, and to derive a functional optimisation criterion, see 
section VD3 (5.26),(5.29),(5.32). Optimal regulators are 
those, that lead to correlation functions as close as pos¬ 
sible to that in the full theory for a given effective cut-off 
scale. The criterion allows for a constructive use, and 
it is applicable to general truncation schemes. It can be 
also used for devising new optimised schemes, for exam¬ 
ples see section VIII, in particular section VIIID,VIIIE. 
The use of optimisation methods becomes crucial in more 
intricate physical problems such as the infrared sector of 
QCD, and can be used to resolve the pending problem of 
full UV-IR flows in QCD. 

Another important structural application concerns 
renormalisation schemes for general functional equations, 
e.g. DSEs and NPl effective actions. The functional flows 
(3.86) can be used for setting up of generalised BHPZ- 
type renormalisation schemes that are, by construction, 
consistent within general truncation schemes, see sec¬ 
tions VIA2,VIB3. Moreover, such subtraction schemes 
are very well adapted for numerical applications. 

The present setting also allows for a concise and flexi¬ 
ble representation of symmetry constraints, which is par¬ 
ticularly relevant in gauge theories. So far, the practical 
implementation of modified Slavnov-Taylor identities was 
restricted to their evaluation for specific momentum val¬ 
ues. The present setting allows for a functional imple¬ 
mentation that possibly adapts more of the symmetry, 
see section VIIB2,VIIB3. This opens a path towards 
improved truncation schemes in gauge theories relevant 
for a more quantitative computation in strongly interact¬ 
ing sectors of QCD. The above analysis also applies the 
Nielsen identities for gauge invariant flows of the geomet¬ 
rical effective action. 

In summary we have presented structural results that 
further our understanding of the Eunctional Renormali¬ 
sation Group. These results can be used to qualitatively 
and quantitatively improve ERG applications. 
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APPENDIX 

A. Metric 


This appendix deals with the non-trivial metric 7 in 
field space in the presence of fermions. The ultra-local 
metric 7 is diagonal in field space for scalars and gauge 
fields and is given by the e-tensor in fermionic space. For 
fa = (■*/') '4’)a the fermionic metric reads 


( 7 “*-) = 




(A.l) 


For raising and lowering indices we use the Northwest- 
Southeast convention, 


</>a = 4>^lba ■ (A.2) 


The metric has the properties 


76“ = 7““7bc = , 

7“6 = 7““7c6 = (-1)“'’<5.“ , (A.3) 


where 


(_l)a6 


— 1 , a and b fermionic 

1 otherwise. 


(A.4) 


Eq. (A.3) extends to indices a = and b = 

bi - ■ - bm with 


(-l)ab 


-1 


1 


a and b contain odd # of 
fermionic indices. 


(A.5) 


otherwise. 


For arbitrary vectors </>, <j) the properties (A.3) lead to 

(j)°-(j)a = = <l>b^°'l^a ■ (A. 6 ) 

Due to the Grassmann nature of the fermionic variables 
the order is important 

We close this appendix with an example. In general a 
(composite) field 4> consists of scalar components, gauge 
fields and fermions, the fundamental field reads in com¬ 
ponents 


{(jf) = (A.7) 
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The contraction of the fundamental (j) with itself leads to 

= «i'h7“Va = j d’^X (^^a„{x)ipa„{x) 


+AAx)AAx) 

1 (A. 8 ) 

where n labels the number of scalar fields, a the gauge 
group, and ^ sums over spinor indices and flavours. The 
current J related to (p is given by 

{da) - {J(pt JAt dlfjj dip) 


{d ) = {difiidA^dlpI dpp) 1 

(A.9) 

which implies schematically 


d°‘(j)a = {dip^ F JaA F dpjip F 'ijjdA- 

(A.IO) 

Moreover 


J>a = Ada = daAAb = <PbdA\ ■ 

(A.ll) 


B. Derivatives 


We deal with derivatives of functionals F[f] w.r.t. 
f{x) = (j){x) or f{x) = J{x). Derivatives are denoted 
as 


FAf] ■■= 


SF[f] 


FAf] ■■= 


dF[f] 


(B.l) 


<5/“ ’ Sfa 

that is, derivatives are always taken w.r.t. the argument 
of the functional F. Eq. (B.l) implies 

F,-[f]=A^Ft[(l>], Fa[f]=7abFAf], (B.2) 


which has to be compared with (A. 6 ). We also take 
derivatives w.r.t. some (logarithmic) scale s, e.g. s = t = 
In k. The total derivative of some functional F splits into 

= d,F[J]+d,rFa[J], 


C. Definition of AS„ 

The part of AS” that contains at least n > 1 derivatives 
w.r.t. the variable x, e.g. x[J] = JAi acting to the right, 
is given by 

^SAx,R] = /^S^,...^Ax.R]^ -(C.la) 

with coefficient 
AS'ai...a„[a;, R] 

= (c.ib) 

i>n ^ ^ 

The coefficients AS'ai -a„ are operators. The functionals 
(AS'a^.-.ai) the coefficients in a Taylor expansion of 
the operator AS' in powers of absorbing n deriva¬ 
tives w.r.t. X of AS[^ -I- (j) , 1?]. We emphasise that 
(ASai...a„ , R]) is a functional, it contains no deriva¬ 
tive operators. If interested in x = J, the expansion 
coefficients (AS'^^ [x , i?]) boil down to the Taylor co¬ 
efficients in an expansion of AS in They are the 
nth right derivatives of AS[x,i?] w.r.t. Xa, evaluated at 

X = A + 


D. Standard IPI flows 

For quadratic regulators (3.3) and a = a the flow (3.55) 
reads more explicitly 

dth + AGRG)bJA^ - [{dtr) (D.l) 

+ \{GRGATA^<^Ad - (l>bR^°^) GaJA = 0, 

where 


= dsF[<l,] + d.Aa FA<f\ , (B.3) 

as 

i.e., dsF[(j>] = ds\^F[(j>] and dsF[J] = ds\jF[J]. Partial 
derivatives w.r.t. the logarithmic infrared scale t = Infc 
we abbreviate with 

F = dtF. (B.4) 

General differential operators are similarly defined as 

DsF[J] = {^s+^A9^^g,+lJ\J'^A^)F[Jl 

DsF[4>] = {ds+-lg'jgidg^+-iAAb-^)F[(l>], (B.5) 

with partial derivatives according to (B.3). The defini¬ 
tions of this appendix directly carry over to the case of 
multi-indices a, b. 


{GRGA = GbaiF'^Gdc. 

For the derivation of (D.l) we have to express AS[^, R] 
in terms of derivatives w.r.t. (j) with the help of (3.50). 
For bosonic variables this is straightforwardly done. If 
fermionic variables are involved the ordering of terms be¬ 
comes important. We shall argue that 

Dab 5 6 

^ 5j^ TF 

_ ^ 5 TDab^ 8 

— ^ac8(t>c^ 

G ryabf 6 /^ \ 6 \ /^ i>ab/-i 5 5 /JA o^ 

8(l>d “T O-ca^ 54>d 54>c ' 

The only non-trivial term is the last one on the right 
hand side. Eq. (3.4) entails that for a being bosonic 
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(fermionic), b is bosonic (fermionic). If either a or c 
or both are bosonic we conclude Gac = G^a- Moreover 
either ^ and Gms^ or both are bosonic and (D.2) 
follows. If a,c both are fermionic, ^ and Gbdg^ are 
fermionic (as b is fermionic) and we have Gac = —Gca¬ 
ll follows that 

= {-^Gbd)-^ - Gbdj^-^ ■ (D.3) 

Inserting (D.3) into (D.2) the right hand side follows. We 
also conclude that for b, c fermionic 

= (D.4) 


The factor 7 ®/ originates in (3.46), Gac(rfc’°*' + i?^°) = 
j'^a- Inserting (D.4) into (D.2) we arrive at 

oab _6_6 _ ^ f^bcf-i 5 5 

^ 5J^ Jj^ ~ ^abJ^ ^cd S(i)^ 5(Pa 

-{GRG)adTk'’^^fl^fGgc^^. (D.5) 
with (G RG)ad = GabR^'^ Gcd- 
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